數(shù)值逼近課程設(shè)計(jì)-其他專(zhuān)業(yè)-資料下載頁(yè)
2025-01-19 11:49本頁(yè)面
【導(dǎo)讀】對(duì)于Runge函數(shù)f=，在[-1,1]上作等距節(jié)點(diǎn)插值，分別取n=4,n=8,n=12,編出程序，畫(huà)出此插值的圖像。求函數(shù)f=在[-1,1]上的二次多項(xiàng)式逼近。其中a0,a1,a2分別為二次多項(xiàng)式零次，一次，二次的系數(shù)。
【正文】 次多項(xiàng)式零次，一次，二次的系數(shù)。 作業(yè) 三 復(fù)化公式 用復(fù)化公式計(jì)算積分 10 dxex?。 程序代碼 (C++實(shí)現(xiàn) )： include include include include double f(double x) //被積分函數(shù) f(x) { return exp(x)。 } int fun(int n) { if(n=0) { if(n==0)return 1。 else return 2*fun(n1)。 } else exit(0)。 } double T(int n) //復(fù)化梯形公式 { double a=0,b=1,sum=0,h=(ba)/n。 long i。 for(i=1。in。i++) { sum+=f(a+i*h)。 } sum+=(f(a)+f(b))/2。 return (h*sum)。 } double S(int n) //復(fù)化 Simpson公式 { double Sn。 Sn=(4*T(2*n)T(n))/3。 return (Sn)。 } double C(int n) //復(fù)化 Cotes公式 { double Cn。 Cn=(16*S(2*n)S(n))/15。 return (Cn)。 } void main() { int i。 cout1復(fù)化梯形公式 ,2復(fù)化 Simpson公式 ,3復(fù)化 Cotes公式 :\n 。 for(i=1。i=10。i++) { coutfun(i) 。 coutsetiosflags(ios::showpoint)setprecision(16)。 coutT(fun(i)) S(fun(i)) C(fun(i))endl。 } }三、運(yùn)行結(jié)果及說(shuō)明：
黨政相關(guān)相關(guān)推薦
鄂ICP備17016276號(hào)-1