index.html

Clear[y, x] &nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&n ... &nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;(* 変数yとxの中身を初期化 *)

A[t_] = -2 + 10 t * UnitStep[0.3 - t] - x[t] ; &nbsp;&nbsp;&nbsp;&nbsp;(* t = 0.3 で変化する微分係数の項 *)

B[t_] = -2 + 10 t * UnitStep[0.4 - t] - y[t] ; &nbsp;&nbsp;&nbsp;&nbsp;(* t = 0.4 で変化する微分係数の項 *)

eq1 = y

eq2 = x

cond1 = y[0] == 1 ; &nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; ... ;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;(* 初期状態 *)
cond2 = x[0] == 1 ;

sol1 = NDSolve[{eq1, eq2, cond1, cond2}, {y, x}, {t, 0, 3}] ; &nbsp;&nbsp;&nbsp;(* 微分方程式の数値解析 *)

y[t_] = First[y[t]/.sol1] ; &nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbs ... nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;(* 計算結果の代入 *)

x[t_] = First[x[t]/.sol1] ; 
&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nb ... ;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;(* グラフ描画 *)

[Graphics:HTMLFiles/index_12.gif]

[Graphics:HTMLFiles/index_13.gif]