Plotting solution of ODE in Maxima ("expression evalues to non-numeric value everywhere")

StackOverflow https://stackoverflow.com/questions/15619231

  •  29-03-2022
  •  | 
  •  

Question

I am solving for w in a 2nd order ODE using Maxima, it gives me some answer, but I am at loss as to how to plot the expression. The session looks like this:

eq: -E*(I0+I1*x)*'diff(w,x,2)=(f1/6)*x^3+(f0/2)*x^2+l*(f0/2+2*f1*l^2/3)*x;
sol2: ode2(eq,w,x);
solbc2: bc2(sol2,x=0,w=0,x=l,w=0);
solNum: subst([E=1,I0=1,I1=1,f0=1,f1=1,l=1],solbc2);

When I try to plot solNum with

plot(solNum,[x,0,1]);

I am getting

plot2d: expression evaluates to non-numeric value everywhere in plotting   range.
plot2d: nothing to plot.

I am Maxima beginner, I apologize if I overlooked something obvious. The tutorials I've crossed only show plot2d with an expression in it.

Was it helpful?

Solution

Just plot RHS of the solution.

eq: -E*(I0+I1*x)*'diff(w,x,2)=(f1/6)*x^3+(f0/2)*x^2+l*(f0/2+2*f1*l^2/3)*x;
sol2: ode2(eq,w,x);
solbc2: bc2(sol2,x=0,w=0,x=l,w=0);
solNum: subst([E=1,I0=1,I1=1,f0=1,f1=1,l=1],solbc2);

plot2d(rhs(solNum),[x,0,1]);
Licensed under: CC-BY-SA with attribution
Not affiliated with StackOverflow
scroll top