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

function [y,X,Y]= Adams_13(X,Y,h,N)
    b=N-1; %4
    c=N-2; %3 
    d=N-3; %2
    a=N-4; %1
    for i=a:b
        q(i)=h*Y(i);
    end;
    for i=a:c
        dq(i)=q(i+1)-q(i);
    end;
    d2q(a)=dq(d)-dq(a);
    d2q(d)=dq(c)-dq(d);
    d3q(a)=dq(d)-d2q(a);
    dy=q(b)+1/2*dq(c)+5/12*d2q(d)+3/8*d3q(a);
    y=Y(b)+dy;
    X(N)=X(1)+h*(N-1);
    Y(N)=y;