多倍長クラスのCLNについてなのですが (ID:67225)

名前

ホームページ(ブログ、Twitterなど)のURL (省略可)

本文 多数の投稿ありがとうございました。 ヤフーで質問もしたのですが、解答を得られず コメントを頂いたのですが、ある程度解決し、そのソースを載せようと思ったのですが、 二つも同じスレをたてるのも、どうかと思い削除してしまい 色々批判を受けて両方削除しました。 決して逃亡したわけでも、ムカついたから削除した等ではないのでご理解ください。 そこでCLNのソースなんですが、 #include <cln/cln.h> using namespace std; using namespace cln; int main() { float_format_t precision = float_format(100); int n = 0; //n次の値を出すためのn double xf = 10; //xをどこまで計算させるか double xia = 0.1; //x start2 // ● cl_F x = "0.1_60"; //x start // ● double delta = 1.0E-2; //x plus // ■ cl_F xd = "1.0E-2_60"; //x plus // ■ double xi; cl_F j[52],scale,sum; //k=50まで計算させる配列,scale cl_F Fk = "50.0_60"; //kと同期させるFk cl_F Fk2 = "50.0_60"; //Fkの初期化用 cl_F b = "1.0_60"; // 数値1 cl_F c = "2.0_60"; // 数値2 j[51] = j[50] = "1.0_60"; cout << "x" << '\t' << "y" << endl; for(xi=xia;xi<=xf;xi=xi+delta){ //xは0.1→40．0まで Fk = Fk2; // Fk=50で初期化している for(int k=50;k>0;k--){ // kは50→0.0まで計算 j[k-1] = (((c * Fk + b)/x) * j[k]) - j[k+1]; // 漸化式 Fk = Fk - b; // Fk と k を同期させている } scale = ((sin(x))/x)/j[0]; // j0=sinx/x sum = j[n]*scale; cout << x << '\t' << sum << endl; x = x + xd; } } このプログラムをすっきりさせる方法はありますでしょうか？ Fkなど二つあるのは、CLNの数値は1.0など加算してしまうと 桁落ちが発生してしまうのを防ぐためにあります。 そこで質問なのですが、y→0 のｘを求めたいのですが（ｘの零点） どのようにしたらいいのでしょうか？ if(sum<xd && sum >0.0){cout << x << endl;} の文を考えてみたのですが、いまいちきれいな形とはいえないと思いますので・・・ 他に何か助言がありましたらお願いします。