help debugging 3-D plot code

3 views (last 30 days)
I'm trying to figure out why this ellipsoid plot isn't working the way it's supposed to.
It seems that the grid direction may be wrong, but every way that I've tried to change it changes the plot in other undesirable ways as well.
I'm not allowed to use the ellipsoid function, I have to use meshgrid, its very close to correct, but isn't.
The code thus far is.
clear variables
clc
[phi,theta] = meshgrid(linspace(0,2*pi,20),(linspace(0,pi,20)));
figure(1)
a=1;
b=.5;
x=a*cos(phi).*cos(theta);
y=b*cos(phi).*sin(theta);
z=b*sin(phi);
surf(x,y,z);
hold on
AZ=-13;
EL=28;
view(AZ,EL)
xlabel('x')
ylabel('y')
zlabel('z')
title('Ellipsoid')
set(gca, 'XTick',[-0.8,-0.6,-0.4,-0.2,0,0.2,0.4,0.6,0.8])
set(gca, 'YTick',[-0.4,-0.2,0,0.2,0.4])
when run the code produces
what I need is
  3 Comments
Brendan Clark
Brendan Clark on 16 Apr 2021
ope, good catch
its supposed to be
x=a*cos(φ)sin(θ);
y=b*cos(φ)sin(θ);
z=b*sin(φ);
when i corrected the cos to sin that was incorrect I ended up with a 2d graph on the 3d plan, what were the corrections that you made?
Clayton Gotberg
Clayton Gotberg on 17 Apr 2021
Check your new equations again! Now, x and y are defined the same.

Sign in to comment.

Accepted Answer

Bjorn Gustavsson
Bjorn Gustavsson on 17 Apr 2021
You have to get your equations for x, y and z right:
x = a*cos(phi).*sin(theta);
y = b*sin(phi).*sin(theta);
z = b*cos(theta);
surf(x,y,z);
HTH
  2 Comments
Brendan Clark
Brendan Clark on 18 Apr 2021
Thank you. This would be the third typo I've found in my homework assignment for this week. What I posted for the equations with x and y being the same was copy and pasted from my hw. The right equation did the trick.
Bjorn Gustavsson
Bjorn Gustavsson on 18 Apr 2021
When I have these sorts of typo-problems I find that printing the code out on paper and annotate it by hand with circling in, crossing over and connecting variables and statements to be helpful. If you see what you think that you've written and not what you've actually written then that helps - sometimes.

Sign in to comment.

More Answers (0)

Community Treasure Hunt

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

Start Hunting!