數(shù)值逼近課程設(shè)計(jì)-閱讀頁(yè)
2025-06-22 22:50本頁(yè)面
【正文】 { if(xxx[n]) { s=remezf(x[n])。 for(k=n2。k) t=t*x[n]+p[k]。 if(s*yy)x[n]=xx。k=n1。 x[n]=xx。 j=j1。 t=p[n1]。k=0。 s=st。 else x[j]=xx。 a=。 eps=。 printf(a0=%e\n,p[0])。 printf(a2=%e\n,p[2])。 printf(MAX(pf)=%e\n,p[3])。 } 三、運(yùn)行結(jié)果及說(shuō)明： 其中 a0,a1,a2 分別為二次多項(xiàng)式零次，一次，二次的系數(shù)。 程序代碼 (C++實(shí)現(xiàn) )： include include include include double f(double x) //被積分函數(shù) f(x) { return exp(x)。 else return 2*fun(n1)。 } double T(int n) //復(fù)化梯形公式 { double a=0,b=1,sum=0,h=(ba)/n。 for(i=1。i++) { sum+=f(a+i*h)。 return (h*sum)。 Sn=(4*T(2*n)T(n))/3。 } double C(int n) //復(fù)化 Cotes公式 { double Cn。 return (Cn)。 cout1復(fù)化梯形公式 ,2復(fù)化 Simpson公式 ,3復(fù)化 Cotes公式 :\n 。i=10。 coutsetiosflags(ios::showpoint)setprecision(16)。 } }三、運(yùn)行結(jié)果及說(shuō)明：
畢業(yè)設(shè)計(jì)相關(guān)推薦
鄂ICP備17016276號(hào)-1