计算方法实验三上机程序.docxVIP

实验题二 2.1简单迭代法： # includestdio.h #in cludemath.h #define xO 3.0 #define MAXREPT 1000 #define EPS IE-6 #define G(x) pow(12*x+sin(x)-lzl-0/3) void main() { int i; double x_k=xO,x_kl=xO; printf(k xk\n); tor(i=0;iMAXREPT;i++) { printf(H%d %g\n,,/i/x_kl); x_kl=G(x_k); if (fabs(x_kl-x_k)EPS) { printffTHE ROOT IS x=%g,k=%d\n,,/x_kl,i); return; } x_k=x_kl; } printf(AFTER %d repeate,no solved.\n,MAXREPT); } -d (4 3?  2.2牛顿迭代法： #in cludestdio.h #in cludemath.h #define xO 3.0 #define MAXREPT 1000 #define EPS IE-6 //define G(x) x-(x*x*x-sin(x)-12*x+l)/(x*x*3-cos(x)-12) void main() { int i; double x_k=xO,x_kl=xO; printfCk xk\nH); for(i=0;iMAXREPT;i++) { printf(u%d %g\n,i/x_kl); x_kl=G(x_k); if (fabs(x_kl-x_k)EPS) { ~ printf(HTHE ROOT IS x=%g/k=%d\n,x_kl,i); return; } x k=x kl; } printf(HAFTER %d repeate,no solved.\MAXREPT); } S M .. ■ b * 2 tfow iRWft BuM looH tH| AN 丄心口 4— C.- ▲ □..- -ox i 3 3.41401 3.41013 3.41013 W ROOT IS s?14IOUkH *r?” any key to continue im.eMe ? erri I ft Cai ? iT EEK Bwaldxltbut \ FX ia 皿―i FW ut Fdt ： \ SQL / 2.3埃特金迭代： #in cludestdio.h #in cludemath.h ttdefine xO 3.0 //define MAXREPT 1000 //define EPS IE-6 #define G(x) (pow(x,3)-sin(x)+l)/12.0) void main() { int i; double x_k=xO,x_kl=xO; double y,z; printf(nk xk\nH); for(i=0;iMAXREPT;i++) { printf(%d %g\nHJ,x_kl); y=G(x_k); z=G(y); x_kl=z-((z-y)*(z-y))/(z-2*y+x_k); if (fabs(x_kl-x_k)EPS) { ~ printf(HTHE ROOT IS x=%g Iterations times k=%d\n,,/x_kl/i); return; } x k=x kl; } printf(HAFTER %d repeate,no solved.\MAXREPT); A Q0 细?U ■ 口蛋害?「 d w 3.1牛顿迭代法 #includestdio.h #in cludemath.h #define xO 1.5 #define MAXREPT 1000 #define EPS IE-6 #define G(x) x-(pow(xz3)+pow(xz2)-3*x-3)/(3*pow(x,2)+2*x-3) void main() { int i; double x_k=xO,x_kl=xO; printf(k xk\n); for(i=0;iMAXREPT;i++) { printf(u%d %g\n,,,i/x_kl); x_kl=G(x_k); if (fabs(x_kl-x_k)EPS) { printf(HTHE ROOT IS x=%g/k=%d\n,x_kl,i); return; } x_k=x_kl; } ~ ~ printf(HAFTER %d repeate,no solved.\n,MAXREPT); 3?2蹲器 snc-udecsidio.hv #inc-ude〈math.hv 芸efine MAXRE