source: CGBLisp/trunk/examples/euler.macsyma

/* Unfinished !! */
if get('geometry2, 'version)=false then load(geometry2);
A:[u1,u2]; B:[u3,u4]; C:[u5,u6]; /* Vertices */
H1:[x1,x2]; H2: [x3,x4]; H3: [x5,x6]; /* Feet of altitudes */
M1:[x7,x8]; M2: [x9,x10]; M3: [x11,x12]; /* Midpoints */
O1:[x13,x14]; /* Centroid */
O2:[x15,x16]; /* Orthocenter */
O3:[x17,x18]; /* Intersection point of altitudes */

vars:listofvars([H1,H2,H3,M1,M2,M3,O1,O2,O3,v]);
params:listofvars([A,B,C]);

eqns:append(
        midpoint(M1,B,C),
        midpoint(M2,A,C),
        midpoint(M3,A,B),
        [
                /* Hi are on the sides */
                collinear(B,C,H1),
                collinear(A,C,H2),
                collinear(A,B,H3),
                /* O3 is intersection of altitudes */
                perpendicular(A,O3,B,C),
                perpendicular(B,O3,A,C),
                perpendicular(C,O3,A,B),
                /* O2 is orthocenter */
                perpendicular(O2,M1,B,C),
                perpendicular(O2,M2,A,C),
                perpendicular(O2,M3,A,B),
                /* O1 is the centroid */
                collinear(O1,M1,A),
                collinear(O1,M2,B),
                /* collinear(O1,M3,C),  */
                /* Conclusion */
                1-v*collinear(O1,O2,O3)
        ]
);

cover:[[[],[]]];

eqns:expand(eqns);
cover:expand(cover);
stringout("euler.out",eqns,vars,params,cover);