2013-03-25 35 views
1

我解决w在二阶ODE使用极大值,它给了我一些答案,但我在如何绘图表达方式。这次会议是这样的:在Maxima中绘制ODE的解决方案(“表达evalues到非数值无处不在”)

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. 

我是千里马初学者,我很抱歉,如果我忽略了一些东西明显。我穿过的教程只显示plot2d并带有一个表达式。

回答

1

只是绘制rhs 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]);