分段线性插值

/**
***分段线性插值***

  计算公式：S1(x)=Φ0（x-x(i)/h(i))）y(i)+Φ1（x-x(i)/h(i)）y(i+1)
 
  步长：h(i)=x(i+1)-x(i)
 
  领域：x(i)≤x≤x(i+1)

  属性：插值计算法

  精度：| f(x) - S1(x) | = h(i)*h(i)/8 *      max     |  f''(x)  |
                                                          x(i)≤x≤x(i+1) 

  《数值分析简明教程》-2 Editon -高等教育出版社 -page 31 -无算法流程图

  代码维护：2005.6.11  DragonLord
**/

#include<iostream.h>
#include<math.h>
#include<stdio.h>
int main()
{
 double x0,x1,y0,y1,u0,u1,x;
 double result;
 while(cin>>x0>>y0>>x1>>y1)
 {
        cin>>x;
        u1=(x-x0)/(x1-x0);
        u0=(x-x1)/(x0-x1);
        result=u1*y1+u0*y0;
  printf("%.5f",result);
 }
 return 0;
}