int samples = 20;
float a1 = 10., a2 = 10., a3 = 10.;
float u1 = 0., u2 = 20., v1 = 0., v2 = 20.;
float dU = (u2 - u1) / samples;
float dV = (v2 - v1) / samples;
float n = 1., e = 1.;
void setup(){
  size(500,500, P3D); //setup the screen
  camera(-15,15,-20,0,0,0,0,0,1); //get a viewpoint
  float u = u1;
  for(int i=0; i<samples; i++){
    float v = v1;
    for(int j=0; j<samples; j++){
      float x = a1 * sqCos (u, n) * sqCos (v, e);
      float y = a2 * sqCos (u, n) * sqSin (v, e);
      float z = a3 * sqSin (u, n);
      point(x,y,z);
      v += dV;
    }
    u += dU;
  }
}

float sign ( float x ) {
  if ( x < 0 )return -1;
  if ( x > 0 )return 1;
  return 0;
}
float sqSin( float v, float n ) {
  return sign(sin(v)) * pow(abs(sin(v)),n);
}
float sqCos( float v, float n ) {
  return sign(cos(v)) * pow(abs(cos(v)),n);
}