контрольные работы 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;