function superhyper(epsilon,a)
  n=50;
  d=.1;
  etamax=pi/2-d;
  etamin=-pi/2+d;
  wmax=3*pi/2-d;
  wmin=pi/2+d;
  deta=(etamax-etamin)/n;
  dw=(wmax-wmin)/n;
  k=0;
  l=0;
  for i=1:n+1
    eta(i)=etamin+(i-1)*deta;
    for j=1:n+1
      w(j)=wmin+(j-1)*dw;
      x(i,j)=a(1)*sign(sec(eta(i)))*abs(sec(eta(i)))^epsilon(1)*sign(sec(w(j)))*abs(sec(w(j)))^epsilon(2);
      y(i,j)=a(2)*sign(sec(eta(i)))*abs(sec(eta(i)))^epsilon(1)*sign(tan(w(j)))*abs(tan(w(j)))^epsilon(2);
      z(i,j)=a(3)*sign(tan(eta(i)))*abs(tan(eta(i)))^epsilon(1);
    endfor;
  endfor;
   mesh(x,y,z);

endfunction;