Matematik A matematik a
Kapitel Overblik Eulers metode
🔮 Matematik A

Eulers metode er en numerisk metode til at approksimere løsninger til differentialligninger. Når der ikke er en analytisk løsning, tager vi små skridt langs den approksimerede løsningskurve.

Du lærer at:
  • Forstå princippet bag Eulers metode: diskrete approksimationstrin
  • Anvende formlen yₙ₊₁ = yₙ + h·f(xₙ, yₙ) til at beregne approksimerede værdier
  • Vurdere nøjagtighed i forhold til skridtlængden h
  • Sammenligne numeriske approksimationer med analytiske løsninger
💡
Intuitionen

"Eulers metode er som at køre en bil i tåge: du kan ikke se langt frem, men du kan se hældningen præcis under dig og tage et lille skridt i den retning. Gentag dette mange gange og du finder vej frem!"

Eulers metode 🎯

Ikke alle differentialligninger kan løses med en pæn formel. Hvad gør vi, når separationsmetoden og andre analytiske teknikker ikke virker? Vi tager computeren (eller lommeregneren) til hjælp og beregner løsningen punkt for punkt — det er essensen af Eulers metode.

Ideen er genial i sin enkelhed: vi kender startpunktet og hældningen (fra differentialligningen), så vi “går” et lille skridt i hældningens retning. Gentag, og vi bygger løsningskurven op stykke for stykke — som at tegne en kurve med linealer!

Mestr Eulers metode og optjen +150 XP!


Teori: Grundideen bag Eulers metode

Betragt en differentialligning af 1. orden:

dydx=f(x,y),y(x0)=y0\frac{dy}{dx} = f(x, y), \quad y(x_0) = y_0

Vi kender altså:

  • Startpunktet (x0,y0)(x_0, y_0)
  • Hældningen i hvert punkt: f(x,y)f(x, y) — differentialligningen fortæller os, hvor stejl kurven er!

Den geniale indsigt: I et lille interval opfører kurven sig næsten som en ret linje. Så vi kan approksimere den næste yy-værdi ved at følge tangentlinjen et lille stykke.

Tangentlinjens ligning i (x0,y0)(x_0, y_0) med hældning f(x0,y0)f(x_0, y_0):

yy0+f(x0,y0)(xx0)y \approx y_0 + f(x_0, y_0) \cdot (x - x_0)

Hvis vi går et skridt af længde hh (dvs. x1=x0+hx_1 = x_0 + h), får vi:

y1y0+f(x0,y0)hy_1 \approx y_0 + f(x_0, y_0) \cdot h

Og ideen gentages: fra (x1,y1)(x_1, y_1) beregner vi (x2,y2)(x_2, y_2), osv.


Teori: Eulers formel

Eulers formel er:

yn+1=yn+hf(xn,yn)\boxed{y_{n+1} = y_n + h \cdot f(x_n, y_n)}

xn+1=xn+hx_{n+1} = x_n + h

hvor:

  • (xn,yn)(x_n, y_n) er den aktuelle position (det kendte punkt)
  • hh er skridtlængden (step size)
  • f(xn,yn)f(x_n, y_n) er hældningen i det aktuelle punkt (givet af differentialligningen)
  • (xn+1,yn+1)(x_{n+1}, y_{n+1}) er det nye approksimerede punkt

Algorithmen i ord:

  1. Start i (x0,y0)(x_0, y_0)
  2. Beregn hældningen: m=f(x0,y0)m = f(x_0, y_0)
  3. Tag et skridt: x1=x0+hx_1 = x_0 + h, y1=y0+hmy_1 = y_0 + h \cdot m
  4. Gentag fra det nye punkt (x1,y1)(x_1, y_1)
  5. Fortsæt til du når den ønskede xx-værdi

Intuition: Det er som at navigere i tåge. Du kan kun se retningen (hældningen) lige der, hvor du står. Du tager et lille skridt i den retning, tjekker retningen igen, tager et nyt skridt — og gradvist finder du frem!

Interaktiv Simulering

Eulers Metode & Retningsfelter

Undersøg hvordan Eulers numeriske metode følger retningsfeltets hældninger. Klik på grafen for at ændre startpunktet (x₀, y₀).

