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
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
|
% D. Roegel, 7 January 2001
% roegel@loria.fr
verbatimtex
%&latex
\documentclass{article}
\usepackage{mathpple}
\begin{document}
etex
input 3danim
input 3dgeom
output_res:=72;
drawing_scale:=13cm;
numeric b,c,d,e;b=70;c=0;d=20;e=50; % parameterized angles
vardef def_T(expr inst)=
new_obj_points(inst,10); % 10 points
set_T_points(inst);
enddef;
vardef set_T_points(expr inst)=
set_point(1)(0,0,0);
set_point(2)(3*cosd(b),3*sind(b),0);
set_point(3)(6*cosd(c),6*sind(c),0);
new_vec(v_a);new_vec(v_b);
vec_def_vec_(v_a,vec_I);
vec_rotate_(v_a,vec_K,d);
vec_prod_(v_b,v_a,vec_K);
vec_rotate_(v_a,v_b,e);
vec_mult_(v_a,v_a,4.5);
vec_sum_(pnt(4),pnt(1),v_a);
free_vec(v_b);free_vec(v_a);
% Determination of I and J:
% I=A + unit(AB)
vec_diff(5,2,1);
vec_unit(5,5);
vec_sum(5,5,1);
% J=A + 4*unit(AC)
vec_diff(6,3,1);
vec_unit(6,6);
vec_mult(6,6,4);
vec_sum(6,6,1);
new_plane(bcd)(2,3,4);
new_line(ij)(5,6);
boolean b;
b:=def_inter_p_l_pl(7)(ij)(bcd); if not b: ... fi;
% K=middle of AD
mid_point(8,1,4);
set_line(ij)(5,8);
b:=def_inter_p_l_pl(9)(ij)(bcd); if not b: ... fi;
set_line(ij)(6,8);
b:=def_inter_p_l_pl(10)(ij)(bcd); if not b: ... fi;
free_line(ij);
free_plane(bcd);
enddef;
vardef draw_T(expr inst)=
drawoptions(dashed evenly);draw_line(2,4);draw_line(8,5);
drawoptions();
draw_line(1,2);draw_line(1,3);draw_line(1,4);
draw_line(2,3);draw_line(3,4);
pickup pencircle scaled 2pt;
draw_line(5,5);draw_line(6,6);draw_line(7,7);draw_line(8,8);
draw_line(9,9);draw_line(10,10);
pickup pencircle scaled .4pt;
draw_line(5,7);draw_line(3,7);
draw_line(8,9);draw_line(4,9);
draw_line(8,10);draw_line(3,10);
draw_line_extra(9,10)(-0.1,1.1);
label_obj.top(btex $A_1$ etex,1);
label_obj.lft(btex $B_2$ etex,2);
label_obj.llft(btex $C_3$ etex,3) shifted (-2.5mm,1mm);
label_obj.top(btex $D_4$ etex,4);
label_obj.lft(btex $I_5$ etex,5);
label_obj.lft(btex $J_6$ etex,6) shifted (-2mm,0);
label_obj.lrt(btex $M_7$ etex,7);
label_obj.urt(btex $K_8$ etex,8);
label_obj.lrt(btex $N_9$ etex,9);
label_obj.lrt(btex $L_{10}$ etex,10);
enddef;
assign_obj("tetra","T");
new_vec(v_a);new_vec(v_b);new_vec(v_c);
vec_diff_(v_a,pnt_obj("tetra",3),pnt_obj("tetra",4)); % vec(DC)
vec_mult_(v_a,v_a,5);
vec_diff_(v_b,pnt_obj("tetra",3),pnt_obj("tetra",2)); % vec(BC)
vec_mult_(v_b,v_b,5);
vec_diff_(v_c,pnt_obj("tetra",1),pnt_obj("tetra",3)); % vec(CA)
vec_mult_(v_c,v_c,3);
vec_sum_(Obs,pnt_obj("tetra",3),v_a); % C + 5DC
vec_sum_(Obs,Obs,v_b);% C + 5DC + 5BC
vec_sum_(Obs,Obs,v_c);% C + 5DC + 5BC + 3CA
free_vec(v_c);free_vec(v_b);free_vec(v_a);
%for i:=0 upto 20:
i=8;
beginfig(200+i);
% Observator
%set_point_(Obs)(20*cosd(3.6*i),20*sind(3.6*i),6);
Obs_phi:=0;Obs_dist:=2;point_of_view_obj("tetra",1,Obs_phi);
draw_obj("tetra");
endfig;
%endfor;
end.
|
