NDSolve[]

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] /. %%