Phương pháp số giải phương trình phi tuyến sẽ giúp các bạn có thêm một số kiến thức tổng hợp về cách viết lệnh, hàm, các thuật toán để giải phương trình phi tuyến tốt hơn. | program giai_phuong_trinh_phi_tuyen; uses crt; var x0,x,eps,ss,fx,dfx: real; i:integer; {----------------------------------------------------------} { viet chuong trinh con tinh ham F(x) } Function F(x: real) : real; var a: real; begin a:= Ln(x)/Ln(10) - Cos(x); F:=a; end; {----------------------------------------------------------} { viet chuong trinh con tinh ham dF(x) } Function dF(x: real) : real; var b: real; begin b:= (1/x)*(1/Ln(10))+Sin(x); dF:=b; end; {----------------------------------------------------------} { viet thuc tuc thuc: hien chuong trinh NEWTON } Procedure NEWTON; { khoi kiem tra dieu kien cua EPSILON } begin repeat write(' Nhap vao Epsilon = '); readln(eps); if (eps1) then begin write(' Nhap sai, hay nhap lai'); writeln; end; until (eps>0) and (eps eps then begin writeln('-------------------------'); writeln('LanLap:',i); writeln('X=',x:8:6); writeln('Fx=',F(x):8:6); writeln('dFx=',dF(x):8:6); end; if ss<= eps then break; end; end; {---------------------------------------------------------} { than chuong trinh chinh} BEGIN CLRSCR; writeln(' Giai phuong trinh phi tuyen '); writeln('-------------------------------------------'); writeln; newton; writeln; writeln(' NGHIEM X = ',x:8:6); write('ket thuc chuong trinh '); writeln; write(' -------- Copyright Nguyen Thanh Que - KTHH2 - K55 --------'); readln; END.
