Creating a Matlab Contour Plot with Iterative Growth - Help Request

[member 428835]

hey pf!

i am trying to plot a contour (topographic) graph the following $f(x,y) = \frac{1}{x^2+y^2} e^{x-(x^2+y^2)^{.5}}$ where $f(x,y) = k$ where $k = .001$ and doubles in size for about 8 iterations. my attempt is here:

x=-3:0.25:3;
y=-3:0.25:3; z=.25:.25:1;
[X,Y]=meshgrid(x,y);
Z=1./(X.^2+Y.^2).*exp(X-(X.^2+Y.^2));
surfc(X,Y,Z); xlabel('x'); ylabel('y'); zlabel('z')

can someone please help me out?

thanks!

 

[kreil]

A solid first shot. I would probably use contour3 instead of surfc (unless you need the surface more than the contour lines). Also, you don't use k at all in your code even though you're trying to solve f(x,y) = k for several values of k.

Try this out:

k = 1e-3;
figure;
hold on
[X,Y]=meshgrid(-3:0.1:3);
for n = 1:8
    k = k*2^(n-1);
    Z = 1./(X.^2+Y.^2).*exp(X-sqrt(X.^2+Y.^2))-k;
    contour3(X,Y,Z)
end

You get a plot like the one attached. This might need more tweaking but should be closer to what you're looking for.

 

[member 428835]

Hey thanks! I really appreciate it!