Det skrå kast med luftmodstand

Det er fordi du løser en forkert differentialligning.

Du løser for. eks. for v_x ligningen v_x'(t) = k*v_x|v| medens det rigtige er v_x'(t) =  k*|v|²

Hvad foreslår du, jeg ændrer, jeg kan ikke lige se problemet? 

De angivede ligninger er altså rigtige nok...

Du har ret. Jeg har sovet i det

Men har du nogle forslag, til hvad er gøres forkert? :D

Ja, jeg har fundet din fejl. Giv mig 10-20 min.

Det tog lidt mere end 20 minutter. Jeg ville lige teste. Jeg kan nu kun mistænke din løsning i Maple for at være upræcis, fordi jeg får det samme med en fejlmargen. Jeg ved umiddelbart ikke hvordan man løser en 2.-ordens to-lignings Runge Kutta 4 i Maples dsolve, så jeg kan ikke tjekke vha. en funktion i Maple, så jeg måtte lave min egen. Jeg får nu de samme værdier som dig.

Det jeg har gjort er at indføre funktionen f, hvor

og så erstattet dine r-værdier med

Dermed kan implementeringen ske korrekt. Du havde derudover divideret med 2 i stedet for 6 for den fjerde indgang i opdateringen af r. Derudover kan du få min Maple-procedure til sammenligning.

Maple-implementering: https://www.dropbox.com/s/sasrinwyj622oik/RK4Studieportalen.mw?dl=0

Jeg er dog usikker på hvorfor dit Excel-dokument kræver flere punkter nu. Jeg har sikkert lavet en copy/paste-fejl i min implementering ud langs kolonnerne. r-værdierne passer dog pr. step, så den burde være god nok. Prøv fx at ændre h til 0.02.

Det er mærkeligt, at ved andre værdier så passer maple grafen og excel-grafen perfekt sammen. men din funktion ser mærkelig ud. hvorfor ganger du med h istedet for k i:
h*r3*(r3^2+r4^2)^0,5

Din metode var højst sandsynligt rigtig nok, jeg kunne bare ikke gennemskue om du brugte r-værdierne rigtigt i k_i-værdierne, så jeg tilføjede f-vektoren :)

Nå, jamen fantastisk! Sådan en Excel-implementering kan være svær at gennemskue, medmindre man selv har lavet det, så jeg kunne ikke umiddelbart overbevise mig om at den grafede det rigtigt, men værdierne passede jo på Maple-implementeringen. :)

Jeg er i øvrigt imponeret over at du har lavet det i Excel!

Må jeg godt kalde dig skat, endelig virker det perfekt !! :D

Haha, jeg er glad for at kunne være til hjælp. :)