| 网站首页 | 业界新闻 | 小组 | 威客 | 人才 | 下载频道 | 博客 | 代码贴 | 在线编程 | 编程论坛
欢迎加入我们,一同切磋技术
用户名:   
 
密 码:  
共有 512 人关注过本帖
标题:如何对非线形方程的二分法求根进行改进?
只看楼主 加入收藏
暖暖BOY
Rank: 1
等 级:新手上路
帖 子:12
专家分:0
注 册:2007-11-8
收藏
 问题点数:0 回复次数:0 
如何对非线形方程的二分法求根进行改进?

在C语言中用二分法求解已知,但如何对非线形方程的二分法求根进行改进?
#include<stdio.h>
#include<math.h>
#include<malloc.h>
void main()
{
int i,j,k;
double x,e,a,b;
double *A,*B,*X,*f;
printf("qing shu ru qu jian de xia jie he shang jie:");
scanf("%lf%lf",&a,&b);
printf("qing shu ru shi yan wu cha:");
scanf("%lf",&e);
k=(int)(log((b-a)/e)/log(2)-1)+1;
printf("k=%d\n",k);
A=(double *)malloc(k*sizeof(double ));
B=(double *)malloc(k*sizeof(double ));
f=(double *)malloc((k-1)*sizeof(double ));
X=(double *)malloc(k*sizeof(double ));
A[0]=a;
B[0]=b;
for(i=0;i<k;i++)
{
X[0]=(A[0]+B[0])/2;
if((pow(X[0],3+0.0)-X[0]-1)==0.0) x=X[0];
else if((pow(A[i],3+0.0)-A[i]-1.0)*(pow(X[i],3+0.0)-X[i]-1.0)<0)
{
A[i+1]=A[i];
B[i+1]=X[i];
X[i+1]=(A[i+1]+B[i+1])/2;
}
else if((pow(B[i],3+0.0)-B[i]-1)*(pow(X[i],3+0.0)-X[i]-1)<0)
{
A[i+1]=X[i];
B[i+1]=B[i];
X[i+1]=(A[i+1]+B[i+1])/2;
}
}
/* for(i=0;i<=k;i++)
{
printf("%lf\n",A[i]);
printf("%lf\n",B[i]);
printf("%lf\n",X[i]);
} */
printf("yong er fen fa qiu de de jie guo wei:\n");
printf("%lf\n",x=X[6]);
}

搜索更多相关主题的帖子: 二分法 方程 非线形 改进 
2007-11-08 19:57
快速回复:如何对非线形方程的二分法求根进行改进?
数据加载中...
 
   



关于我们 | 广告合作 | 编程中国 | 清除Cookies | TOP | 手机版

编程中国 版权所有,并保留所有权利。
Powered by Discuz, Processed in 0.026620 second(s), 9 queries.
Copyright©2004-2024, BCCN.NET, All Rights Reserved