0011223344xy
Global fejl ved x = 4.0:
|y_Euler - y_Eksakt| = 0.0288

Kurven søger mod linjen y = x - 1. Dejlig stabiliserende effekt.

Parametre

Skridtlængde (h):0.5
0.1 (Meget præcis)1.0 (Store skridt)
Start x₀:0.00
Start y₀:1.00
Hvordan virker det?Eulers metode tager den lokale hældning y' = x - y i punktet (xn,yn)(x_n, y_n), og går en skridtlængde hh frem i den retning for at finde yn+1=yn+hcdotf(xn,yn)y_{n+1} = y_n + h \\cdot f(x_n, y_n).
Prøv dette: Sæt h=1,0h = 1{,}0 og derefter h=0,1h = 0{,}1 for at se, hvor hurtigt fejlen formindskes!

Iterativ tabelberegning (Hold musen over en række for at highlighte trinnet på grafen)

nnAktuelt punkt (xn,yn)(x_n, y_n)Hældning f(xn,yn)f(x_n, y_n)Trin-stigning hcdotfh \\cdot fNæste punkt yn+1=yn+hcdotfy_{n+1} = y_n + h \\cdot f
0(0.00, 1.0000)-1.00000.51.000=0.50000.5 \cdot -1.000 = -0.50000.5000
1(0.50, 0.5000)0.00000.50.000=0.00000.5 \cdot 0.000 = 0.00000.5000
2(1.00, 0.5000)0.50000.50.500=0.25000.5 \cdot 0.500 = 0.25000.7500
3(1.50, 0.7500)0.75000.50.750=0.37500.5 \cdot 0.750 = 0.37501.1250
4(2.00, 1.1250)0.87500.50.875=0.43750.5 \cdot 0.875 = 0.43751.5625
5(2.50, 1.5625)0.93750.50.938=0.46880.5 \cdot 0.938 = 0.46882.0313
6(3.00, 2.0313)0.96880.50.969=0.48440.5 \cdot 0.969 = 0.48442.5156
7(3.50, 2.5156)0.98440.50.984=0.49220.5 \cdot 0.984 = 0.49223.0078
8(4.00, 3.0078)0.99220.50.992=0.00000.5 \cdot 0.992 = 0.00003.0078

Vis Eksempel: Euler med y’ = 2y, tabelform ⚡

Opgave: Brug Eulers metode til at approksimere løsningen til y=2yy' = 2y, y(0)=1y(0) = 1 fra x=0x = 0 til x=1x = 1 med skridtlængde h=0,25h = 0{,}25.

Løsning:

Her er f(x,y)=2yf(x, y) = 2y (hældningen afhænger kun af yy).

Vi udfylder tabellen trin for trin:

Trin 0: (x0,y0)=(0,1)(x_0, y_0) = (0, 1)

f(x0,y0)=21=2y1=y0+hf(x0,y0)=1+0,252=1,5x1=0+0,25=0,25\begin{aligned} f(x_0, y_0) &= 2 \cdot 1 = 2 \\ y_1 &= y_0 + h \cdot f(x_0, y_0) = 1 + 0{,}25 \cdot 2 = 1{,}5 \\ x_1 &= 0 + 0{,}25 = 0{,}25 \end{aligned}

Trin 1: (x1,y1)=(0,25,  1,5)(x_1, y_1) = (0{,}25, \; 1{,}5)

f(x1,y1)=21,5=3y2=1,5+0,253=2,25x2=0,25+0,25=0,5\begin{aligned} f(x_1, y_1) &= 2 \cdot 1{,}5 = 3 \\ y_2 &= 1{,}5 + 0{,}25 \cdot 3 = 2{,}25 \\ x_2 &= 0{,}25 + 0{,}25 = 0{,}5 \end{aligned}

Trin 2: (x2,y2)=(0,5,  2,25)(x_2, y_2) = (0{,}5, \; 2{,}25)

f(x2,y2)=22,25=4,5y3=2,25+0,254,5=3,375x3=0,5+0,25=0,75\begin{aligned} f(x_2, y_2) &= 2 \cdot 2{,}25 = 4{,}5 \\ y_3 &= 2{,}25 + 0{,}25 \cdot 4{,}5 = 3{,}375 \\ x_3 &= 0{,}5 + 0{,}25 = 0{,}75 \end{aligned}

