... còn phương phápNewton‐Raphsonđòihỏitìmđạohàmdθ/dt.DovậycúngtasẽdùngphươngphápBrent.Tómlạithuậtoán giải bàitoángiátrịbiêngồmcácbước:‐Môtảgiátrịu1vàu2vâynghiệmucủa(3)‐Dùng phương phápBrenttìmnghiệmucủa(3).Chúýlàmỗibướclặpđòihỏitínhθ(u)bằngcách giải phương trình vi phân nhưlàbàitoánđiềukiệnđầu.‐Khiđãcóu, giải phương trình vi phân lầnnữađểtìmnghiệmTaxâydựnghàmbvp2shoot()để giải phương trình bậc2:function[t,x]=bvp2shoot(f,t0,tf,x0,xf,n,tol,kmax)%Giaiphuongtrinh:[x1,x2]ʹ=f(t,x1,x2)voix1(t0)=x0,x1(tf)=xfifnargin<8kmax=10;endifnargin<7tol=1e‐8;endifnargin<6n=100;enddx0(1)=(xf‐x0)/(tf‐t0);%chogiatridaucuaxʹ(t0)y0=[x0dx0(1)]ʹ;[t,x]=rungekutta(f,t0,tf,y0,n);%khoiganbgRK4 ... k3(4)k1 k k k 13 11ff3f fh132 6yyhf hf2h4 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 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 ... 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)=β...