How to plot this function

2 views (last 30 days)
Dinh Le Dung
Dinh Le Dung on 1 Aug 2022
Commented: Star Strider on 2 Aug 2022

Accepted Answer

Star Strider
Star Strider on 2 Aug 2022
Force the NaN value at (0,0) to be 0
[X,Y] = ndgrid(linspace(-5,5,50));
Z = (X.^3.*Y - X.*Y.^3)./(X.^2+Y.^2);
Z(isnan(Z)) = 0;
figure
surfc(X,Y,Z)
colormap(turbo)
.
  2 Comments
Star Strider
Star Strider on 2 Aug 2022
As always, my pleasure!

Sign in to comment.

More Answers (2)

Sam Chak
Sam Chak on 1 Aug 2022
Think it should look like this:
[X, Y] = meshgrid(-0.5:1/40:0.5);
Z = X.*Y.*(X.^2 - Y.^2)./(X.^2 + Y.^2);
surf(X, Y, Z), xlabel('x'), ylabel('y'), zlabel('f(x, y)')
  3 Comments
Dinh Le Dung
Dinh Le Dung on 2 Aug 2022
Oh it look great!! Thanks a lot!!

Sign in to comment.


Issa
Issa on 1 Aug 2022
Hi!
Since you need to multiply by the truth, maybe this below:
f = @(x,y)((((x.^3).*y - x.*(y.^3)) ./ (x.^2 + y.^2) ).*( x ~= 0 & y~=0) + 0.*( x == 0 & y == 0))
f = function_handle with value:
@(x,y)((((x.^3).*y-x.*(y.^3))./(x.^2+y.^2)).*(x~=0&y~=0)+0.*(x==0&y==0))
[X, Y] = meshgrid([1:.5:10], [0:.1:10]);
F = f(X, Y) ;
surf(X, Y, F)
colorbar
Hope this helps!
  1 Comment
Dinh Le Dung
Dinh Le Dung on 2 Aug 2022
I try to zoom in the interval [-0.5, 0.5] but it seems like it does not cover the case f(0,0) = 0

Sign in to comment.

Community Treasure Hunt

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

Start Hunting!

Translated by