#include #include double g(double alpha,double x){ return tan(x)-alpha*exp(x); } double dgdx(double alpha,double x){ double sec=1/cos(x); return sec*sec-alpha*exp(x); } double newton(double alpha,double x){ return x-g(alpha,x)/dgdx(alpha,x); } double f(double x,double alpha){ for(int i=1;i<10;i++){ x=newton(alpha,x); } return x; } int main(){ double x=0.8766; for(int i=0;i<16;i++){ double alpha=0.5+0.1*i; x=f(x,alpha); printf("f(%.17g)=%.15g\n",alpha,x); } return 0; }