程序師世界是廣大編程愛好者互助、分享、學習的平台,程序師世界有你更精彩!
首頁
編程語言
C語言|JAVA編程
Python編程
網頁編程
ASP編程|PHP編程
JSP編程
數據庫知識
MYSQL數據庫|SqlServer數據庫
Oracle數據庫|DB2數據庫
 程式師世界 >> 編程語言 >> JAVA編程 >> 關於JAVA >> 數值計算,利用梯形法和拋物線法計算圓周率

數值計算,利用梯形法和拋物線法計算圓周率

編輯:關於JAVA

公式pai=4*∫10dx/(1+x2)

把[0.1]分成10等份,用兩中方法分別計算

=======================================================

class work2 { public static void main (String args[]){ double h = 0.1;//參數h double m1 =0 ;//復化梯形求和求和 double n =10; double I=0;//最後結果 double H_2 =0.2; double m =5 ;//由於n=2m.n=10,所以m=5 double s1=0;//復化拋物線第一個數組求和 double s2=0;//復化拋物線第二個數組求和 double ss[] = new double[10];//復化梯形公式計算中的求和單元 double ss1[]=new double[5];//復化拋物線中用到的第一個數組 double ss2[]=new double[4];//復化拋物線中用到的第二個數組

double k[]={0,0.1,0.2,0.3,0.4,0.5,0.6,0.7,0.8,0.9};//10等分[0,1] for(int i=0;i<10;i++){ ss[i] = 1/(1+k[i]*k[i]); } //teat pass for(int i =0;i<10;i++){ m1 = m1 +ss[i]; } //test pass I = H_2*(0.5+2*m1);//0.5是f(b)

// --------------------test code System.out.println("===復化梯形公式計算(包含測試及結果)==="); System.out.println("復化梯形公式單項結果"); for(int i =0;i<10;i++){ System.out.println(ss[i]);

} System.out.println("結果求和 m is "+ m1); System.out.println("結果是 "+ I); System.out.println("==================================");// ------------------------------------------------ //==開始復化拋物線公式計算====== System.out.println("復化拋物線公式計算"); //由於拋物線運算對數組的操作比較復雜,所以我重新定義數組k[] double k_2[] = {0,0,0.1,0.2,0.3,0.4,0.5,0.6,0.7,0.8,0.9}; //這樣數組k就和實際中的對應起來了,哦哈哈 for(int i =1;i<=m;i++){ int n1 = 2*i-1; ss1[i-1] = 1/(1+k_2[n1]*k_2[n1]); System.out.println(ss1[i-1]);//測試計算是否正確, } System.out.println("運算正確"); System.out.println("第二個數組計算"); for(int i = 1;i<=(m-1);i++){ int n1 = 2*i; ss2[i-1]=1/(1+k_2[n1]*k_2[n1]); System.out.println(ss2[i-1]);//測試計算是否正確, } System.out.println("運算正確"); //對2個數組求和 for(int i =0;i<5;i++){ s1 = s1 +ss1[i]; } System.out.println("ss1求和是"+ s1); for(int i =0;i<4;i++){ s2 = s2 +ss2[i]; } System.out.println("ss2求和是"+ s2); System.out.println("運算正確"); I = (h/3)*(1+0.5+4*s1+2*s2); System.out.println("結果是 "+ I); //拋物線法不知道正不正確,書上沒有答案 } } ===========================================

  1. 上一頁:
  2. 下一頁:
Copyright © 程式師世界 All Rights Reserved