контрольные работы 7 вариант / 13 / RungeKutt_13

function [X,Y]= RungeKutt_13(y0,x0,x1,h)
    N=(x1-x0)/h;
    x(1)=x0;
    y(1)=y0;
    for i=2:N+1
        x(i)=x(1)+h*(i-1);
        k1=h*F(x(i-1),y(i-1));
        k2=h*F(x(i-1)+h/2,y(i-1)+k1/2);
        k3=h*F(x(i-1)+h/2,y(i-1)+k2/2);
        k4=h*F(x(i-1)+h,y(i-1)+k3);
        y(i)=y(i-1)+1/6*(k1+2*k2+2*k3+k4);
    end;
    X=x;
    Y=y;
    
function z=F(x,y)
    z=x+cos(y/exp(1));