甜橘子，专注成都网站制作网站设计与营销型网站建设与优化
QQ :532337155
电话：028-86922220
我只看了你的第二个
主要从事网页设计、PC网站建设（电脑版网站建设）、wap网站建设（手机版网站建设）、自适应网站建设、程序开发、微网站、小程序开发等，凭借多年来在互联网的打拼，我们在互联网网站建设行业积累了丰富的网站设计、成都做网站、网络营销经验，集策划、开发、设计、营销、管理等多方位专业化运作于一体，具备承接不同规模与类型的建设项目的能力。
你定义的f函数只有一个形参x，所以在main函数中输入的你所认为的实参abcd 的值无法转给f函数 所以运行时会出错 只要赋给abcd确定的值或增加f函数的形参
scanf("%f,%f,%f,%f",d,d,d,d);这句也错了应该是scanf("%f,%f,%f,%f",a,b,c,d);
S*cos(α+0.01410087457*S*cosα)=K
令b=0.01410087457，u=α+b*S*cosα，
则S=(u-α)/(bcosα)，原方程变为 (u-α)cosu=bKcosα
令f(u)=(u-α)cosu，则f'(u)=cosu-(u-α)sinu
∵f(u)=(u-α)cosu=bKcosα=常数，∴f'(u)=cosu-(u-α)sinu=0
联立上述两个方程，得 sin²u+bKcosα*sinu-1=0
∵△=(bKcosα)^2+4≥4，∴√△≥2
∴sinu=(-bKcosα+√△)/2 (另一根≤-1舍弃)
∴u=arcsin[(-bKcosα+√△)/2]
（若你的解不是唯一，则应取u=nπ±arcsin[(-bKcosα+√△)/2]，n∈Z）
∴S=(u-α)/(bcosα)={arcsin[(-bKcosα+√△)/2]-α}/(bcosα)
（其中△=(bKcosα)^2+4，b=0.01410087457；其余α，K均为已知）
用这个公式就可以直接算出你要求的S的值了
#includestdio.h
#includemath.h
void main()
{
float x0,x1,x2,fx0,fx1,fx2;
do
{
printf("enter x1 x2:");
scanf("%f,%f",x1,x2); //输入在哪个区间寻找方程的根
fx1=x1*((2*x1-4)*x1+3)-6; //计算在两个端点，函数的值
fx2=x2*((2*x2-4)*x2+3)-6;
}while(fx1*fx20); //必须要函数在两个端点，函数的值异号，才能保证区间内有根
do
{
x0=(x1+x2)/2; //x0为区间的中点
fx0=x0*((2*x0-4)*x0+3)-6; //中点的函数值
if((fx0*fx1)0) //将中点修改为函数值与中点函数值同号的端点
{
x2=x0;
fx2=fx0;
}
else
{
x1=x0;
fx2=fx0;
}
}while(fabs(fx0)=1e-5); //直到中点的函数值与0足够接近
printf("x=%6.2f\n",x0);
}
可以使用二分法，弦截法，牛顿法，以前写的，希望对LZ有帮助。
①弦截法：
：
用弦截法求方程F(x)=X^3-5X^2+16X-80=0的根。
方法如下：
（1）取两个不同的点x1,x2,如果f(x1)和f(x2)符号相反，则（x1,x2）区间内必有一个根。如果f(x1)和f(x2)同符号，则应改变x1,x2，直到f(x1)和f(x2)异号为止。注意x1和x2的值不应差太大，以保证（x1,x2）区间内只有一根。
（2）连接（x1,f(x1)）和（x2,f(x2)）两点，此线交x轴于x，x的坐标可以用下式求出：x=(x1*f(x2)-x2*f(x1))/f(x2)-f(x1),
再从x求出f(x);
（3）若f(x)和f(x1)同符号，则根必在(x,x2)内，此时将x作为新的x1;如果f(x)和f(x2)同符号，则表示根在(x1,x)，此时将x作为新的x2。
（4）重复步骤（2）和（3），知道|f(x)|e,e是一个很小的数，例如10^-6;
程序如下：
#includestdio.h
#includemath.h
float f(float x)
{
float y;
y=((x-5.0)*x+16.0)*x-80.0;
return y;
}
float xpoint(float x1,float x2)
{//求(x1,f(x1)和(x2,f(x2))连线与x轴的交点
float y;
y=(x1*f(x2)-x2*f(x1))/(f(x2)-f(x1));
return y;
}//xpoint函数
float root(float x1,float x2)
{
float x,y,y1;
y1=f(x1);
do
{
x=xpoint(x1,x2);
y=f(x);
if(y*y10)
{
y1=y;
x1=x;
}
else x2=x;
}while (fabs(y)=0.00001);
return x;
}//root函数
void main()
{
float x1,x2,f1,f2,x;
do
{
printf("input x1,x2:\n");
scanf("%f%f",x1,x2);
f1=f(x1);
f2=f(x2);
} while(f1*f2=0);
x=root(x1,x2);
printf("a root of equation is %8.4f\n",x);
getch();
}
②：牛顿迭代法：
牛顿迭代法
牛顿迭代法又称牛顿切线法，它采用以下方法求根：先任意设定一个与真实的根接近的值x0作为第一个近似根，由x0求出f(x0)，过(x0,f(x0))点做f(x)的切线，交x轴于x1，把它作为第二次近似根，再由x1求出f(x1)，再过(x1,f(x1))点做f(x)的切线，交x轴于x2，再求出f(x2)，再作切线……如此继续下去，直到足够接近真正的x为止。
其中f'(X0)是函数在X0处的斜率，也就是在X0处的导数。
代码如下：
#includestdio.h
#includemath.h
float f(float a,float b,float c,float d,float x)
{
float f;
f=((a*x+b)*x+c)*x+d;
return f;
}
float f1(float a,float b,float c,float x)
{
float f;
f=(x*3*a+2*b)*x+c;
return f;
}
float root(float a,float b,float c,float d)
{
float x0,x1=1;
do
{
x0=x1;
x1=x0-f(a,b,c,d,x0)/f1(a,b,c,x0);
}while(fabs(x1-x0)=1e-6);
return x0;
}
void main()
{
float a,b,c,d,x;
printf("input four float numbers:\n");
scanf("%f%f%f%f",a,b,c,d);
x=root(a,b,c,d);
printf("%.1fX^3+%.1fX^2+%.1fX+%.1f=0 its root near x=1.5 is :%.4f\n",a,b,c,d,x);
getch();
}
求sin的：参考下 #includestdio.h void main() { double x,a,b,sum=0; printf("请输入x的弧度值：\n"); scanf("%lf",x); int i,j,count=0; for(i=1;;i+=2) { count++; a=b=1; for(j=1;j=i;j++) { a*=x; b*=(double)j; } if(a/b0.0000001) break; else { if(count%2==0) sum-=a/b; else sum+=a/b; } } printf("%lf\n",sum); }
要利用根与根之差的绝对值=1这个条件。
下面是刚学的时候写的，你可以看看
#includestdio.h
#includemath.h
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-x110e-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;i3;i++)
{
if (x[i]==x1)
{
return;
}
if (x[i]==-200)
{
count++;
x[i]=x1;
return;
}
}
}
else if (y2==0)
{
for (i=0;i3;i++)
{
if (x[i]==x2)
{
return;
}
if (x[i]==-200)
{
count++;
x[i]=x2;
return;
}
}
}
else if (y3==0)
{
for (i=0;i3;i++)
{
if (x[i]==x3)
{
return;
}
if (x[i]==-200)
{
count++;
x[i]=x3;
return;
}
}
}
else if (y1*y30)
{
x2=x3;
}
else
{
x1=x3;
}
}
for (i=0;i3;i++)
{
if (x[i]==x1)
{
return;
}
if (x[i]==-200)
{
count++;
x[i]=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;x1100;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;i3;i++)
{
printf("%.2lf ",x[i]);
}
printf("\n");
return 0;
}