v=linspace(0,2*pi,120); [u,v]=meshgrid(u,v); R=0.5; %set large radius of torus r=abs(cos(v/2))*0.2+0.05; %set small radius to vary with angle %parametrization...
