Runge-Kutta method: Difference between revisions
no edit summary
SqrtNegInf (talk | contribs) m (→{{header|Perl 6}}: replace zip with hyper in loop control) |
No edit summary |
||
Line 1,632:
10 675.99994901671 5.09832902935159E-05
</pre>
=={{header|Pure Basic}}==
{{trans|BBC Basic}}
<lang PureBasic>EnableExplicit
Define.i i
Define.d y=1.0, k1=0.0, k2=0.0, k3=0.0, k4=0.0, t=0.0
If OpenConsole()
For i=0 To 100
t=i/10
If Not i%10
PrintN("y("+RSet(StrF(t,0),2," ")+") ="+RSet(StrF(y,4),9," ")+#TAB$+"Error ="+RSet(StrF(Pow(Pow(t,2)+4,2)/16-y,10),14," "))
EndIf
k1=t*Sqr(y)
k2=(t+0.05)*Sqr(y+0.05*k1)
k3=(t+0.05)*Sqr(y+0.05*k2)
k4=(t+0.10)*Sqr(y+0.10*k3)
y+0.1*(k1+2*(k2+k3)+k4)/6
Next
Print("Press return to exit...") : Input()
EndIf
End</lang>
{{out}}
<pre>y( 0) = 1.0000 Error = 0.0000000000
y( 1) = 1.5625 Error = 0.0000001457
y( 2) = 4.0000 Error = 0.0000009195
y( 3) = 10.5625 Error = 0.0000029096
y( 4) = 25.0000 Error = 0.0000062349
y( 5) = 52.5625 Error = 0.0000108197
y( 6) = 100.0000 Error = 0.0000165946
y( 7) = 175.5625 Error = 0.0000235177
y( 8) = 289.0000 Error = 0.0000315652
y( 9) = 451.5625 Error = 0.0000407232
y(10) = 675.9999 Error = 0.0000509833
Press return to exit...</pre>
=={{header|Python}}==
| |||