1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
|
size(8cm);
import smoothcontour3;
// Erdos lemniscate of order n:
real erdos(pair z, int n) { return abs(z^n-1)^2 - 1; }
real h = 0.12;
// Erdos lemniscate of order 3:
real lemn3(real x, real y) { return erdos((x,y), 3); }
// "Inflate" the order 3 (planar) lemniscate into a
// smooth surface:
real f(real x, real y, real z) {
return lemn3(x,y)^2 + (16*abs((x,y))^4 + 1) * (z^2 - h^2);
}
// Draw the implicit surface on a box with diagonally opposite
// corners at (-3,-3,-3), (3,3,3).
draw(implicitsurface(f,a=(-3,-3,-3),b=(3,3,3),overlapedges=true),
surfacepen=material(diffusepen=gray(0.5),emissivepen=gray(0.4),
specularpen=gray(0.1)));
|
