20 lines
483 B
Plaintext
20 lines
483 B
Plaintext
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
|