數(shù)值計(jì)算(二分法、簡單迭代法、newton迭代法、弦截法(割線法、雙點(diǎn)弦法))(編輯修改稿)
2025-07-25 14:26 本頁面
【文章內(nèi)容簡介】 //根的個(gè)數(shù)int x_count。 double precision=。 //精度要求//函數(shù)的表達(dá)式（x^53x^3+x1）double f(double x){ return (pow(x,5)3*pow(x,3)+x1)。}void init(){ //根據(jù)函數(shù)圖像確定根的區(qū)間和迭代初值 r[0]=。 r[1]=1。 r[2]=。 rn[0]=。 rn[1]=。 rn[2]=。}//尋找根的區(qū)間 void search(){ //若沒有給出區(qū)間和初值，進(jìn)行逐步搜索有根區(qū)間 for(int i=0。i*foot88。i++){ if(f(i*foot8)*f((i+1)*foot8)0){ rn[m]=i*foot8。 m++。 } } }//=====================二分法==========================double Dichotomy (double a,double b){double mid=0。int i=0。while(fabs(ba)precision){ mid = (a+b)/2。 if(f(a)*f(mid)=0) b=mid。 //判斷與端點(diǎn)函數(shù)值得符號(hào) else a=mid。 coutmidendl。 } r[x_count++]=mid。 return mid。 //返回最終結(jié)果 }//================簡單迭代法========================= //構(gòu)造迭代公式double fitera(double x){ double result=0。 double xx=3*pow(x,3)x+1。 if(xx=0){ xx=xx。 return pow(xx,)*(1)。} else return pow(xx,)。 }//簡單迭代double itera(double x0){ coutx0endl。 double x1=fitera(x0)。 while(fabs(x1x0)precision){ x0=x1。 x1=fitera(x0)。 //沒有到達(dá)精度要求繼續(xù)迭代 coutx1endl。 } re
醫(yī)療健康相關(guān)推薦
鄂ICP備17016276號(hào)-1