DECLARE SUB PolyCurveFit (indvar!(), depvar!(), numobs%, order%, coef!(), yest!(), resid!(), see!, coefsig!(), rsq!, r!, cferror%)
n% = 5
DIM indvar!(0 TO n% - 1)
DIM depvar!(0 TO n% - 1)
order% = 2
DIM coeff!(0 TO order%)
DIM yest!(0 TO n% - 1)
DIM resid!(0 TO n% - 1)
DIM coefferr!(0 TO order%)

indvar!(0) = 1
indvar!(1) = 2
indvar!(2) = 3
indvar!(3) = 4
indvar!(4) = 5

depvar!(0) = 2 + 1
depvar!(1) = 2 + 4
depvar!(2) = 2 + 9
depvar!(3) = 2 + 16
depvar!(4) = 2 + 25

CALL PolyCurveFit(indvar!(), depvar!(), n%, order%, coeff!(), yest!(), resid!(), see!, coefferr!(), rsqval!, r!, cferror%)

PRINT
FOR i% = 0 TO order%
 PRINT "c("; i%; ") = "; coeff!(i%)
NEXT i%

END