; File: fit_xy.pro
; Author: Erik Brisson


N = 10
xy = fltarr(2,N)

openr, 1, 'dat/ex_xy.dat'
readf, 1, xy
close, 1

x = xy(0,*)
y = xy(1,*)
print, 'x'
print, x
print, 'y'
print, y
plot, psym=1, x, y

ypoly = poly_fit(x, y, 3)
yapprox = ypoly(0) + ypoly(1)*x + ypoly(2)*x^2 + ypoly(3)*x^3
print, 'yapprox'
print, yapprox
oplot, x, yapprox

end