定積分的montecarlo計(jì)算方法的實(shí)現(xiàn)-資料下載頁(yè)
2025-05-15 07:07本頁(yè)面
【正文】 e s,t。 double i。 long ko1,ko2。 ko1=ko2=1。 for(i=0。in。i++) { t=ran2(amp。ko1)。 s=ran2(amp。ko2)。 x=a+(ba)*t。 y=(1)*fMAX+2*fMAX*s。 if(y=0) { if(y==0) sy++。 if(y=(*p)(x)) sx++。 } else { if(y=(*p)(x)) vx++。 } } return ((sxvx)/n)。 } main() { double n,fMAX,a,b,Mont,A,B,h,i=0,N。 long time0,time1。 printf(Enter a=:)。scanf(%lf,amp。a)。 printf(\nEnter b=:)。scanf(%lf,amp。b)。 printf(\nEnter N=:)。scanf(%lf,amp。N)。 /*N表示區(qū)間的個(gè)數(shù) */ printf(\nEnter n=:)。scanf(%lf,amp。n)。 printf(\nEnter fMAX(此函數(shù)的界 )=:)。scanf(%lf,amp。fMAX)。 h=(ba)/N。 time0=clock()。 for(i=0。iN。i++) { A=a+i*h。 B=A+h。 Mont+=Monte(n,A,B,fun1,fMAX)*(BA)*2*fMAX。 } printf(\nintegration=%lf,Mont)。 time1=clock()。 printf(\n耗時(shí)： %ldms,time1time0)。 getch()。
環(huán)評(píng)公示相關(guān)推薦
鄂ICP備17016276號(hào)-1