[1] | 1 | /* Unfinished !! */
|
---|
| 2 | if get('geometry2, 'version)=false then load(geometry2);
|
---|
| 3 | A:[u1,u2]; B:[u3,u4]; C:[u5,u6]; /* Vertices */
|
---|
| 4 | H1:[x1,x2]; H2: [x3,x4]; H3: [x5,x6]; /* Feet of altitudes */
|
---|
| 5 | M1:[x7,x8]; M2: [x9,x10]; M3: [x11,x12]; /* Midpoints */
|
---|
| 6 | O1:[x13,x14]; /* Centroid */
|
---|
| 7 | O2:[x15,x16]; /* Orthocenter */
|
---|
| 8 | O3:[x17,x18]; /* Intersection point of altitudes */
|
---|
| 9 |
|
---|
| 10 | vars:listofvars([H1,H2,H3,M1,M2,M3,O1,O2,O3,v]);
|
---|
| 11 | params:listofvars([A,B,C]);
|
---|
| 12 |
|
---|
| 13 | eqns:append(
|
---|
| 14 | midpoint(M1,B,C),
|
---|
| 15 | midpoint(M2,A,C),
|
---|
| 16 | midpoint(M3,A,B),
|
---|
| 17 | [
|
---|
| 18 | /* Hi are on the sides */
|
---|
| 19 | collinear(B,C,H1),
|
---|
| 20 | collinear(A,C,H2),
|
---|
| 21 | collinear(A,B,H3),
|
---|
| 22 | /* O3 is intersection of altitudes */
|
---|
| 23 | perpendicular(A,O3,B,C),
|
---|
| 24 | perpendicular(B,O3,A,C),
|
---|
| 25 | perpendicular(C,O3,A,B),
|
---|
| 26 | /* O2 is orthocenter */
|
---|
| 27 | perpendicular(O2,M1,B,C),
|
---|
| 28 | perpendicular(O2,M2,A,C),
|
---|
| 29 | perpendicular(O2,M3,A,B),
|
---|
| 30 | /* O1 is the centroid */
|
---|
| 31 | collinear(O1,M1,A),
|
---|
| 32 | collinear(O1,M2,B),
|
---|
| 33 | /* collinear(O1,M3,C), */
|
---|
| 34 | /* Conclusion */
|
---|
| 35 | 1-v*collinear(O1,O2,O3)
|
---|
| 36 | ]
|
---|
| 37 | );
|
---|
| 38 |
|
---|
| 39 | cover:[[[],[]]];
|
---|
| 40 |
|
---|
| 41 | eqns:expand(eqns);
|
---|
| 42 | cover:expand(cover);
|
---|
| 43 | stringout("euler.out",eqns,vars,params,cover);
|
---|
| 44 |
|
---|