MATLAB programmet formler och delsresultat till Uppgift 4.
OVNODE(t,y) = [y(2) ,  - (1+t) y(2)^2 ]
%%%%%%%%%%%%%%%%%%%%
h = 0.1 ;   t = 0 ;  y = [1;0.5] ;
for i = 1:2  
t=t 
k1 = OVNODE(t,y) ;     k2 =  OVNODE(t+h,y+h*k1)  ;
['     y     ' , '   k1    '   ,   '      k2'  ] 
[ y  k1  k2 ] 
y = y + 0.5*h*(k1+k2) ;
t = t + h ;
end ;
t = t
y = y
%%%%%%   delsresultat :
t = 0 ;
k1 = [             1   
           - (1+0) ( 0.5)^2 ]

k2 =  [           1 + 0.1*1
                 - (1+0+0.1) ( 0.5+0.1*(-0.25))^2 ]
t =  0
                  y                              k1                              k2

            1.0000                   0.5                              0.4750
         0.5                       -0.2500                          -0.2482

t = 0.10000000000000

y =  [    1 +    0.5*0.1(0.5+0.4750)
          0.5 +   0.5*0.1(-0.2500-0.2482) ]

                   y                              k1                                k2
   1.0488                                0.4751                       0.4503
   0.4751                               -0.2483                       -0.2433

t = 0.2
y = [  1.0488  + 0.5*0.1* (0.4751 + 0.4503)
         0.4751 + 0.5*0.1 (-0.2483    - 0.2433) ]

        y 
   1.0950
   0.4505