... k3(4)k1 k k k13 11ff3f fh132 6 yyhf hf2h4 360 CHƯƠNG 7: CÁC PHƯƠNG TRÌNH VI PHÂN THƯỜNG §1.BÀITOÁNCAUCHY Một phương trình vi phân cấp 1 có thể vi t dưới dạng giảiđược′=yf(x,y)màtacóthểtìmđượchàmytừđạohàmcủanó.Tồntạivôsốnghiệmthoảmãn phương trình trên.Mỗinghiệm ... 384ifnargin(f)>1F(i‐1,:)=feval(f,t(i),y(i,:));elseF(i‐1,:)=feval(f,t(i));endendfori=4:np=y(i‐3,:)+(4*h/3)*(2*F(1,:)‐F(2,:)+2*F(3,:));%Pt.(1)ifnargin(f)>1F(4,:)=f(t(i+1),p);%Pt.(2)elseF(4,:)=f(t(i+1));endy(i+1,:)=y(i‐1,:)+(h/3)*(F(2,:)+4*F(3,:)+F(4,:));%Pt.(3)F(1,:)=F(2,:);F(2,:)=F(3,:);ifnargin(f)>1F(3,:)=f(t(i+1),y(i+1,:));elseF(3,:)=f(t(i+1));endendĐểgiải phương trình tadùng chương trình ctmilne.m:clearall,clca=0;b=1;y=@f2;ya=1;n=10;[t,y]=milne(y,a,b,ya,n);plot(t,y)§10.BÀITOÁNGIÁTRỊBIÊN1.Kháiniệmchung:Taxét bài toán tìmnghiệmcủa phương trình: 385′′ ′==α=βyf(x,y,y)y(a) ,y(b)Đâylà bài toán tìmnghiệmcủa phương trình vi phân khi biếtđiều kiện biên vàđược gọi là bài toán giá trị biênhaiđiểm.Trong bài toán giá trị đầu tacó ... sốphương trình phituyếnnếu phương trình vi phân làphituyến.Các phương trình nàyđượcgiảibằng phương pháplặpnênrấttốnthờigiantính toán. Vìvạy vi cgiải bài toán biênphituyếnrấtkhó. Ngoàira, đối với phương pháplặp, vi c chọn giá trị đầu ...