Trin 3: (x3,y3)=(0,75,  3,375)(x_3, y_3) = (0{,}75, \; 3{,}375)

f(x3,y3)=23,375=6,75y4=3,375+0,256,75=5,0625x4=0,75+0,25=1,0\begin{aligned} f(x_3, y_3) &= 2 \cdot 3{,}375 = 6{,}75 \\ y_4 &= 3{,}375 + 0{,}25 \cdot 6{,}75 = 5{,}0625 \\ x_4 &= 0{,}75 + 0{,}25 = 1{,}0 \end{aligned}

Tabeloversigt:

nnxnx_nyny_nf(xn,yn)=2ynf(x_n, y_n) = 2y_nhfh \cdot fyn+1y_{n+1}
00120,51,5
10,251,530,752,25
20,52,254,51,1253,375
30,753,3756,751,68755,0625

Eulers approksimation: y(1)5,0625y(1) \approx 5{,}0625.

Eksakt løsning: y=e21=e27,389y = e^{2 \cdot 1} = e^2 \approx 7{,}389.

Fejlen: 7,3895,063=2,326|7{,}389 - 5{,}063| = 2{,}326, svarende til en relativ fejl på ca. 31%31\%. Det er en stor fejl — men med mindre hh bliver det bedre! ✅


Teori: Beregning i tabelform

I praksis (og til eksamen) opstilles Eulers metode systematisk i en tabel. Her er skabelonen:

nnxnx_nyny_nf(xn,yn)f(x_n, y_n)hf(xn,yn)h \cdot f(x_n, y_n)yn+1=yn+hfy_{n+1} = y_n + h \cdot f
0x0x_0y0y_0beregnberegnberegn
1x0+hx_0 + hfra forrigeberegnberegnberegn

Tips til tabelberegning:

  1. Skriv alle kolonner op — det gør det let at tjekke fejl
  2. Afrund ikke undervejs — behold alle decimaler til slutresultatet
  3. Dobbelttjek f(xn,yn)f(x_n, y_n) — den mest almindelige fejl er at beregne hældningen forkert
  4. Antal skridt = xslutx0h\frac{x_{\text{slut}} - x_0}{h} — tjek at det giver et helt tal!

Vis Eksempel: Euler med f(x,y) der afhænger af x og y ⚡

Opgave: Brug Eulers metode til at approksimere y(1)y(1) for differentialligningen:

y=x+y,y(0)=1y' = x + y, \quad y(0) = 1

med skridtlængde h=0,5h = 0{,}5.

Løsning:

Her er f(x,y)=x+yf(x, y) = x + y — hældningen afhænger af både xx og yy!

Antal skridt: 100,5=2\frac{1 - 0}{0{,}5} = 2.

Trin 0: (x0,y0)=(0,1)(x_0, y_0) = (0, 1)

f(0,1)=0+1=1y1=1+0,51=1,5x1=0+0,5=0,5\begin{aligned} f(0, 1) &= 0 + 1 = 1 \\ y_1 &= 1 + 0{,}5 \cdot 1 = 1{,}5 \\ x_1 &= 0 + 0{,}5 = 0{,}5 \end{aligned}

Trin 1: (x1,y1)=(0,5,  1,5)(x_1, y_1) = (0{,}5, \; 1{,}5)

f(0,5,  1,5)=0,5+1,5=2y2=1,5+0,52=2,5x2=0,5+0,5=1,0\begin{aligned} f(0{,}5, \; 1{,}5) &= 0{,}5 + 1{,}5 = 2 \\ y_2 &= 1{,}5 + 0{,}5 \cdot 2 = 2{,}5 \\ x_2 &= 0{,}5 + 0{,}5 = 1{,}0 \end{aligned}

Tabel:

nnxnx_nyny_nf(xn,yn)=xn+ynf(x_n, y_n) = x_n + y_nhfh \cdot fyn+1y_{n+1}
00110,51,5
10,51,521,02,5

Eulers approksimation: y(1)2,5y(1) \approx 2{,}5.

