C/C++の宿題を片付けます 100

[1] 授業単元：プログラミング
[2] 問題文(含コード&リンク)：
多倍長精度の計算をするためにgmpライブラリを使用して、
ボールウェイン(Bolwein)の４次式で円周率を求めるプログラムを作成せよ。
アルゴリズム
1. 初期化 a=6 - 4√2 ,y=√2 - 1 ,x=2
2. yが0にならない間以下を繰り返す
(1) u = √ { √(1-y^4) }
(2) y =(1-u)/(1+u)
(3) v =(1+y)
(4) w =v*v*v*v
(5) x = x*4
(6) a=a*w - y*(1+y+y*y)*x
3. 円周率π = 1/a

この長精度整数用アルゴリズムを構成するには、上記公式の中の「定数1」に対応する整数Mを設定しなければならない。
このMを基準として、各値を計算する。
長精度型の変数をA,Y,U,V,W,　通常の整数型の変数をxとおくとき、Mで上の公式を書き直し、テキストに記せ。
ヒント：aの初期化では A=6・M-4・√（M・M・2）　となる。　

[3] 環境
　[3.1] OS： XP　pro
　[3.2] gcc
　[3.3] 言語： C
[4] 期限： 12月6日
[5] その他の制限： 特になし
よろしくお願いします。