In[]:=
Needs["VisualDSolve`VisualDSolve`"]
PhasePlot[{Derivative[1][x][t]==y[t],Derivative[1][y][t]==–x[t]–y[t]/20},{x[t],y[t]},{t,0,5π},{x,–1,1},{y,–1,1},InitialValues{1/2,1/2},ShowInitialValuesTrue,FrameLabel{x,y},RotateLabelFalse]
Out[]=
Leave a Reply