Aの一番右が解になっている。sys.stat.hがないと機能しない。その場合は手動でanimeのフォルダを作ってコメントアウトする。"\t%d\t%f\n",i,uw[i]);}fclose(file);}まあ、だいたいコメントに書いておきました。今回は2次元配列を使っています。2次元配列に、行列の係数と右辺を格納しています。A[i][j]は(0から数えて)i行j列目です。その2次元配列をヤコビの反復法で解きます。式と色々見比べてみてください。Cでトーマス法で書いたのもgibhubに上げておいたので、生成物を見比べて見ると面白いと思います。max_iter=100とかにすると、かなり雑で誤差がだいぶ出てきます。コード置いときますまあ前回と全く同じプロセスですが、計算を回してみます。-g./a.outで実行できるかと思います。オプションはデバッグ用の情報を追加で記載する設定です。./a.outとかでデバッカーを使えます。普通gdbなのですが、Macがサポートしてないので、lldbです。実行すると、animeというフォルダが作成されて、そこにファイルが書き出されると思います。前回作った、anime.gpを使ってアニメーションも書いてみてください。前回と同じですが、一応アニメーションのプログラムを。今回はC++じゃなくて、Cverを貼っときます。以下のプログラムをheat_animetion_make.Cとでもして保存してください。メッシュ数に合わせる。メッシュ数に合わせる。1\n");}}fclose(file);}こいつで、anime.gpといった奴を作り、それを使ってgnuplotで簡単にアニメーションを書きます。anime.gp陰解法の方が精度が高いので、メッシュ数とかdtとか色々いじってみると面白いです。コメントをキャンセルメールアドレスが公開されることはありません。