目前某叶编的最难的程序了......感觉算是跨越吧,之前最难的是一元二次方程求解,虽然是最“难”的,只是因为最长,但是写起来还是很轻松的
不过一元三次方程可一点都不轻松,很累,因为没学过一元三次方程解法,所以需要百度,不过在C语言学习的书上找到了解法,要用牛顿迭代法求......
百度娘处充电,勉勉强强算是能看懂牛顿迭代法了,可能根本没看懂,不过至少知道怎么算了
解法思路:
先把求与X轴交点坐标公式放着免得忘记了
x= x1f(x2)-x2f(x1)/f(x2)-f(x1)
之后比较x1的y1值和x2的y2值,如果两个为异号,那么两个x之间一定有方程的根
如果同号,那么继续输入直到异号为止
这个时候用求交点坐标公式求出交点坐标x,它的y值同样代入求出
再次比较y与y1值,如果异号那么x与x1之间必有方程根
如果同号那么x与x2之间必有方程根
循环以上直到y的绝对值小于一个非常小的数,也就近似为0的时候,输出x值既为方程根......
理好了解法思路就不怕了,因为这个程序的关键就是我没学过解法,现在解法清楚了就开始编了
写了很多的函数,对我来说很复杂,所以注释比较多,省的以后自己都看不懂了......
书上给的例子方程式x^3-5x^2+16x-80=0给出的一个解x=5,图2反应出来了......
#include <stdio.h>
#include <math.h>
#include <conio.h>
float a,b,c,d; //定义外部变量,使全局可以调用
float f(float x) //x函数
{
floaty;
y=a*x*x*x+b*x*x+c*x+d;
returny;
}
float xpoint(float x1,float x2) //求弦与x轴交点坐标
{
floaty;
y=(x1*f(x2)-x2*f(x1))/(f(x2)-f(x1));
returny;
}
float root(float x1,float x2) //求根函数
{
floatx,y,y1;
y1=f(x1);//y1为x1纵坐标
do
{
x=xpoint(x1,x2); //求x1与x2之间弦与x轴交点赋值于x
y=f(x); //代入方程中求得y
if(y*y1>0) //判断y与y1是否同号
{
x1=x;
y1=y;