鎷夋牸鏈楁棩鎻掑煎叕寮忥細
聽聽聽聽聽聽聽聽聽聽聽聽聽聽聽聽聽聽聽聽聽聽 聽n聽聽聽聽聽聽聽聽聽n
聽 Pn(x(i))=聽聽聽聽聽 鈭戙斅犅?鈭徛?錛坸-x(j)錛?錛坸(k)-x(j)錛?銆晊(k)
聽聽聽聽聽聽聽聽聽聽聽聽聽聽聽聽聽聽聽聽聽聽k=0聽聽聽 聽j=0
聽聽聽聽聽聽聽聽聽聽聽聽聽聽聽聽聽聽聽聽聽聽聽聽聽聽聽聽聽聽聽聽 j鈮爇
聽
聽 灞炴э細鎻掑艱綆楁硶
聽聽聽聽聽聽聽聽聽聽聽聽聽聽聽聽聽聽聽聽聽聽聽聽聽聽聽聽聽聽聽聽聽聽聽聽聽聽聽聽聽聽聽聽聽聽聽聽聽聽聽聽聽聽聽聽聽聽聽聽聽聽聽聽聽聽聽聽聽聽聽聽聽聽聽聽聽聽聽聽聽聽聽聽聽聽聽聽聽聽 n
聽 綺懼害(灞閮ㄦ埅鏂宸?錛殀 f(x) - Pn(x) | = [f(蔚)] / (n+1)!聽 鈭徛?( x - x(k) )聽聽 (娉細鍏朵腑[f(蔚)]涓篺(蔚)絎琻+1嬈℃眰瀵肩殑琛ㄨ揪寮?
聽聽聽聽聽聽聽聽聽聽聽聽聽聽聽聽聽聽聽聽聽聽聽聽聽聽聽聽聽聽聽聽聽聽聽聽聽聽聽聽聽聽聽聽聽聽聽聽聽聽聽聽聽聽聽聽聽聽聽聽聽聽聽聽聽聽聽聽聽聽聽聽聽聽聽聽聽聽聽聽聽聽聽聽聽聽聽聽聽k=0
聽聽聽聽聽聽聽聽聽聽聽聽聽聽聽聽聽聽聽聽聽聽聽聽聽聽聽聽聽聽聽聽聽聽聽聽聽聽聽聽聽
聽銆婃暟鍊煎垎鏋愮畝鏄庢暀紼嬨?2 Editon -楂樼瓑鏁欒偛鍑虹増紺?-page 18 -綆楁硶嫻佺▼鍥?
聽 浠ｇ爜緇存姢錛?005.6.11聽 DragonLord
**/
#include<iostream.h>
int main()
{
聽float x;//鎻掑?br />聽float p[10][2];//宸茬煡(x0,y0),(x1,y1)...
聽int n;//杈撳叆宸茬煡鎻掑肩粍鏁?br />聽float y,t;
聽int k;
聽
聽cout<<"杈撳叆鎻掑肩粍鏁?"<<endl;
聽cin>>n;
聽cout<<"杈撳叆"<<n<<"緇勫凡鐭ユ彃鍊兼暟錛圶,Y錛?<<endl;
聽聽聽 for(int i=0;i<n;i++)
聽{
聽聽cin>>p[i][0]>>p[i][1];
聽}
聽聽聽 cout<<"杈撳叆鎻掑?"<<endl;
聽cin>>x;
聽
聽y=0;
聽k=0;
聽
聽for(k=0;k<n;k++)
聽{
聽聽t=1;
聽聽for(i=0;i<n;i++)
聽聽{
聽聽聽if(i!=k)
聽聽聽{
聽聽聽聽t=t*(x-p[i][0])/(p[k][0]-p[i][0]);聽 //key step
聽聽聽}
聽聽}
聽聽y=y+p[k][1]*t;
聽}
聽cout<<"鎻掑肩粨鏋滐細"<<y<<endl;
return 0;
}
]]>