I am solving a set of differential equations using the following command, numericalsol:=dsolve({ODEsys, ics}, numeric, method=rkf45) and I'm plotting my numerical solution using the following command, plots[odeplot](numericalsol,[t,w[1](t)],0..20); But when the solution is plotted the graph looks jagged and not smooth at all. I'm assuming that I am not getting enough points in my numerical solution and that's whats causing the problem. Is there any way to get more points so that my plot will look smooth?

I need help making a 3D plot of a numerical solution to Euler's Equation for Rigid Bodies(Non-linear Differential Equation). The solution to this equation are the angular velocity components w1,w2,w3. I am able to plot each component using the following commands, plots[odeplot](numericalsol,[t,w[1](t)],0..20); plots[odeplot](numericalsol,[t,w[2](t)],0..20); plots[odeplot](numericalsol,[t,w[3](t)],0..20); However, how do I make a 3D plot using these numerical solutions?? Any help is greatly appreciated.

I need to numerically solve a 2nd order ode in maple and compare the numerically solved DE with the analytical solution by finding the error and plotting it as a function of time. I know how to solve the DE numerically in maple using dsolve(......, numeric), for example, sol :=dsolve({D(y)(t)=t^2-sin(y(t)),y(0)=1},y(t),type=numeric); and then out comes the solution as a procedure, sol := proc(rkf45_t) ... end sol(2); [t=2, y(t)=2.156829]