# The (equation of the ) Conic(s) through the inputed Points
ConicThroughPoints:=proc() local Conic,A,B,C,D,E,F,i:
Conic:=A*x^2+B*x*y+C*y^2+D*x+E*y+F:
end:
#It plugs in the inputed points, args[1],args[2], .., args[nargs],
# into a general Conic, solves the resulting
# system of equations and plugs back into the Conic
subs(solve({seq(subs({x=args[i][1],y=args[i][2]},Conic),
i=1..nargs)},{A,B,C,D,E,F}),Conic):
end: