Trial>> syms s;
[t,z]=meshgrid(-10:0.5:10);
arfa=sin(s);beita=sin(z);
rou=sqrt(arfa./beita);
eita=0.5;v=0;miu=1;w0=0;z0=0;c1=1;
w=c1.*int(arfa,s,0,s)+w0;
w1=subs(w,s,z);
Z0=0.5.*arfa;
Z=int(Z0,s,0,s)+z0;
Z1=subs(Z,s,z);
T=t+w1;
A=2.*eita.*(T-miu);
B=4.*eita.^2.*(T.^2+16.*eita.^2.*Z1.^2);
u=4.*eita.*rou.*((8.*eita.^2.*Z1-i).*cosh(A)+2.*i.*eita.*T.*sinh(A))./(cosh(A).^2+B);
u=abs(double(u));
mesh(t,z,u)
错误使用 symengine
DOUBLE cannot convert the input expression into a double array.
出错 sym/double (line 643)
Xstr = mupadmex('symobj::double', S.s, 0);