plot - 在 Maxima 中绘制 ODE 的解（“表达式 evalues 无处不在”）

Question

我正在w使用 Maxima 求解二阶 ODE，它给了我一些答案，但我不知道如何绘制表达式。会话如下所示：

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);

当我尝试solNum用

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

我正进入（状态

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

我是 Maxima 初学者，如果我忽略了一些明显的事情，我深表歉意。我看过的教程只显示plot2d了一个表达式。

Answer 1

只需绘制解决方案的 RHS 即可。

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]);