求助牛顿下山法解方程的根

cllc1985

x3-x-1=0 x0=0.6 精确到小数点后８位　　　哪位高人给个程序　谢谢啦

assist

找 数值计算方法（VC++版）里面有现成的代码

yelv123

原帖由 cllc1985 于 2007-12-24 14:46 发表

                               
登录/注册后可看大图

x3-x-1=0 x0=0.6 精确到小数点后８位　　　哪位高人给个程序　谢谢啦

这个是牛顿迭代的程序，自己再改一点就行了..
牛顿法求ln(x+sinx)=0的根

1. 
   
2. #include "stdio.h"
   
3.   #include "math.h"
   
4.   int newt(x,eps,js,f)
   
5.   int js;
   
6.   double *x,eps;
   
7.   void (*f)();
   
8.   { int k,l;
   
9.     double y[2],d,p,x0,x1;
   
10.     l=js; k=1; x0=*x;
    
11.     (*f)(x0,y);
    
12.     d=eps+1.0;
    
13.     while ((d>=eps)&&(l!=0))
    
14.       { if (fabs(y[1])+1.0==1.0)
    
15.           { printf("err\n"); return(-1);}
    
16.         x1=x0-y[0]/y[1];   /*迭代计算  x1＝x0－f(x0)/f'(x0)*/
    
17.         (*f)(x1,y);     /*计算 f(x1)与f'(x1)*/
    
18.         d=fabs(x1-x0); p=fabs(y[0]);
    
19.         if (p>d) d=p;
    
20.         x0=x1; l=l-1;
    
21.       }
    
22.     *x=x1;
    
23.     k=js-l;
    
24.     return(k);
    
25.   }
    

复制代码

主函数：：
参数说明：
x
双精度实型变量指针。在该指针指向的变量中存放迭代初值；返回时在该指针指向的变量中存放终值。
eps
双精度实型变量。控制精度要求。
js

整型变量。最大迭代次数
f

无类型函数指针变量。指向计算方程左端函数 f(x)值及其导数值f'(x )的函数。
该函数由用户自编，其形式为

void f(x，y)


double x，y[2]；


{ y[0]＝f(x)的表达式；


y[1]＝f'(x)的表达式；



return；}

本函数返回一个整型标志值。若返回的标志值小于0，则表示在出现了 的情况，函数中将输出信息“err”；若返回的标志值等于最大迭代次数，则表示还未满足精度要求，此时返回的实根终值只作为参考；若返回的标志值大于或等于0且小于最大迭代次数，则表示正常返回。

1.   #include "math.h"
   
2.   #include "stdio.h"
   
3.   #include "dnewt.c"
   
4.   main()
   
5.   { int js,k;
   
6.     double x,eps;
   
7.     eps=0.000001; js=60;
   
8.     x=0.15;/*迭代初值*/
   
9.     k=dnewt(&x,eps,js);
   
10.     if (k>=0)
    
11.       printf("k=%d  x=%13.7e\n",k,x);
    
12.     printf("\n");
    
13.   }
    
14. 
    
15.   void newtf(x,y)
    
16.   double x,y[2];
    
17.   { y[0]=log(x+sin(x));
    
18.     y[1]=(1+cos(x))/(x+sin(x));
    
19.     return;
    
20.   }

复制代码

风花雪月

其实这些问题自己稍微花点时间就能找到

dgl0611

计算方法很多，随便一本数值分析的书上就有源程序