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

#define PI 3.141592654
using namespace std;

void gaussian_vals(double x[], double vals[], int length, double mean, double sig);

double gaussian(double x, double mean , double sig);

void main()
{
	double values,x[100],mean,sigma,vals[100]={0};



	cout<<"enter mean of gaussian"<<endl;
	cin>>mean;

	cout<<"enter standard deviation of gaussian"<<endl;
	cin>>sigma;

	cout<<"how many values to be entered? (max 100)"<<endl;
	cin>>values;

	cout<<"enter values"<<endl;

		for(int i=0;i<values;i++)
		{
			cin>>x[i];
		}

		gaussian_vals(x,vals,values,mean,sigma);

		for(int j=0;j<values;j++)
			cout<<"gauss("<<x[j]<<")"<<setw(21)<<"="<<vals[j]<<endl;

}

double gaussian(double x , double mean , double sig)
{


	double g=0;
	double r =((x-mean)/sig)*((x-mean)/sig);
	g = 1/(sig*sqrt(2*PI))*exp(-0.5*r);

	return g;
}

void gaussian_vals(double x[], double vals[], int length, double mean, double sig)
{

    for (int i=0;i<length;i++)

        vals[i]=gaussian(x[i],mean,sig);
}