Discharge gsl_root_fsolve in IsotonicContraction

Replace gsl_root_fsolve with differential equations of motion in
IsotonicContraction module.
IsotonicContractionモデル中で、hSMLを代数方程式の求根で求めるのを廃止。代わりに運動方程式を微分方程式で解くようにした。バ
ネの質量に相当するパラメータKmは従来のモデルのシミュレーション結果から概算して 0.01 としているが、最適ではないかもしれない。