{{{このプログラムは陽解法、陰解法、クランク・ニコルソン法で解を求めるプログラムです。クランク・ニコルソン法の関数を追加したらオーバーフローしました。youkaihou(fai,t,dt,x,dx,xmax)!t秒後の温度を求める。名付けて、コンパイル。2.5884155414774054E-0024.1881457553708118E-0026.7706545758653772E-0020.108356796100975230.168428254787954400.249453271232849690.348310092710075290.457943076854455840.569644427034318170.675269324610401340.768402528453419720.844511207117318860.900597510572817890.934789482293640230.946056061403698800.934054932793998320.899061626359431880.841932242118456630.764072095966407130.360083586218275460.138334122223380870.273262007598271850.401461306479500790.519775399736733480.625291232254801540.715411352069695530.787918722078526200.841033530157858640.873462420817231130.884442579973002130.873785951751065330.841931580980350810.790013415168563380.719940398872163880.634457665769337150.537112080177459150.432005905338776640.323245143005442550.214133764527224460.10641095232193158一応、動いてるようですが。知らないんで、的外れかもしれませんが。http://ameblo.jp/mratassyu/手提鞄屋魚有店(てさげかばんやうおありてん)レスがついていないものを優先して、レスしています。時々、見当外れなレスをします。サブルーチンのcranknicholson(fai,t,dt,x,dx,xmax)関数内で、例えばa(1,1)=rを書き足してコンパイルすると、オーバーフローします。2次元配列を書き足すとオーバーフローします。わかる方お願いします。16億バイトの領域を必要とします。4ギガバイトで、約43億。ユーザプログラムに与えられるメモリはその半分以下かもしれません。環境やエラーメッセージなどの正確な情報をお願いします。16億バイトの領域を必要とします。引数なので呼び出し元の領域が参照されるだけで大きな領域は必要ないでしょう。1年前その通りでした。ありがとうございます。イメージが、ごっちゃになっていた。メガデモを作ることにしました。十進演算ランタイムが圧縮できねえよ。http://ameblo.jp/mratassyu/手提鞄屋魚有店(てさげかばんやうおありてん)レスがついていないものを優先して、レスしています。時々、見当外れなレスをします。ミクプラは自由にリンクできます。