blob: 7fc4029591b2416e537c8ede5500df13acc3a20e (
plain)
1
2
3
4
5
6
7
8
9
10
11
|
// example 6.13,
// caption: solve the IVP by euler midpoint method,
// with h=0.2,
// u'=f(t,u)
deff('[z]=f(t,u)','z=-2*t*u^2');
deff('[w]=fp(t,u)','w=-2*u^2-4*u*t');
[u] = eulermidpoint(1,0,1,.2,f,fp) // h=0.2;
|