Bemærk: Denne differentialligning er ikke separabel (x+yx + y kan ikke skrives som f(x)g(y)f(x) \cdot g(y)), så Eulers metode er her et af de bedste værktøjer! ✅


Vis Eksempel: Finer skridtlængde — h = 0,1 ⚡

Opgave: Løs y=2yy' = 2y, y(0)=1y(0) = 1 fra x=0x = 0 til x=0,5x = 0{,}5 med h=0,1h = 0{,}1 og sammenlign med den eksakte løsning.

Løsning:

Antal skridt: 0,50,1=5\frac{0{,}5}{0{,}1} = 5.

nnxnx_nyny_nf=2ynf = 2y_nhfh \cdot fyn+1y_{n+1}
001,00002,00000,20001,2000
10,11,20002,40000,24001,4400
20,21,44002,88000,28801,7280
30,31,72803,45600,34562,0736
40,42,07364,14720,41472,4883

Eulers approksimation: y(0,5)2,4883y(0{,}5) \approx 2{,}4883.

Eksakt løsning: y(0,5)=e20,5=e1=2,7183y(0{,}5) = e^{2 \cdot 0{,}5} = e^1 = 2{,}7183.

Fejl: 2,71832,4883=0,2300|2{,}7183 - 2{,}4883| = 0{,}2300, svarende til en relativ fejl på ca. 8,5%8{,}5\%.

Sammenlign med h=0,25h = 0{,}25: Med kun 2 skridt ville vi få y(0,5)=2,25y(0{,}5) = 2{,}25 (fejl 17%17\%). Halvering af hh halverer omtrent fejlen! ✅


Teori: Skridtlængde og fejlanalyse

Skridtlængden hh er den afgørende parameter i Eulers metode. Den bestemmer balancen mellem nøjagtighed og beregningsmængde.

Lokal afskæringsfejl:

I hvert enkelt trin begår Euler en fejl af størrelsen:

lokal fejl12h2y(xn)\text{lokal fejl} \approx \frac{1}{2} h^2 \cdot y''(x_n)

Denne fejl er proportional med h2h^2 — halverer vi hh, reduceres den lokale fejl med en faktor 4.

Global fejl:

Den samlede fejl efter N=xslutx0hN = \frac{x_{\text{slut}} - x_0}{h} skridt er:

global fejlh\text{global fejl} \propto h

Halverer vi hh, halveres den globale fejl (ca.), men vi skal tage dobbelt så mange skridt. Det er en afvejning!

Tommelfingerregler:

Skridtlængde hhNøjagtighedAntal skridtBeregningstid
Stor hhLavKort
Lille hhHøjMangeLang
Halver hh\approx dobbelt så godDobbeltDobbelt

Konvergens: Eulers metode konvergerer — dvs. når h0h \to 0, nærmer approksimationen sig den eksakte løsning. Metoden er en 1. ordens metode.


Vis Eksempel: Fejlanalyse — effekten af h ⚡

Opgave: For y=yy' = y, y(0)=1y(0) = 1 (eksakt løsning: y=exy = e^x), sammenlign Eulers approksimation af y(1)y(1) for forskellige hh-værdier.

Løsning:

Vi ved, at den eksakte værdi er y(1)=e12,7183y(1) = e^1 \approx 2{,}7183.

For y=yy' = y med Eulers metode bliver hvert trin:

yn+1=yn+hyn=yn(1+h)y_{n+1} = y_n + h \cdot y_n = y_n(1 + h)

Efter N=1/hN = 1/h trin:

yN=(1+h)1/hy_N = (1 + h)^{1/h}

hhAntal skridtyEuler(1)=(1+h)1/hy_{\text{Euler}}(1) = (1+h)^{1/h}FejlRelativ fejl
11(1+1)1=2,000(1+1)^1 = 2{,}0000,71826,4%
0,52(1,5)2=2,250(1{,}5)^2 = 2{,}2500,46817,2%
0,254(1,25)4=2,441(1{,}25)^4 = 2{,}4410,27710,2%
0,110(1,1)10=2,594(1{,}1)^{10} = 2{,}5940,1244,6%
0,01100(1,01)100=2,705(1{,}01)^{100} = 2{,}7050,0130,5%

