抛物线插值.txt

该压缩文件有一些经典插值c语言版的源程序

#include<stdio.h>
#include<math.h>



main()
{
float a,b,t,s,q,fa,fb,ft,fs,m,n;
clrscr();
printf("input the a ,b,t,q:\n");
scanf("%f,%f,%f,%f",&a,&b,&t,&q);
while(1)
{
fa=8*a*a*a-2*a*a-7*a+3;
fb=8*b*b*b-2*b*b-7*b+3;
ft=8*t*t*t-2*t*t-7*t+3;
s=((t*t-b*b)*fa+(b*b-a*a)*ft+(a*a-t*t)*fb)/(2*((t-b)*fa+(b-a)*ft+(a-t)*fb));
fs=8*s*s*s-2*s*s-7*s+3;
m=s-t;
n=fabs(m);
printf("s=%f\t t=%f\t m=%f\t f(t)=%f\t f(s)=%f\n",s,t,m,ft,fs);
if(n>q)
{
if(m>0)
{
if(fs<=ft)
{
a=t;
t=s;
}
else b=s;
}
else
{
if(fs<=ft)
{
b=t;
t=s;
}
else a=s;
}
}

else break;

}
printf("the result is:\n");
printf("s=%f\n",s);
printf("f(s)=%f",fs);
}