1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
|
/* -*-ePiX-*- */
/* levelset.xp -- level sets */
#include "epix.h"
using namespace ePiX;
const int MAX=3;
const double k=0.5;
P F1(double u, double v)
{
return P(u+v/u, u-v/u, 2*k*v);
}
double f(double x, double y)
{
return 0.5*k*(x*x-y*y);
}
int main()
{
bounding_box(P(-MAX,-MAX), P(MAX, MAX));
unitlength("0.5in");
picture(4,4);
// offset(4,0);
begin();
domain R1(P(-MAX,-MAX), P(MAX,MAX), mesh(4,4), mesh(60,60));
domain R2(P(-4,-4), P(4,4), mesh(1,24), mesh(120,120));
viewpoint(1,-3,5); camera.range(40);
clip_box(P(MAX,MAX,2*MAX));
clip();
red();
plot(f, R1);
bold(); rgb(0.5,0.5,1);
line(P(MAX,-MAX,0), P(-MAX, MAX,0));
plot(F1, R2);
end();
}
|
