パラメトリック微分方程式 

パラメトリックの初期条件．

In[1]:=

X pfun = ParametricNDSolveValue[{y''[t] + 2 y[t] == Sin[.4 t] + Cos[y[t]], y[0] == a, y'[0] == b}, y, {t, 0, 10}, {a, b}]; Grid[Plot[Evaluate@Table[pfun[a, 1][t], {a, -2, 2, .5}], {t, 0, 10}, PlotRange -> All], Plot[Evaluate@Table[pfun[1, b][t], {b, -2, 2, .5}], {t, 0, 10}, PlotRange -> All]]

Out[1]=