контрольные работы 7 вариант / 14 / Eiler_14

function [T,P1,P2,P3,P4]=Eiler_14(t0,t1,p10,p20,p30,p40,N)
dt=(t1-t0)/N;
t(1)=t0;
p1(1)=p10;
p2(1)=p20;
p3(1)=p30;
p4(1)=p40;
for i=1:N
t(i+1)=t(1)+dt*i;
p1(i+1)=p1(i)+dt*F1_14(t(i),p1(i),p2(i),p3(i),p4(i));
p2(i+1)=p2(i)+dt*F2_14(t(i),p1(i),p2(i),p3(i),p4(i));
p3(i+1)=p3(i)+dt*F3_14(t(i),p1(i),p2(i),p3(i),p4(i));
p4(i+1)=p4(i)+dt*F4_14(t(i),p1(i),p2(i),p3(i),p4(i));
end;
T=t;
P1=p1;
P2=p2;
P3=p3;
P4=p4;

function z=F1_14(t,p1,p2,p3,p4)
z=-0.7*p1+2.3*p2;

function z=F2_14(t,p1,p2,p3,p4)
z=0.7*p1+1.8*p2+4.6*p3;

function z=F3_14(t,p1,p2,p3,p4)
z=0.7*p2+1.8*p3+6.9*p4;

function z=F4_14(t,p1,p2,p3,p4)
z=0.7*p3-6.9*p4;