このディレクトリの索引
#  
#  常微分方程式 dy/dx = f'(x,y) = y + sin(x)
#  
#  初期値 (x0,y0) = (0,0.5) の時の解は f(x,y) = exp(x) - (sin(x) + cos(x)) * 0.5 
#  


f(_x,_y,_z) :-
        _y is exp(_x) - ( sin(_x) + cos(_x)) * 0.5).

'オイラー法'(_f,_x,_y,_h,X) :-
        f(_x,_y,_z),
        X is _y + _z * _h.

'オイラー法を用いて常微分方程式の解を数値的に求める'(_分割数,_xの初期値,_xの最終値,_yの初期値,_x,_y) :-
        _刻み is ( _xの最終値 - _xの初期値) / _分割数,
        繰り返し(0,_分割数,_刻み,_xの初期値,_yの初期値,_x,_y).

'オイラー法を用いて常微分方程式の解を数値的に求める'(_i,_分割数,_刻み,_x_1,_y,_x,_y) :-
        _i =< _分割数,
        _x is _x_1 + _i * _刻み.
'オイラー法を用いて常微分方程式の解を数値的に求める'(_i,_分割数,_刻み,_x_1,_y_1,_x,_y) :-
        _i < _分割数,
        _x_2 is _x_1 + _i * _刻み,
        'オイラー法'(_x_2,_y_1,_刻み,_y_2),
        _i_2 is _i + 1,
        'オイラー法を用いて常微分方程式の解を数値的に求める'(_i_2,_分割数,_刻み,_x_2,_y_2,_x,_y).