数値的に解くための代表的な離散化方法の一つとして有限要素法がある。有限要素法では、解析モデルの物体や空間にメッシュを生成して有限要素に分割し、有限要素の中での変数分布として一次式等の単純な変数分布を仮定することにより、問題を線形方程式の求解に帰着させる。この線形方程式の近似解を行列計算により求め、解析モデルの複雑な変数分布を数値シミュレーションすることができる。データを入力する入力部と、入力部に入力されるデータに基づいて有限要素法による解析を行う解析部と、解析部による解析結果を目的に応じた形式で可視化するなどして出力する出力部と、から構成される。シミュレーションシステムにおける計算コストの大部分は、線形方程式の求解に係る行列計算に費やされ、解析モデルの規模の増大に伴って行列計算の計算コストは指数的に増大する。そのため、行列計算を高速化するための種々の技術が開発されている。非特許文献1には、解に近づくベクトルを逐次計算する反復法による非対称行列の解法が記載されている。非特許文献2には、反復法によるスパース行列の解法の収束性を改善する技術が記載されている。加えて、近年の計算機性能の向上や、マウスを用いた直感的な形状入力が可能なGUIやメッシュの自動生成機能のユーザビリティの向上により、有限要素法によるシミュレーションシステムのユーザは、専用の大型計算機や熟練したオペレータを有する大企業や研究所等に限定されないものとなりつつある。解析プログラムに関する。測ベクトルとして利用する第1のレベルと、(2)前記再接続メッシュ生成工程が生成したメッシュと、当該形状変更前の解析モデルに対する解析において用いた係数行列Aを反復法で解くための前処理行列Bの非零要素の場所を表す非零パターン行列Bpと、該前処理行列Bの非零要素の値を表す値リスト行列Bvと、を再利用するとともに、当該形状変更前の解析モデルに対する解析において求めた解を反復法の初期推測ベクトルとして利用する第2のレベルと、(3)当該形状変更前の解析モデルに対する解析において用いたメッシュ、行列及び解を再利用しない第3のレベルと、の3つの再利用のレベルのうちから一のレベルを選択することにより、前記再利用する部分を決定することを特徴とする有限要素法解析プログラム。