| 查看: 2058 | 回复: 5 | ||
[求助]
求高手帮我编一个一元三次方程求解的小程序！重谢哦！
|
| 求哪位程序高手，能够帮我编写一个C语言程序解一元三次方程！！重谢！ |
回复此楼» 猜你喜欢
332求调剂
已经有8人回复
-
材料工程085601，270求调剂
已经有29人回复
-
材料与化工调剂
已经有14人回复
-
291 求调剂
已经有17人回复
-
机械专硕270求调剂，接受跨专业
已经有11人回复
-
求调剂
已经有8人回复
-
本科211 工科085400 280分求调剂 可跨专业
已经有8人回复
-
化学工程与技术324调剂
已经有23人回复
-
材料085601调剂
已经有15人回复
-
269求调剂
已经有10人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 求高手帮忙鉴定一下植物
已经有14人回复
- 请帮忙翻译一句话，谢谢！！！
已经有2人回复
- 求解 2N HCl,是什么意思？
已经有22人回复
- 物化考博七题求解，谢谢
已经有10人回复
- 核磁/质谱解谱高手帮忙看看这个化合物结构怎么解析
已经有12人回复
- 镁基复合材料的XRD图 高手帮忙看一下啊 急！！！
已经有12人回复
- 高手帮忙翻译下！谢谢！
已经有1人回复
- 关于xp安装ipv6的问题，盼高手帮忙~~~
已经有6人回复
- 求高手精修晶体并得到氢键表，本人非此专业。烦请高手求解
已经有9人回复
- 新虫求解文献里的一小段英文，求高手翻译．．．
已经有12人回复
- 先谢过了(求解一个物理的问题，急需数学大虾的帮助）
已经有7人回复
- 写了一个fortran90的小程序，编译通不过，请大侠帮忙
已经有59人回复
- 【求助】50金币求解一个晶体，高手进
已经有22人回复
- 【求助】求助：vb编程中用牛顿迭代解三次方程为什么只得到一个根？【已完成】
已经有15人回复
- 高手求解将 -S- 氧化成亚砜怎么做？
已经有19人回复
锐利的碎片
木虫 (正式写手)
star watcher
- 应助: 136 (高中生)
- 金币: 3637.1
- 散金: 252
- 红花: 22
- 帖子: 988
- 在线: 1224.9小时
- 虫号: 961933
- 注册: 2010-03-05
- 专业: 凝聚态物性 II ：电子结构
2楼2011-10-27 12:52:11
3楼2011-10-27 15:26:40
noryark
木虫 (正式写手)
- 应助: 0 (幼儿园)
- 贵宾: 0.5
- 金币: 3658.4
- 散金: 1338
- 红花: 20
- 帖子: 873
- 在线: 479.2小时
- 虫号: 558815
- 注册: 2008-05-15
- 专业: 放射医学
4楼2011-10-27 18:46:46
capital_j
木虫 (小有名气)
- 应助: 7 (幼儿园)
- 金币: 3599.4
- 散金: 141
- 帖子: 116
- 在线: 68.7小时
- 虫号: 1720641
- 注册: 2012-03-27
- 性别: GG
- 专业: 可再生与替代能源利用中的
【答案】应助回帖
★
感谢参与，应助指数 +1
jjdg: 金币+1, 感谢参与 2012-04-02 03:38:34
感谢参与，应助指数 +1
jjdg: 金币+1, 感谢参与 2012-04-02 03:38:34
|
#include #include using namespace std; int isFloatsEqual(float a,float b,float c=1e-2) { if(fabs(a-b) } int func(float a,float b,float c,float d) { if(a==0) return 1; float A = b*b-3*a*c; float B = b*c-9*a*d; float C = c*c-3*b*d; float derta = B*B - 4*A*C; cout< cout<<"A =B 公式一：有三个重根\n X1,2,3 = "<<-b/3/a< } if(derta>0){ float Y1 =A*b+1.5*a*(-B+pow(derta,0.5)); float Y2 =A*b+1.5*a*(-B-pow(derta,0.5)); cout<<"Δ>0 公式二：有一个实根和两个共轭虚根\n"; cout< float vir = pow(3,0.5)*(pow(Y1,1.0/3)-pow(Y2,1.0/3))/6/a; cout<<"X2 = "< } if(isFloatsEqual(derta,0)) { if(isFloatsEqual(A,0)) return 1; float K = B/A; cout<<"Δ=0 公式三：有三实根包括两个重根\n"; cout<<"X1 = "<<-b/a+K< } if(derta<0){ float T = (2*A*b-3*a*B)/(2*pow(A,1.5)); float degree = acos(T); cout<<"Δ<0 公式四：有三不相等实根\n"; float squ1 = (-b+pow(A,0.5)*cos(degree/3))/3/a; float squ2 = pow(3*A,0.5)*sin(degree/3)/3/a; cout<<"X1 = "<<(-b-2*pow(A,0.5)*cos(degree/3))/(3*a)< } } int main() { int a1=1,a2=67.4,a3=1417.92,a4=-9539.712; while(a1) { func(a1,a2,a3,a4); cout<<"= = = = = = = = = = = = = = ="< } return 0; } 利用盛金公式。参见百度百科。但是没有严格的得到学术界的认可。 程序没有经过优化，只是将盛金公式以C++的形式展现出来。 主要就看个思路吧 |
5楼2012-04-02 00:10:04
emanlee
木虫 (小有名气)
- 应助: 28 (小学生)
- 金币: 3521.6
- 散金: 100
- 红花: 2
- 帖子: 116
- 在线: 125.6小时
- 虫号: 1466309
- 注册: 2011-10-29
- 性别: GG
- 专业: 计算机应用技术
【答案】应助回帖
★
感谢参与，应助指数 +1
dubo: 金币+1, 欢迎讨论 2012-04-06 19:01:20
感谢参与，应助指数 +1
dubo: 金币+1, 欢迎讨论 2012-04-06 19:01:20
|
#include #include int a,b,c,d; int count=0; double x[3]={-200,-200,-200}; void fun(double x1,double x2) { double x3,y1,y2,y3; int i=0; while (x2-x1>10e-6) { x3=(x1+x2)/2; y1=a*x1*x1*x1+b*x1*x1+c*x1+d; y2=a*x2*x2*x2+b*x2*x2+c*x2+d; y1=a*x1*x1*x1+b*x1*x1+c*x1+d; y3=a*x3*x3*x3+b*x3*x3+c*x3+d; if (y1==0) { for (i=0;i<3;i++) { if (x==x1) { return; } if (x==-200) { count++; x=x1; return; } } } else if (y2==0) { for (i=0;i<3;i++) { if (x==x2) { return; } if (x==-200) { count++; x=x2; return; } } } else if (y3==0) { for (i=0;i<3;i++) { if (x==x3) { return; } if (x==-200) { count++; x=x3; return; } } } else if (y1*y3<0) { x2=x3; } else { x1=x3; } } for (i=0;i<3;i++) { if (x==x1) { return; } if (x==-200) { count++; x=x1; return; } } } int main() { double x1,x2,y1,y2,x3,x4; scanf("%d,%d,%d,%d",&a,&b,&c,&d); for (x1=-100,x2=x1+1;x1<100;x1+=1,x2+=1) { y1=a*x1*x1*x1+b*x1*x1+c*x1+d; y2=a*x2*x2*x2+b*x2*x2+c*x2+d; x3=x1; x4=x2; if (y1*y2<=0) { fun(x3,x4); if (count==3) { break; } } } for(int i=0;i<3;i++) { printf("%.2lf ",x); } printf("\n"  ;return 0; } |
6楼2012-04-03 08:42:01
5