Mønsteret er klart: halvering af hh halverer omtrent fejlen.

Og her er en sjov observation: udtrykket (1+h)1/h(1 + h)^{1/h} er faktisk den klassiske definition af ee for h0h \to 0! Eulers metode og Eulers tal ee hænger sammen. ✅


Teori: Hvornår bruges Eulers metode?

Eulers metode er et værktøj, du bruger, når:

  1. Differentialligningen ikke kan løses analytisk — dvs. separation og andre teknikker ikke virker (f.eks. y=x2+y2y' = x^2 + y^2)

  2. Du skal finde en numerisk tilnærmelse — til eksamen eller i praksis, hvor et tal er nok

  3. Du vil visualisere en løsningskurve — ved at plotte de beregnede punkter

  4. Computersimulering — Eulers metode er grundlaget for mange numeriske metoder i software

Fordele:

  • Simpel at forstå og implementere
  • Virker for enhver differentialligning af 1. orden
  • Kræver kun basale regneoperationer

Ulemper:

  • Relativt unøjagtig (1. ordens metode)
  • Kan give store fejl, hvis hh er for stor
  • Kræver mange skridt for god nøjagtighed

Bedre alternativer (uden for pensum i Matematik A):

  • Heuns metode (forbedret Euler)
  • Runge-Kutta-metoder (4. orden — bruges i praksis)

Vis Eksempel: Euler til en ikke-separabel ligning ⚡

Opgave: Brug Eulers metode til at approksimere y(2)y(2) for:

y=x2y,y(0)=0y' = x^2 - y, \quad y(0) = 0

med h=0,5h = 0{,}5.

Løsning:

Her er f(x,y)=x2yf(x, y) = x^2 - y. Ligningen er ikke separabel, da x2yx^2 - y ikke kan skrives som f(x)g(y)f(x) \cdot g(y).

Antal skridt: 200,5=4\frac{2 - 0}{0{,}5} = 4.

Trin 0: (x0,y0)=(0,0)(x_0, y_0) = (0, 0)

f(0,0)=020=0y1=0+0,50=0x1=0,5\begin{aligned} f(0, 0) &= 0^2 - 0 = 0 \\ y_1 &= 0 + 0{,}5 \cdot 0 = 0 \\ x_1 &= 0{,}5 \end{aligned}

Trin 1: (x1,y1)=(0,5,  0)(x_1, y_1) = (0{,}5, \; 0)

f(0,5,  0)=0,250=0,25y2=0+0,50,25=0,125x2=1,0\begin{aligned} f(0{,}5, \; 0) &= 0{,}25 - 0 = 0{,}25 \\ y_2 &= 0 + 0{,}5 \cdot 0{,}25 = 0{,}125 \\ x_2 &= 1{,}0 \end{aligned}

Trin 2: (x2,y2)=(1,0,  0,125)(x_2, y_2) = (1{,}0, \; 0{,}125)

f(1,0,  0,125)=10,125=0,875y3=0,125+0,50,875=0,5625x3=1,5\begin{aligned} f(1{,}0, \; 0{,}125) &= 1 - 0{,}125 = 0{,}875 \\ y_3 &= 0{,}125 + 0{,}5 \cdot 0{,}875 = 0{,}5625 \\ x_3 &= 1{,}5 \end{aligned}

Trin 3: (x3,y3)=(1,5,  0,5625)(x_3, y_3) = (1{,}5, \; 0{,}5625)

f(1,5,  0,5625)=2,250,5625=1,6875y4=0,5625+0,51,6875=1,40625x4=2,0\begin{aligned} f(1{,}5, \; 0{,}5625) &= 2{,}25 - 0{,}5625 = 1{,}6875 \\ y_4 &= 0{,}5625 + 0{,}5 \cdot 1{,}6875 = 1{,}40625 \\ x_4 &= 2{,}0 \end{aligned}

Tabel:

nnxnx_nyny_nf(xn,yn)=xn2ynf(x_n, y_n) = x_n^2 - y_nhfh \cdot fyn+1y_{n+1}
000000
10,500,250,1250,125
21,00,1250,8750,43750,5625
31,50,56251,68750,84381,4063

Eulers approksimation: y(2)1,406y(2) \approx 1{,}406.

Uden Eulers metode ville vi have svært ved at finde denne værdi, da ligningen ikke er separabel! ✅


Vis Eksempel: Geometrisk fortolkning ⚡

Opgave: Forklar geometrisk, hvorfor Eulers metode typisk underestimerer løsningen til y=yy' = y (en konveks kurve).

Løsning:

Betragt differentialligningen y=yy' = y med løsningen y=exy = e^x (en konveks kurve, dvs. y>0y'' > 0).

Geometrisk argument:

  1. I punktet (xn,yn)(x_n, y_n) beregner Euler hældningen f(xn,yn)=ynf(x_n, y_n) = y_n
  2. Euler følger tangentlinjen fra dette punkt et stykke hh fremad
  3. For en konveks kurve (y>0y'' > 0) ligger tangentlinjen under kurven
  4. Derfor giver Euler en yy-værdi, der er for lav

Det er som at køre langs en vej, der svinger opad. Hvis du kører ligeud (langs tangenten) i stedet for at følge vejen, ender du lavere end vejen!

For konkave kurver (y<0y'' < 0) er det omvendt: tangenten ligger over kurven, og Euler overestimerer.

Konklusion: Eulers systematiske fejl afhænger af kurvens krumning (yy''). Jo kraftigere krumning, jo større fejl — og jo vigtigere er det at bruge lille hh. ✅


Teori: Opsummering — Eulers metode trin for trin

Her er din komplette tjekliste til at anvende Eulers metode:

Forberedelse:

  1. Identificer f(x,y)f(x, y) fra differentialligningen y=f(x,y)y' = f(x, y)
  2. Aflæs startbetingelsen (x0,y0)(x_0, y_0)
  3. Vælg skridtlængden hh (eller brug den givne)
  4. Beregn antal skridt: N=xslutx0hN = \frac{x_{\text{slut}} - x_0}{h}

Beregning (gentag NN gange): 5. Beregn hældningen: mn=f(xn,yn)m_n = f(x_n, y_n) 6. Beregn næste punkt: yn+1=yn+hmny_{n+1} = y_n + h \cdot m_n og xn+1=xn+hx_{n+1} = x_n + h

Præsentation: 7. Opstil resultaterne i en tabel 8. Angiv den approksimerede slutværdi 9. Sammenlign med eksakt løsning, hvis mulig 10. Kommenter på fejlens størrelse og retning


🏋️ Træningsopgaver

Opgave 1: Brug Eulers metode til at approksimere y(1)y(1) for y=3yy' = 3y, y(0)=2y(0) = 2 med:

a) h=0,5h = 0{,}5 (2 skridt)

