#include "iostream.h"
#include "math.h"
#include "conio.h"

int round(double x)
{
	int temp;
	double temp1;
	temp=int(x);
	temp1=x-temp;
	if(temp1>0.5)
	{
		return temp=temp+1;
	}
	else
	{
		return temp;
	}
}

void main()
{
    //program d13r1
    //driver for routine moment
    const double pi = 3.14159265;
    int npts = 10000;
    int nbin = 100;
    int ndat = npts + nbin;
    double data[100101];
    int nlim,j,k,i = 1;
	double x;
    for (j = 1; j<=nbin; j++)
	{
		double temp=0;
        x = pi * j / nbin;
        nlim = round(sin(x)*pi/2.0*(double)npts/(double)nbin);
        for (k = 1; k<=nlim; k++)
		{
            data[i] = x;
            i = i + 1;
        }
    }
	data[npts]=0.0;
	cout.setf(ios::fixed|ios::left);
	cout.precision(4);
    cout<<endl;
    cout<<"Moments of a sinusoidal distribution"<<endl;
	double adev,sdev,var,ave,skew,curt;
    moment(data, npts, ave, adev, sdev, var, skew, curt);
    cout<<"               "<<"Calculated     Expected"<<endl;
    cout<<"Mean :                 ";
	cout.width(14);	   cout<<ave;
	cout.width(14);    cout<<pi / 2.0<<endl;
    cout<<"Average Deviation :    ";
	cout.width(14);    cout<<adev;
	cout.width(14);    cout<<0.5708<<endl;
    cout<<"Standard deviation :   ";
	cout.width(14);    cout<<sdev;
	cout.width(14);    cout<<0.6837<<endl;
    cout<<"Variance :             ";
	cout.width(14);    cout<<var;
	cout.width(14);    cout<<0.4674<<endl;
    cout<<"Skewness :             ";
	cout.width(14);    cout<<skew;
	cout.width(14);    cout<<0<<endl;
    cout<<"Kurtosis :             ";
	cout.width(14);    cout<<curt;
	cout.width(14);    cout<<-0.8062;
	cout<<endl;
}