mathematica 画微分方程数值解

tt1=NDSolve[{D[y[t],t]\[Equal]y[t]^2+t^2,y[0]\[Equal]0},y,{t,-2,2}]
{{y -> InterpolatingFunction[{{-2., 2.}}, "<>"]}}
Plot[y[t] /. tt1, {t, -2, 2}]

该方程的解在多处发散。

但你很幸运，方程有解析解，为
(t BesselJ[3/4, t^2/2])/BesselJ[-(1/4), t^2/2]

其中BesselJ表示第一类贝塞尔函数

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
如果你要数值求解的话，t的范围是
{-2.003147359426885`， 2.003147359426885`}

f = NDSolve[{y'[t] == y[t]^2 + t^2, y[0] == 0}, y, {t, -20, 20}];
Plot[Evaluate[y[t] /. f[[1]]], {t, -20, 20}]

你如果不隐身我还可以给你答案的.但是......