... ʹ(x)=f(x,y(x))hayy’=f(x)ifn<2n=2;endh=(xf‐xo)/n;X=zeros(n+1,1);M=max(size(yo));%sophuongtrinh(socotcuamatranY)Y=zeros(n+1,M);%datdieukiendaux=xo;X(1)=x;y=yo;Y(1,:)=yʹ;fori=1:nifnargin(fxy)>1k1=h*feval(fxy,x,y);elsek1=h*feval(fxy,x);endy=y+k1;x=x+h;X(i+1)=x;Y(i+1,:)=yʹ;endfunctiondy=f1(t,y)dy=zeros(3,1);dy(1)=y(2)*y(3);dy(2)=‐y(1)*y(3);dy(3)=‐0.51*y(1)*y(2);Đểgiải phương trình chobởihàmf1(x,y)tadùng chương trình cteuler.m:clearall,clca=0; 360CHƯƠNG 7: CÁC PHƯƠNG TRÌNH VI PHÂN THƯỜNG §1.BÀITOÁNCAUCHY Một phương trình vi ... 361Nếu phương trình vi phân cóbậccaohơn(n),nghiệmsẽphụthuộcvàonhằngsốtuỳý.Đểnhậnđượcmộtnghiệmriêng,taphảichonđiềukiệnđầu.Bàitoánsẽcógiátrịđầunếuvớigiátrịxođãchotachoy(xo),y′(xo),y″(xo), Một phương trình vi phân bậcncóthểđưavềthànhmộthệ phương trình vi phân cấp1.Vídụnếutacó phương trình vi phân cấp2:⎩⎨⎧β=′α=′=′′)a(y,)a(y)y,y,x(fyKhiđặtu=yvàv=y′tanhậnđượchệ phương trình vi phân cấp1:⎩⎨⎧=′=′)v,u,x(gvvuvớiđiềukiệnđầu:u(a)=αvàv(a)=β ... chỉnhbằng phương phápngoạisuyRichardson.Khuyếtđiểmnàyđượckhácphụctrong phương phápBurlisch‐Stör.Ýtưởngcủa phương phápnàylàápdụng phương phápđiểmgiữatrêntừngđoạn.Taxâydựnghàmburlischstoer()đểthựchiệnthuậttoánnày:function[xout,yout]=burlischstoer(f,x,x1,y,H,tol)%PhuongphapBulirsch‐Stoergiaiphuongtrinhyʹ=F(x,y)hayy’=f(x).%[x,x1]lakhoangtimnghiem.%H=dotangsaumoilantinhifsize(y,1)>1y=yʹ;end%yphailavectohangifnargin<6tol=1.0e‐8;endn=length(y);xout=zeros(2,1);yout=zeros(2,n);xout(1)=x;yout(1,:)=y;k=1;whilex<x1 3695 1−115452− 70 27 35 27 ∗0 277 143366 7 8163155296 175 512 575 1382444 275 1105922534096512 177 114Saisốsẽlà: E(h)=y5(x+h)‐y4(x+h)==−∑6iiii1(C...