NDSolve[{y'[x]==Sin[y[x]],y[0]==1}, y, {x,0,4}] Plot[Evaluate[y[x] /. %], {x,0,4}]
(この NDSolve[] は数値的に微分方程式を解くコマンドであり、結果 は離散的な点での関数値を近似的に求めたもののリストである。 計算していない点での値は補間により簡略計算する。 こういうものを InterpolatingFunction という。)
y[1.5] /. %%
