... ʹ(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 ... Một phương trình viphânbậcncóthểđưavềthànhmộthệ phương trình viphâncấp1.Vídụnếutacó phương trình viphâncấp 2: ⎩⎨⎧β=′α=′=′′)a(y,)a(y)y,y,x(fyKhiđặtu=yvàv=y′tanhậnđượchệ phương trình viphâncấp1:⎩⎨⎧=′=′)v,u,x(gvvuvớiđiềukiệnđầu:u(a)=αvàv(a)=β Các phương pháp giải phương trình vi phânđược trình bày trong chương nàylà các phương pháprờirạc:đoạn[a,b]đượcchiathànhnđo ... ương trình đầutiên trong (2)dùngcôngthứcEulerđểthaycho phương phápđiểmgiữa. Các phương trình kháclà các côngthứcđiểmgiữa.Kếtquảcuốicùnglàtrungbìnhcộngcủayn trong (2)vàtacó:[]−+=...