■例題の計算■
ここで取り上げている問題では、三角形要素で用いた時と同じ計算領域と境界条件を使います。三角形要素との精度の比較ができますからね。念のためにもう一度、領域と境界条件を下に示します。

それでは早速、プログラムを実行してみましょう。その前に、入力データを作成する必要がありますので、 SETFEM4.FOR(コンパイル忘れずに)を走らせてデータファイル FEM4DATA.QQQ を作成して下さい。
プログラムの実行中、NEXとNEYの入力を要求してきますので、共に2をキーインして下さい。NEXとNEYは、それぞれ、x軸方向とy軸方向の要素の数を意味します。NEX=2とNEY=2とすると、下図に示す様に、領域は4つの要素に分割されます。

その後、FEM4Q(コンパイル忘れずに)とキーインすると、 SOLUTION.FEM と言う名前の解析結果ファイルが出来上がります。このファイルには、入力データのエコーと節点での温度が計算結果として出力されています。