怎样用mathematic画常微分方程的图像

你是指mathematica吧
s = NDSolve[{y'[x] == y[x] Cos[x + y[x]],y[0] == 1},y,{x,0,30}](*此行求数值解*)
Plot[Evaluate[y[x] /.s],{x,0,30},PlotRange -> All](*此行画图*)
解析解
DSolve[{y''[x] + y[x] == 0,y[0] == 1,y'[0] == 1/3},y ,x] （*求解*）
Plot[y[x] /.%,{x,0,17}]（*画图*）
具体可参考
DSolve
和
NDSolve