b) h=0,25h = 0{,}25 (4 skridt)

c) Sammenlign begge med den eksakte løsning y(1)=2e3y(1) = 2e^3.

Opgave 2: Løs y=xyy' = x - y, y(0)=2y(0) = 2 med Eulers metode fra x=0x = 0 til x=1x = 1 med h=0,25h = 0{,}25. Opstil tabel.

Opgave 3: En differentialligning y=y(1y)y' = y(1 - y) med y(0)=0,1y(0) = 0{,}1 beskriver logistisk vækst med M=1M = 1. Brug Eulers metode med h=0,5h = 0{,}5 til at approksimere y(2)y(2).

Opgave 4: Forklar med egne ord, hvorfor Eulers metode bliver mere nøjagtig, når hh gøres mindre.

Opgave 5: Beregn y(0,4)y(0{,}4) for y=2x+1y' = 2x + 1, y(0)=0y(0) = 0 med h=0,1h = 0{,}1. Sammenlign med den eksakte løsning (som du finder ved direkte integration).

Opgave 6 (Boss-kamp!): En kemisk reaktion beskrives af y=0,5y2y' = -0{,}5y^2, y(0)=4y(0) = 4. Denne ligning kan løses analytisk (prøv separation!), men brug Eulers metode med h=0,5h = 0{,}5 til at approksimere y(2)y(2) og sammenlign de to resultater.


Quiz – Test din forståelse

Matematik Boss-Kamp ⚔️

Løs opgavesættet

Op til +50 XP
Hvad er Eulers formel for det næste trin?