fourier.c running 
fourier series of sin2.dat for all terms 
npts=0, x=0.000000, y=0.000000 
npts=1, x=1.000000, y=0.382680 
npts=2, x=2.000000, y=0.707100 
npts=3, x=3.000000, y=0.923880 
npts=4, x=4.000000, y=1.000000 
npts=5, x=5.000000, y=0.923880 
npts=6, x=6.000000, y=0.707100 
npts=7, x=7.000000, y=0.382680 
npts=8, x=8.000000, y=0.000000 
npts=9, x=9.000000, y=-0.382680 
npts=10, x=10.000000, y=-0.707100 
npts=11, x=11.000000, y=-0.923880 
npts=12, x=12.000000, y=-1.000000 
npts=13, x=13.000000, y=-0.923880 
npts=14, x=14.000000, y=-0.707100 
npts=15, x=15.000000, y=-0.382680 
npts=16, x=16.000000, y=0.000000 
npts=17, x=17.000000, y=0.382680 
npts=18, x=18.000000, y=0.707100 
npts=19, x=19.000000, y=0.923880 
npts=20, x=20.000000, y=1.000000 
npts=21, x=21.000000, y=0.923880 
npts=22, x=22.000000, y=0.707100 
npts=23, x=23.000000, y=0.382680 
npts=24, x=24.000000, y=0.000000 
npts=25, x=25.000000, y=-0.382680 
npts=26, x=26.000000, y=-0.707100 
npts=27, x=27.000000, y=-0.923880 
npts=28, x=28.000000, y=-1.000000 
npts=29, x=29.000000, y=-0.923880 
npts=30, x=30.000000, y=-0.707100 
npts=31, x=31.000000, y=-0.382680 
npts=32, x=32.000000, y=0.000000 
33 points read, using 15 terms 
xmin=0.000000, xmax=32.000000, ymin=-1.000000, ymax=1.000000 
coefficients 
a[0]=-0.000000, b[0]=0.000000 
a[1]=-0.000000, b[1]=0.000000 
a[2]=0.000000, b[2]=0.999997 
a[3]=-0.000000, b[3]=-0.000000 
a[4]=0.000000, b[4]=-0.000000 
a[5]=-0.000000, b[5]=0.000000 
a[6]=-0.000000, b[6]=-0.000004 
a[7]=-0.000000, b[7]=0.000000 
a[8]=0.000000, b[8]=0.000000 
a[9]=-0.000000, b[9]=-0.000000 
a[10]=0.000000, b[10]=0.000001 
a[11]=-0.000000, b[11]=-0.000000 
a[12]=0.000000, b[12]=0.000000 
a[13]=0.000000, b[13]=-0.000000 
a[14]=-0.000000, b[14]=0.000002 

y[0]=0.000000, approx=-0.000000, err=1.22884e-16 
y[1]=0.382680, approx=0.382680, err=-4.996e-16 
y[2]=0.707100, approx=0.707100, err=2.22045e-16 
y[3]=0.923880, approx=0.923880, err=-5.55112e-16 
y[4]=1.000000, approx=1.000000, err=8.88178e-16 
y[5]=0.923880, approx=0.923880, err=-3.33067e-16 
y[6]=0.707100, approx=0.707100, err=-7.77156e-16 
y[7]=0.382680, approx=0.382680, err=1.27676e-15 
y[8]=0.000000, approx=0.000000, err=-1.12888e-15 
y[9]=-0.382680, approx=-0.382680, err=3.33067e-16 
y[10]=-0.707100, approx=-0.707100, err=2.22045e-16 
y[11]=-0.923880, approx=-0.923880, err=-4.44089e-16 
y[12]=-1.000000, approx=-1.000000, err=0 
y[13]=-0.923880, approx=-0.923880, err=-1.11022e-15 
y[14]=-0.707100, approx=-0.707100, err=6.66134e-16 
y[15]=-0.382680, approx=-0.382680, err=6.10623e-16 
y[16]=0.000000, approx=0.000000, err=-1.28365e-15 
y[17]=0.382680, approx=0.382680, err=-1.27676e-15 
y[18]=0.707100, approx=0.707100, err=1.88738e-15 
y[19]=0.923880, approx=0.923880, err=-7.77156e-16 
y[20]=1.000000, approx=1.000000, err=4.44089e-16 
y[21]=0.923880, approx=0.923880, err=1.33227e-15 
y[22]=0.707100, approx=0.707100, err=-1.55431e-15 
y[23]=0.382680, approx=0.382680, err=8.32667e-16 
y[24]=0.000000, approx=0.000000, err=-1.27665e-15 
y[25]=-0.382680, approx=-0.382680, err=8.32667e-16 
y[26]=-0.707100, approx=-0.707100, err=1.33227e-15 
y[27]=-0.923880, approx=-0.923880, err=-9.99201e-16 
y[28]=-1.000000, approx=-1.000000, err=6.66134e-16 
y[29]=-0.923880, approx=-0.923880, err=-7.77156e-16 
y[30]=-0.707100, approx=-0.707100, err=5.55112e-16 
y[31]=-0.382680, approx=-0.382680, err=1.55431e-15 
y[32]=0.000000, approx=-0.000000, err=6.12728e-16 
avgerr=8.23769e-16, rmserr=1.06662e-16, maxerr=1.88738e-15 
fourier.c finished