Oracle中實現圓周率計算（三）

今天兩個同事用JAVA實現圓周率一百位小數的實現。一個同事問我要不要試試，由於很長時間沒有寫過JAVA程式碼，而且本身JAVA的水平就很差，於是打算用ORACLE實現。

這篇給出一個真正的演算法。

Oracle中實現圓周率計算（一）：http://yangtingkun.itpub.net/post/468/468799

Oracle中實現圓周率計算（二）：http://yangtingkun.itpub.net/post/468/468870

 

 

上一篇文章提到，Oracle的NUMBER型別精度只有38位，因此想要確保100位的精度就無法直接使用Oracle的NUMBER型別來實現了。

如果在Oracle希望實現這個功能，那麼最簡單的辦法莫過於使用JAVA儲存過程了，下面演算法不變，將程式碼遷移到JAVA儲存過程中：

SQL> create or replace and compile java source named "pi" as
  2  import java.math.BigDecimal;
  3  public class pi extends Object
  4  {
  5   public static String pi (int n)
  6   {
  7    BigDecimal result = new BigDecimal(1);
  8    for (int i=n;i>0;i--)
  9    {
 10     BigDecimal temp = new BigDecimal(i);
 11     result = new BigDecimal(2).add
 12      (result.multiply
 13       (temp.divide(new BigDecimal(1).add
 14        (temp.multiply
 15         (new BigDecimal(2))
 16        ),102,BigDecimal.ROUND_HALF_UP)
 17       )
 18      );
 19    }
 20    return result.toString().substring(0,102);
 21   }
 22  }
 23  /

Java 已建立。

建立儲存過程後，在建立一個PL/SQL函式介面：

SQL> create or replace function f_pi (p_n in number) return varchar2 as
  2  language java name 'pi.pi(int) return String';
  3  /

函式已建立。

下面就可以呼叫這個JAVA儲存過程進行運算了：

SQL> SELECT F_PI(10) FROM DUAL;

F_PI(10)
--------------------------------------------------------------------------------------------------------
3.1408420956408572507643715074055631331173126838761513684114303309349748978231950368173278389996656250

SQL> SELECT F_PI(100) FROM DUAL;

F_PI(100)
--------------------------------------------------------------------------------------------------------
3.1415926535897932384626433832792952864908441267910666216977616133327809721178539595244159254337299498

在多項式N達到330的時候，返回的圓周率精度到了100位。

SQL> SELECT F_PI(329) FROM DUAL;

F_PI(329)
--------------------------------------------------------------------------------------------------------
3.1415926535897932384626433832795028841971693993751058209749445923078164062862089986280348253421170678

SQL> SELECT F_PI(330) FROM DUAL;

F_PI(330)
--------------------------------------------------------------------------------------------------------
3.1415926535897932384626433832795028841971693993751058209749445923078164062862089986280348253421170679