c语言公式法求解一元三次方程

在这里插入图片描述

#include<stdio.h>
#include<math.h>
int main() {
	float a,b,c,d,x1,x2,x3,p,q,i,w,s1,s2,s3,X2,X3;
	printf("请输入一元三次方程标准形式ax^3+bx^2+cx+d=0\n");
	printf("a b c d的值,用空格隔开\n");
	scanf("%f%f%f%f",&a,&b,&c,&d);
	p=(3*a*c-b*b)/(3*a*a);
	q=(27*a*a*d-9*a*b*c+2*b*b*b)/(27*a*a*a);
	w=-0.5;
	s1=-b/(3*a);
	if(-q/2+sqrt(q*q/4+p*p*p/27)>0) {
		s2=pow((-q/2+sqrt(q*q/4+p*p*p/27)),1.0/3);
	} else {
		s2=-pow(-(-q/2+sqrt(q*q/4+p*p*p/27)),1.0/3);
	}
	if(-q/2-sqrt(q*q/4+p*p*p/27)>0) {
		s3=pow((-q/2-sqrt(q*q/4+p*p*p/27)),1.0/3);
	} else {
		s3=-pow(-(-q/2-sqrt(q*q/4+p*p*p/27)),1.0/3);
	}
	x1=s1+s2+s3;
	x2=s1+w*s2+w*s3;
	X2=sqrt(3)/2*s2+(-sqrt(3)/2)*s3;
	x3=s1+w*s2+w*s3;
	X3=sqrt(3)/2*s3+(-sqrt(3)/2)*s2;
	printf("第一个解为:%f\n",x1);
	if(X2>0) {
		printf("第二个解为:%f+%f*i\n",x2,X2);
	} else {
		printf("第二个解为:%f%f*i\n",x2,X2);
	}
	if(X3>0) {
		printf("第三个解为:%f+%f*i\n",x3,X3);
	} else {
		printf("第三个解为:%f%f*i\n",x3,X3);
	}
}
  • 6
    点赞
  • 26
    收藏
    觉得还不错? 一键收藏
  • 4
    评论

“相关推荐”对你有帮助么?

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论 4
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值