二次插值法C语言相关程序

#include<stdio.h>
#include<math.h>
float f(float t)
{ float z;
  z=t*t-10*t+36;
  return z;
}
void main()
{float e=0.01,c1,c2,ap,yp,a,y;
 float h0=0.03,h=h0,a1=0,a2=h,y1,y2,a3,y3;
 y1=f(a1);
 y2=f(a2);
 if(y2>y1)
{
 h=-h;a3=a1;y3=y1;
 a1=a2;y1=y2;a2=a3;y2=y3;
}
 a3=a2+h;y3=f(a3);
while(y3<y2)
{
 h=2*h;
 a1=a2;y1=y2;a2=a3;y2=y3;
 a3=a2+h;y3=f(a3);
}
printf("f(%f)=%f\nf(%f)=%f\nf(%f)=%f",a1,y1,a2,y2,a3,y3);

 y1=f(a1);
 y2=f(a2);
 y3=f(a3);
 c1=(y3-y1)/(a3-a1);