The coneplot function is a velocity plot.
If you want to plot a single cone, this will get you started:
N = 10;
r = linspace(0, 1, N);
[X,Y,Z] = cylinder(r, N);
figure(1)
h = surf(X, Y, Z);
axis square
rotate(h, [-1 1 0], 90);
I don’t know exactly what you want to do, so you will have to experiment with this. Other functions to consider are cart2sph and sph2cart to produce the rotation you want.
