# include<iostream.h>
# include<math.h>
# include<iomanip.h>

void main()
{
    //program d2r8
    //driver for routine polcof
    int nfunc,i,j,np = 5;
    double x,f,sum,xa[6], ya[6], coeff[6],pi = 3.141593;
    for (nfunc = 1; nfunc<=2; nfunc++)
	{
        if (nfunc == 1)
		{
            cout<<endl;
            cout<<"sine function from 0 to pi"<<endl;
            for (i = 1; i<=np; i++)
			{
                xa[i] = i * pi /double(np);
                ya[i] = sin(xa[i]);
            }
		}
        else 
		{
			if (nfunc ==2)
			{
            cout<<endl;
            cout<< "Exponential function from 0 to 1"<<endl;
            for (i = 1;i<=np;i++)
				{
					xa[i] = 1.0 * i / double(np);
					ya[i] = exp(xa[i]);
				}
			}
            else
            return;
        }
        polcof(xa, ya, np, coeff);
        cout<<"coefficients"<<endl;
        for (i = 1;i<=np;i++)
            cout<<setw(12)<<coeff[i];
		cout<<endl;
        cout<< "           x          f[x]       polynomial"<<endl;
        for (i = 1; i<=10; i++)
		{
            if (nfunc == 1 )
			{
                x = (-0.05 + i / 10.0) * pi;
                f = sin(x);
			}
            else 
			{
				if (nfunc == 2)
				{
					x = (-0.05 + double(i) / 10.0);
					f = exp(x);
				}
            }
            sum = coeff[np];
            for (j = np - 1; j>=1; j--)
				{
					sum = coeff[j] + sum * x;
				}
            cout<< setw(14)<<x;
            cout<< setw(14)<<f;
            cout<< setw(14)<<sum<<endl;
		}
        cout<<  "**********************************"<<endl;
    }
}