geo.tst
5.88 KB
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
114
115
116
117
118
119
120
121
122
123
124
125
126
127
128
129
130
131
132
133
0,0,0,1,0,[1e-10,1e-17],12,[1,50,0,25],0,0,0,"Warning: some commands like subs might change arguments order",
-10.0,10.0,-10.0,10.0,-10.0,10.0,-10.0,10.0,-10.0,10.0,-1.4,1.1,1,0.0,1.0,1,
"Done",
parameter([Ax,-10.0,10.0,-1,0]),
parameter([Ay,-10.0,10.0,-1,0]),
pnt(pnt[Ax+i*Ay,[536870968],A]),
parameter([Bx,-10.0,10.0,0.5,0]),
parameter([By,-10.0,10.0,1,0]),
pnt(pnt[Bx+i*By,[536870968],B]),
parameter([Cx,-10.0,10.0,1,0]),
parameter([Cy,-10.0,10.0,-1,0]),
pnt(pnt[Cx+i*Cy,[536870968],C]),
pnt(pnt[group[Ax+i*Ay,Bx+i*By,Cx+i*Cy,Ax+i*Ay],56]),
pnt(pnt[line[Ax+i*Ay,(Bx+i*By+Cx+i*Cy)/2],56,a]),
pnt(pnt[line[Bx+i*By,(Cx+i*Cy+Ax+i*Ay)/2],56,b]),
pnt(pnt[line[Cx+i*Cy,(Ax+i*Ay+Bx+i*By)/2],56,c]),
pnt(pnt[1/3*Ax+i/3*Ay+1/3*Bx+i/3*By+1/3*Cx+i/3*Cy,56,M]),
pnt(pnt[1/3*Ax+i/3*Ay+1/3*Bx+i/3*By+1/3*Cx+i/3*Cy,56,N]),
0,
"Done",
pnt(pnt[line[-By+Ay+(Ax+Bx)/2+i*(-Ax+Bx+(Ay+By)/2),By-Ay+(Ax+Bx)/2+i*(Ax-Bx+(Ay+By)/2)],56,a]),
pnt(pnt[line[-Cy+By+(Bx+Cx)/2+i*(-Bx+Cx+(By+Cy)/2),Cy-By+(Bx+Cx)/2+i*(Bx-Cx+(By+Cy)/2)],56,b]),
pnt(pnt[line[-Ay+Cy+(Cx+Ax)/2+i*(-Cx+Ax+(Cy+Ay)/2),Ay-Cy+(Cx+Ax)/2+i*(Cx-Ax+(Cy+Ay)/2)],56,c]),
pnt(pnt[((-i)*Ax^2*Bx+Ax^2*By+i*Ax^2*Cx-Ax^2*Cy+i*Ax*Bx^2+i*Ax*By^2+(-i)*Ax*Cx^2+(-i)*Ax*Cy^2+(-i)*Ay^2*Bx+Ay^2*By+i*Ay^2*Cx-Ay^2*Cy-Ay*Bx^2-Ay*By^2+Ay*Cx^2+Ay*Cy^2+(-i)*Bx^2*Cx+Bx^2*Cy+i*Bx*Cx^2+i*Bx*Cy^2+(-i)*By^2*Cx+By^2*Cy-By*Cx^2-By*Cy^2)/(2*Ax*By-2*Ax*Cy-2*Ay*Bx+2*Ay*Cx+2*Bx*Cy-2*By*Cx),56,M]),
pnt(pnt[((-i)*Ax^2*Bx+Ax^2*By+i*Ax^2*Cx-Ax^2*Cy+i*Ax*Bx^2+i*Ax*By^2+(-i)*Ax*Cx^2+(-i)*Ax*Cy^2+(-i)*Ay^2*Bx+Ay^2*By+i*Ay^2*Cx-Ay^2*Cy-Ay*Bx^2-Ay*By^2+Ay*Cx^2+Ay*Cy^2+(-i)*Bx^2*Cx+Bx^2*Cy+i*Bx*Cx^2+i*Bx*Cy^2+(-i)*By^2*Cx+By^2*Cy-By*Cx^2-By*Cy^2)/(2*Ax*By-2*Ax*Cy-2*Ay*Bx+2*Ay*Cx+2*Bx*Cy-2*By*Cx),56,N]),
0,
"Done",
pnt(pnt[0,536870968,A]),
pnt(pnt[1,[536870968],B]),
parameter([x,-10.0,10.0,1.2,0]),
parameter([y,-10.0,10.0,-1,0]),
pnt(pnt[x+i*y,[536870968],C]),
pnt(pnt[group[0,1,x+i*y,0],56]),
pnt(pnt[line[0,1+(x+i*y)/sqrt(x*x+y*y)],56,a]),
pnt(pnt[line[1,x+i*y-sqrt((1-x)*(1-x)+y*y)],56,b]),
pnt(pnt[line[x+i*y,(-x-i*y+1)*sqrt(x*x+y*y)/sqrt((x-1)*(x-1)+y*y)],56,c]),
pnt(pnt[rootof([[2*i*x-2*y-i,(-2*i)*x+i,(-8*i)*x^3+8*x^2*y+12*i*x^2+(-8*i)*x*y^2-8*x*y+(-6*i)*x+8*y^3+4*i*y^2+4*y+i,8*i*x^3+(-12*i)*x^2+8*i*x*y^2+4*x*y+6*i*x+(-4*i)*y^2-2*y-i],[1,0,-4*x^2+4*x-4*y^2-2,0,4*x^2-4*x+1]])/(8*x*y-4*y),56,M]),
pnt(pnt[rootof([[2*i*x-2*y-i,(-2*i)*x+i,(-8*i)*x^3+8*x^2*y+12*i*x^2+(-8*i)*x*y^2-8*x*y+(-6*i)*x+8*y^3+4*i*y^2+4*y+i,8*i*x^3+(-12*i)*x^2+8*i*x*y^2+4*x*y+6*i*x+(-4*i)*y^2-2*y-i],[1,0,-4*x^2+4*x-4*y^2-2,0,4*x^2-4*x+1]])/(8*x*y-4*y),56,N]),
0,
assume[real,1.2],assume[DOM_int,-1],
"Done",
parameter([a,-10.0,10.0,-1,0]),
parameter([b,-10.0,10.0,-1,0]),
pnt(pnt[0,536870968,A]),
pnt(pnt[2,[536870968],B]),
pnt(pnt[a+i*b,[536870968],C]),
pnt(pnt[group[0,2,2*(sqrt(3)*i+1)/2,0],green,T1]),
pnt(pnt[group[2,a+i*b,2+(a+i*b-2)*(sqrt(3)*i+1)/2,2],green,T2]),
pnt(pnt[group[a+i*b,0,a+i*b+(-a-i*b)*(sqrt(3)*i+1)/2,a+i*b],green,T3]),
pnt(pnt[circle(pnt[[0,2*(i*sqrt(3)+3)/3],0,2*pi]),green]),
pnt(pnt[circle(pnt[[2,2*((i*sqrt(3)+3)/6*a+(-sqrt(3)+3*i)/6*b+((-i)*sqrt(3)+3)/3)-2],0,2*pi]),green]),
pnt(pnt[circle(pnt[[a+i*b,-a-i*b+2*(((-i)*sqrt(3)+3)/6*a+(sqrt(3)+3*i)/6*b)],0,2*pi]),green]),
pnt(pnt[group[0,2],56,AB]),
pnt(pnt[group[0,a+i*b],56,AC]),
pnt(pnt[group[a+i*b,2],56,CB]),
pnt(pnt[(i*sqrt(3)+3)/3,536870968,M]),
pnt(pnt[(i*sqrt(3)+3)/6*a+(-sqrt(3)+3*i)/6*b+((-i)*sqrt(3)+3)/3,536870968,N]),
pnt(pnt[((-i)*sqrt(3)+3)/6*a+(sqrt(3)+3*i)/6*b,536870968,O]),
pnt(pnt[group[(i*sqrt(3)+3)/3,(i*sqrt(3)+3)/6*a+(-sqrt(3)+3*i)/6*b+((-i)*sqrt(3)+3)/3],red,MN]),
pnt(pnt[group[(i*sqrt(3)+3)/6*a+(-sqrt(3)+3*i)/6*b+((-i)*sqrt(3)+3)/3,((-i)*sqrt(3)+3)/6*a+(sqrt(3)+3*i)/6*b],red,NO]),
pnt(pnt[group[((-i)*sqrt(3)+3)/6*a+(sqrt(3)+3*i)/6*b,(i*sqrt(3)+3)/3],red,OM]),
0,
0,
assume[DOM_int,-1],assume[DOM_int,-1],
"Done",
parameter([xa,-10.0,10.0,1/2,0]),
parameter([ya,-10.0,10.0,1,0]),
pnt(pnt[xa+i*ya,[536870968],A]),
pnt(pnt[-1,[536870968],B]),
pnt(pnt[group[xa+i*ya,-1],56,AB]),
pnt(pnt[1,[536870968],C]),
pnt(pnt[group[xa+i*ya,1],56,AC]),
pnt(pnt[group[1,-1],56,CB]),
pnt(pnt[xa+i*ya+(-xa-i*ya-1)/2,56,E]),
pnt(pnt[xa+i*ya+(-xa-i*ya+1)/2,56,F]),
pnt(pnt[0,56,G]),
pnt(pnt[circle(pnt[[xa+i*ya+(-xa-i*ya-1)/2,-xa-i*ya-(-xa-i*ya-1)/2+2*((-i)*xa^2+2*xa*ya+i*ya^2+i)/(4*ya)],0,2*pi]),56,c]),
pnt(pnt[line[xa+i*ya,xa+i*ya+2*i],56,H]),
pnt(pnt[xa,56,I]),
"Done",
1,
"Done",
"Done",
"Done",
0,
"Done",
"Done",
pnt(pnt[rootof([[i*xa-ya,(-2*i)*xa,(-4*i)*xa^3+4*xa^2*ya+(-4*i)*xa*ya^2+4*ya^3+4*ya,8*i*xa^3+8*i*xa*ya^2],[1,0,-4*xa^2-4*ya^2-4,0,16*xa^2]])/(8*xa*ya),56,P]),
(2*xa*ya/4/ya-rootof([[-ya,0,4*xa^2*ya+4*ya^3+4*ya,0],[1,0,-4*xa^2-4*ya^2-4,0,16*xa^2]])/8/xa/ya)^2+((-xa^2+ya^2+1)/4/ya-rootof([[xa,-2*xa,-4*xa^3-4*xa*ya^2,8*xa^3+8*xa*ya^2],[1,0,-4*xa^2-4*ya^2-4,0,16*xa^2]])/8/xa/ya)^2,
rootof([[1,-2,-4*xa^2-4*ya^2,8*xa^2+8*ya^2],[1,0,-4*xa^2-4*ya^2-4,0,16*xa^2]])*1/8/ya,
rootof([[xa^2+ya^2-1,-2*xa^2-ya^2+2,-4*xa^4-8*xa^2*ya^2+4*xa^2-4*ya^4,8*xa^4+12*xa^2*ya^2-8*xa^2+4*ya^4],[1,0,-4*xa^2-4*ya^2-4,0,16*xa^2]])/(4*ya^2),
"Done",
0,
"Done",
0,
0,
1,
0,
parameter([xc,-10.0,10.0,-1,0]),
parameter([yc,-10.0,10.0,1,0]),
pnt(pnt[0,[536870968],A]),
pnt(pnt[1,[536870968],B]),
pnt(pnt[xc+i*yc,[536870968],C]),
pnt(pnt[group[0,1,xc+i*yc,0],56]),
pnt(pnt[circle(pnt[[0,2*(i*xc^2+(-i)*xc+i*yc^2+yc)/(2*yc)],0,2*pi]),56,G]),
parameter([x,-10.0,10.0,1,0]),
parameter([y,-10.0,10.0,1,0]),
pnt(pnt[x+i*y,[536870968],M]),
pnt(pnt[1+((x-1)*(xc-1)+y*yc)/((xc-1)*(xc-1)+yc*yc)*(xc+i*yc-1),56,P]),
pnt(pnt[xc+i*yc+(-(x-xc)*xc-(y-yc)*yc)/(xc*xc+yc*yc)*(-xc-i*yc),56,Q]),
pnt(pnt[x,56,R]),
-yc^2*(xc^2*y-xc*y-y^2*yc+y*yc^2-yc*x^2+yc*x),
-x^2*yc+x*yc+xc^2*y-xc*y-y^2*yc+y*yc^2,
(x-yc/(2*yc))^2+(y-(xc^2-xc+yc^2)/(2*yc))^2=((xc^4-2*xc^3+2*xc^2*yc^2+xc^2-2*xc*yc^2+yc^4+yc^2)/(4*yc^2)),
x^2*yc-x*yc-yc^2*y+yc*y^2-y*xc^2+y*xc,
-1,
"Done",
parameter([a1,-10.0,10.0,0.3,0]),
parameter([a2,-10.0,10.0,0.4,0]),
pnt(pnt[0,536870968,A]),
pnt(pnt[1,536870968,B]),
"Done",
"Done",
"Done",
"Done",
"Done",
"Done",
"Done",
"Done",
"Done",
"Done",
"Done",
0,0