2

Maximaを使用して2次常微分方程式で解きますw。それは私にいくつかの答えを与えますが、式をプロットする方法について私は途方に暮れています。セッションは次のようになります。

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、式が含まれているだけです。

4

1 に答える 1

2

ソリューションの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]);
于 2013-03-25T22:10:52.300 に答える