I obtained different figure after using fsurf and surf

1 view (last 30 days)
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 Comment
Ameer Hamza
Ameer Hamza on 27 Apr 2020
In your code, the 2nd equation is reciprocal of the first. Shouldn't that change the output.

Sign in to comment.

Accepted Answer

Toder
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 Comments
Toder
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.

Sign in to comment.

More Answers (0)

Products

Community Treasure Hunt

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

Start Hunting!