# I obtained different figure after using fsurf and surf

1 view (last 30 days)
Casey lv on 27 Apr 2020
Commented: Toder on 29 Apr 2020
Here is the code, I think these two commands will generate the same figure but it did not.
% syms x t;
k = 0.6;
sigma=-0.5;
r =0.5;
s = 1.2;
x=10;t=0;
lambda=k.*sqrt(2.*r.^2./-sigma);
A=lambda.*((k.*jacobiSN(r.*x-s.*t,k.^2).*jacobiCN(r.*x-s.*t,k.^2))./(jacobiDN(r.*x-s.*t,k.^2)));
AA=lambda.*((k.*jacobiSN(r.*x-s.*t,k.^2).*jacobiCN(r.*x-s.*t,k.^2))./(jacobiDN(r.*x-s.*t,k.^2)));
B=4.*r.*s.*ellipticE(k.^2).*ellipticK(k.^2).^(-1)+(-2).*r.*s.*jacobiDN(s.*t+(-1).*r.*x,k.^2).^(-2).*(1+(-1).*k.^2+jacobiDN(s.*t+(-1).*r.*x,k.^2).^4);
A;
B;
lambda;
[x,t]=meshgrid(-5:0.2:5);
bb=k.*jacobiSN(r.*x-s.*t,k.^2).*jacobiCN(r.*x-s.*t,k.^2)./(jacobiDN(r.*x-s.*t,k.^2));
figure;
mesh(x,t,abs(bb));
syms c d;
bb=k.*jacobiSN(r.*c-s.*d,d.^2).*jacobiCN(r.*c-s.*d,k.^2).*(jacobiDN(r.*c-s.*d,k.^2)).^(-1);
figure;
fsurf(abs(bb),'MeshDensity',10);
g=400-300+500;
g ##### 1 CommentShowHide None
Ameer Hamza on 27 Apr 2020
In your code, the 2nd equation is reciprocal of the first. Shouldn't that change the output.

Toder on 28 Apr 2020
It appears you may have a typo in your second definition of bb. Your first call to jacobiSN in this line currently reads as
jacobiSN(r.*c-s.*d,d.^2)
Based on your first definition of bb, it looks like it should instead be
jacobiSN(r.*c-s.*d,k.^2)
##### 2 CommentsShowHide 1 older comment
Toder on 29 Apr 2020
Sorry, I don't know much about Mathematica. My only recommendation is to check that the axes are the same; if they're different than the surf in Matlab the plot could look very different.

### Community Treasure Hunt

Find the treasures in MATLAB Central and discover how the community can help you!

Start Hunting!