clear all;close all;X=-10:1:10;T=-10:1:10;C1=2;C2=1;b=-1.5;[x,t]=meshgrid(X,T);x1=C1*(C1*exp(sqrt(2)*b*x)+C2*exp(-sqrt(2)*b*x)).*sin(1/2*sqrt(-1i.^4*b.^2+4)*t).*(C1*exp(sqrt(2)*b*x)+C2*exp(-sqrt(2)*b*x))+C2*(C1*exp(sqrt(2)*b*x)+C2*exp(-sqrt(2)*b*x)).*cos(1/2*sqrt(-1i.^4*b.^2+4)*t).*(C1*exp(sqrt(2)*b*x)+C2*exp(-sqrt(2)*b*x));y1=-1/2*C1*(C1*sqrt(2)*exp(sqrt(2)*b*x)-C2*sqrt(2)*exp(-sqrt(2)*b*x)-C1*exp(sqrt(2)*b*x)-C2*exp(-sqrt(2)*b*x)).*cos(1/2*sqrt(-1i.^4*b^2+4)*t).*(C1*sqrt(2)*exp(sqrt(2)*b*x)-C2*sqrt(2)*exp(-sqrt(2)*b*x)-C1*exp(sqrt(2)*b*x)-C2*exp(-sqrt(2)*b*x))*sqrt(-1i*(C1*sqrt(2)*exp(sqrt(2)*b*x)-C2*sqrt(2)*exp(-sqrt(2)*b*x)-C1*exp(sqrt(2)*b*x)-C2*exp(-sqrt(2)*b*x)).^4*b(C1*sqrt(2)*exp(sqrt(2)*b*x)-C2*sqrt(2)*exp(-sqrt(2)*b*x)-C1*exp(sqrt(2)*b*x)-C2*exp(-sqrt(2)*b*x)).^2+4)+1/2*C2*(C1*sqrt(2)*exp(sqrt(2)*b*x)-C2*sqrt(2)*exp(-sqrt(2)*b*x)-C1*exp(sqrt(2)*b*x)-C2*exp(-sqrt(2)*b*x)).*sin(1/2*sqrt(-1i.^4*b^2+4)*t).*(C1*sqrt(2)*exp(sqrt(2)*b*x)-C2*sqrt(2)*exp(-sqrt(2)*b*x)-C1*exp(sqrt(2)*b*x)-C2*exp(-sqrt(2)*b*x))*sqrt(-1i.*(C1*sqrt(2)*exp(sqrt(2)*b*x)-C2*sqrt(2)*exp(-sqrt(2)*b*x)-C1*exp(sqrt(2)*b*x)-C2*exp(-sqrt(2)*b*x)).^4*b*(C1*sqrt(2)*exp(sqrt(2)*b*x)-C2*sqrt(2)*exp(-sqrt(2)*b*x)-C1*exp(sqrt(2)*b*x)-C2*exp(-sqrt(2)*b*x)).^2+4)-1/2*1i*(C1*sqrt(2)*exp(sqrt(2)*b*x)-C2*sqrt(2)*exp(-sqrt(2)*b*x)-C1*exp(sqrt(2)*b*x)-C2*exp(-sqrt(2)*b*x)).^2*b*(C1*sqrt(2)*exp(sqrt(2)*b*x)-C2*sqrt(2)*exp(-sqrt(2)*b*x)-C1*exp(sqrt(2)*b*x)-C2*exp(-sqrt(2)*b*x))*C1*(C1*sqrt(2)*exp(sqrt(2)*b*x)-C2*sqrt(2)*exp(-sqrt(2)*b*x)-C1*exp(sqrt(2)*b*x)-C2*exp(-sqrt(2)*b*x))*sin(1/2*sqrt(-1i.^4*b^2+4)*t).*(C1*sqrt(2)*exp(sqrt(2)*b*x)-C2*sqrt(2)*exp(-sqrt(2)*b*x)-C1*exp(sqrt(2)*b*x)-C2*exp(-sqrt(2)*b*x))-1/2*1i*(C1*sqrt(2)*exp(sqrt(2)*b*x)-C2*sqrt(2)*exp(-sqrt(2)*b*x)-C1*exp(sqrt(2)*b*x)-C2*exp(-sqrt(2)*b*x)).^2*b*(C1*sqrt(2)*exp(sqrt(2)*b*x)-C2*sqrt(2)*exp(-sqrt(2)*b*x)-C1*exp(sqrt(2)*b*x)-C2*exp(-sqrt(2)*b*x))*C2*(C1*sqrt(2)*exp(sqrt(2)*b*x)-C2*sqrt(2)*exp(-sqrt(2)*b*x)-C1*exp(sqrt(2)*b*x)-C2*exp(-sqrt(2)*b*x))*cos(1/2*sqrt(-1i.^4*b^2+4)*t).*(C1*sqrt(2)*exp(sqrt(2)*b*x)-C2*sqrt(2)*exp(-sqrt(2)*b*x)-C1*exp(sqrt(2)*b*x)-C2*exp(-sqrt(2)*b*x));A=-2*1i*b*x1.*conj(y1);B=x1.*conj(x1)+y1.*conj(y1);R=(A./B);surf(x,t,abs(R))
When I evaluate it, error occur with comment "Subscript indices must either be real positive integers or logicals".
How to remove it?
Best Answer