summaryrefslogtreecommitdiff
path: root/Master/texmf-dist/doc/metapost/featpost/example/standard/hiddensurfaces.mp
blob: df698f038bf5f34e87585eae2cc0354d288c4d7a (plain)
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
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
% hiddensurfaces.mp
% L. Nobre G.
% 2003

%input featpost3Dplus2D;

def zsurface( expr xc, yc ) = 
	begingroup
		numeric zout;
		zout = yc**4-xc*yc*(1-yc**2);
		if zout > 1.5:
			zout := 1.5;
		fi;
		( zout )
	endgroup 
enddef;

f := (2.5,2,1.25);
Spread := 127;

beginfig(1); 
	numeric i, j, nx, ny, lowx, higx, lowy, higy, refz
	        posx, posy, posz, facz, counter;
	color poi[][];
	nx = 30;
	ny = 30;
	lowy = -2.5;
	higy = 1;
	lowx = 1.5;
	higx = 3;
	refz = 0;
	facz = 5;
	counter = 0;
	for i=0 upto nx:
	    for j=0 upto ny:
		posx := lowx + i*( higx - lowx )/nx;
		posy := lowy + j*( higy - lowy )/ny;
		posz := zsurface( posx, posy );
		posx := bracket(lowx,posx,higx);
		posy := bracket(lowy,posy,higy);
		posz := posz/facz + refz;
		poi[i][j] := ( posx, posy, posz );
	    endfor;
	endfor;
	for i=1 upto nx:
	    for j=1 upto ny:
		counter := incr(counter);
		npf[counter] := 3;
		F[counter]p1 := poi[i-1][j-1];
		F[counter]p2 := poi[i][j-1];
		F[counter]p3 := poi[i][j];
		counter := incr(counter);
		npf[counter] := 3;
		F[counter]p1 := poi[i-1][j-1];
		F[counter]p2 := poi[i-1][j];
		F[counter]p3 := poi[i][j];
	    endfor;
	endfor;
	NF := counter;
	drawarrow rp((-1,-1,refz))--rp((1,-1,refz));
	drawarrow rp((-1,-1,refz))--rp((-1,1,refz));
	fill_faces( withcolor blue );
endfig;

end;