Скачиваний:
13
Добавлен:
23.11.2018
Размер:
372 б
Скачать
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;
Соседние файлы в папке 13