RosettaCodeData/Task/Polynomial-regression/Emacs-Lisp/polynomial-regression.el

(let ((x '(0 1 2 3 4 5 6 7 8 9 10))
      (y '(1 6 17 34 57 86 121 162 209 262 321)))
  (calc-eval "fit(a*x^2+b*x+c,[x],[a,b,c],[$1 $2])" nil (cons 'vec x) (cons 'vec y)))