Thread Subject: plotting troubles

I'm trying to take a 2D plot and rotate it about an axis to make it a 3D plot. I have used

plot(x,((BESSELJ(0,x)).^2).*((sin(pi*x)).^2),'c+-');

and am trying to rotate this about the x-axis so that is makes a 3D image of it. I have tried a 3D approach using

x = (BesselJ(0,(x.^2+y.^2).^(1/2))).^2.*(sin(pi*t)).^2;
y = (BesselJ(0,(x.^2+y.^2).^(1/2))).^2.*(sin(pi*t)).^2;
z = t;
plot3(x, y, z);

but this gives me a bunch of screwy lines that does not at all look like what I want/need. Any help/suggestions?


~Josh

Hi,

try this :

%-------------
x = 0:.1:2;
y = (besselj(0,x).^2).*(sin(pi*x).^2);

th = linspace(0,2*pi,100);

yrev = cos(th(:))*y;
zrev = sin(th(:))*y;
xrev = repmat(x,numel(th),1);

figure
p=plot(x,y,'r-');
set(p,'linewidth',3)
hold on
surf(xrev,yrev,zrev)
grid on
axis equal vis3d
%-------------

Jerome