(in-package :division)
(use-package :symbolic-polynomial)

(let ((f (string->poly "x^2" '(x y)))
      ;;(y-sq (string->poly "y^2" '(x y)))
      (fl (cdr (string->poly "[x+y,x-2*y]" '(x y))))
      ;;(quotients (cdr (string->poly "[x-y,0]" '(x y))))
      )
  (multiple-value-bind (f c division-count)
      (normal-form f fl)
    (format t "Quotient:        ~S~%" (->sexp f))
    (format t "Coefficient:     ~S~%" c)
    (format t "Division count:  ~S~%" division-count)))