Paraang Euler

Sa matematika at komputasyonal na agham, ang paraang Euler ay isang unang-order na paraang numerikal para sa paglutas ng mga ekwasyong ordinaryong diperensiyal(ODE) na may isang ibinigay na inisyal na halaga. Ito ang pinaka basikong paraang hayagan para sa numerikal na integrasyon ng mga ekwasyong ordinaryong diperensiyal at ang pinaka simpleng paraang Runge-Kutta. Ang paraang Euler ay ipinagangalan kay Leonard Euler na trumato nito sa kanyang aklat na Institutionum calculi integralis (na inilimbag noong 1768–70). Ang paraang Euler ay isang paraang unang order na nangangahulugan ang lokal na pagkakamali(pagkakamali kada hakbang) ay proporsiyonal sa kwadrado ng sukat ng hakbang at ang pagkakamaling global(pagkakamali sa isang ibinigay na panahon) ay proporsiyonal na sukat ng hakbang. Ito ay dumaranas rin mula sa mga problemang pagiging matatag. Para sa mga dahilang ito, ang paraang Euler ay hindi kadalasang ginagamit sa pagsasanay. Ito ay nagsisilbing basehan upang lumikha ng mas komplikadong mga paraan.

Inpormal na paglalarawang heometrikal
Isaalang alang ang problema ng pagkukwenta ng hugis ng isang hindi alam na kurba na nagsisimula sa isang ibinigay na punto at sumasapat sa isang ibinigay na ekwasyong diperensiyal. Dito, ang isang ekwasyong diperensiyal ay maaaring akalain bilang isang pormula kung saan ang lihis ng linyang tangent sa kurba ay maaaring kwentahin sa anumang punto sa kurba kapag ang posisyon ng puntong ito ay nakwenta na. Ang ideya ay bagaman ang kurba ay hindi alam sa simula, ang pasimulang punto nito na ating tutukuying  ay alam. Pagkatapos, mula sa ekwasyong diperensiyal, ang lihis sa kurba sa  ay maaaring kwentahin at kaya ang linyang tangent. Kumuha ng isang maliit na hakbang sa kahabaan ng linyang tangent hanggang sa isang puntong  Sa kahabaan ng maliit na hakbang na ito, ang lihis ay hindi labis na nagbabago kaya ang  ay magiging malapit sa kurba. Kung tayo'y magpapanggap na ang  ay nasa kurba pa rin, ang parehong pangangatwiran gaya ng sa puntong  sa itaas ay maaring gamitin. Pagkatapos ng ilang mga hakbang, ang isang kurbang poligonal na  ay kukwentahin. Sa pangkalahatan, ang kurbang ito ay hindi lumilihis ng labis na malayo mula sa orihinal na hindi alam na kurba at ang pagkakamali sa pagitan ng dalawang mga kurba ay maaaring gawing maliit kung ang sukat ay sapat na maliit ang ang interbal ng komputasyon ay may hangganan.

Pormulasyon ng paraan
Ipagpalagay na nais nating tantiyahin ang solusyo ng problemang inisyal na halaga na 

Pumili ng isang halagang  para sa sukat ng bawat hakbang at itakda ang  . Ngayon, ang isang hakaban ng paraang Euler mula sa   to  ay 

ng halaga ng  ay isang aproksimasyon ng solusyon sa ODE sa panahong  : . Ang paraang Euler ay hayagan, i.e. ang solusyong  ay isang hayagang punsiyon ng  for .

Bagaman ang paraang Euler ay nag-iintegrado ng isang unang order na ODE, ang anumang ODE ng order na N ay maaaring ikatawan bilang isang unang order na ODE: upang tratuhin ang ekwasyong
,

ating ipakikilala ang mga sumusuportang baribulong  at makamit ang katumbas na ekwasyong 

Ito ang isang sistemang unang order sa bariabulong  at maaaring hawakan ng paraang Euler o sa katotohanan ng anumang ibang skema para sa mga sistemang unang order.

Halimbawa

Sa ibinigay na problemang inisyal na halagang 

nais nating gamitin ang paraang Euler upang tantiyahin ang .

Paggamit ng sukat ng hakbang na katumbas ng 1

Ang paraang Euler ay 

kaya dapat nating kwentahin ang  . Sa simpleng ekwasyong diperensiyal na ito, ang punsiyong  ay inilalarawan ng . Mayroon tayong

Sa paggawa ng hakbang sa itaas, ating natagpuan ang lihis ng linya na tangent sa kurbang solusyon sa puntong . Tandaan na ang lihis ay inilalarawan bilang pagbabago sa  na hinati ng pagbabago sa , o .  

Ang susunod na hakban ay paramihin ang halaga sa itaas sa sukat ng hakbang na  na ating kukuning katumbas sa isang narito: 

Dahil ang sukat ng hakbang ang pagbabago sa , kung ating pararamihin ang sukat ng hakban at ang lihis ng tangent, ating makukuha ang pagbabago sa  halagang  value. Ang halagang ito ay idadagdag naman sa inisyal na halagang  upang makamit ang susunod na halagang gagamit para sa mga pagkukwenta. 

Ang mga hakbang sa itaas ay dapt ulitin upang mahanap ang ,  at .

Dahil sa paulit ulit na kalikasan ng algoritmong ito, makakatulong na ayusin ang mga pagkukwenta sa isang anyong tsart gaya ng makikita sa ibaba upang maiwasan ang paggawa ng mga pagkakamali. 

{| class="wikitable"
|-
!  !!  !!  !! !!  !!  !! 
|-
| 0 || 1 || 0 || 1 || 1 || 1 || 2
|-
| 1 || 2 || 1 || 2 || 1 || 2 || 4
|-
| 2 || 4 || 2 || 4 || 1 || 4 || 8
|-
| 3 || 8 || 3 || 8 || 1 || 8 || 16
|}

Ang konklusyon ng pagkukwentang ito ang  . Ang eksaktong solusyon ng ekwasyong diperensiyal ay  kaya ang . Kaya ang aproksimasyon ng paraang Euler ay hindi napaka buti sa kasong ito. Gayunpaman, gaya ng ipinapakita ng pigura, ang pag-aasal nito ay kwalitatibong tama.

Paggamit ng ibang mga sukat ng hakbang 

Gaya ng iminungkahi sa introduksiyon, ang paraang Euler ay mas tiyak kung ang mga sukat ng hakbang na  ay mas maliit. Ang tabla sa ibaba ay nagpapakita ng resultang may iba't ibang mga sukat ng hakbang. Ang taas na row ay tumutugon sa halimbawa ng nakaraang seksiyon at ang ikalawang row ay ipinapakita sa pigura.

{| class="wikitable"
|-
! sukat ng hakbang!! resulta ng paraang Euler !! pagkakamali
|-
| 1 || 16 || 38.598
|-
| 0.25 || 35.53 || 19.07
|-
| 0.1 || 45.26 || 9.34
|-
| 0.05 || 49.56 || 5.04
|-
| 0.025 || 51.98 || 2.62
|-
| 0.0125 || 53.26 || 1.34
|}

Ang itinalang pagkakamali sa huling column ng tabla ang pagkakaiba sa pagitan ng eksaktong solusyon sa   at ang aproksimasyong Euler. Sa ilalim ng tabla, ang sukat ng hakban ay kalahati ng sukat ng hakban sa nakaraang row at ang pagkakamali ay tinatayang kalahati rin ng pagkakamali ng nakaraang row. Ito ay nagmumungkahing ang pagkakamali ay mga proporsiyonal sa sukat ng hakban, kahit papano para sa katamtamang maliit na mga halaga ng sukat ng hakbang. Ito ay totoo sa pangkalahatanat para sa ibang mga ekwasyon. Ang ibang mga paraan gaya ng paraang gitnang punto na ipinapakita rin sa mga pigura ay umaasal ng mas mapapaboran. Ang pagkakamali ng paraang gitnang punto ay mga proporsional sa kwadrado ng bawat sukat ng hakbang. Dahil dito, ang paraang Euler ay sinasabing isang paraang unang order samantalang ang paraang gitnang punto ay ikalawang order. Pwede nating iekstrapola mula sa taas na tabla na ang bawat sukat ng hakban na kailangan upang makuha ang sagot na tama sa mga lugar na tatlong desimal ay tinatayang 0.00001 na ngangahulugan kailangan nating mga 400,000 hakbang. Ang malaking bilang na ito ng mga hakbang ay nag aatas ng isang mataas na gastos na komputasyon. Dahil dito, ang mga tao ay karaniwang gumagamit ng alternatibong mga paraang mataaas na order gaya ng mga paraang Runge-Kutta o mga paraang linyar maraming hakbang lalo na kung ang isang mataas na pagiging tumpak ay ninanais.

Mga sanggunian

Mga sanggunian
 .
 .
 .
 .
 
 .
 .

Mga ekwasyong numerikal diperensiyal
Mga paraang Runge-Kutta
Mga paraang unang order