e=0.6;
for (i = [0:1.0: N]){
theta=i * 360.0/N;
R=a1*sqrt(1-e*e*sin(theta)*sin(theta));
r1=0.1*R;
rotate( theta, [0, 0, 1])
translate([r1,0, w*i])
cube([R,0.1,W], center = false);
}
a=25.0;
M=8;
e2=0.99;
for(j=[0:1:M]){
phi=j*180/M;
translate([0,0,40.0*j])
rotate(phi,[0,0,1])
for (i = [0:1.0: N]){
theta=i * 360.0/N;
R=5*a/j*sqrt(1-e2*e2*sin(theta)*sin(theta));
r1=0.1*R;
rotate( theta, [0, 0, 1])
translate([r1,0, w*i])
cube([R,0.1,W], center = false);
}
}