-
function superquadric(epsilon,a) n = 50;
etamax = pi/2;
etamin = -pi/2; wmax = pi; wmin = -pi; deta = (
etamax-etamin)/n; dw = (wmax-wmin)/n; [i,j] = meshgrid(1:n+1...
- supertorus:
function supertoroid(epsilon,a) n=50; d=.1;
etamax=pi; etamin=-pi; wmax=pi; wmin=-pi; deta=(
etamax-etamin)/n; dw=(wmax-wmin)/n; k=0; l=0; for i=1:n+1...