#define PI 3.1415926535897932384626433 

double PNormal(double a, double b, double mean, double var)
{
	double acc=1000;	//精度。値を大きくすれば精度は良くなりますが、計算時間はかかります
	double c,i, p=0 ,var2=var*var;
	
	//b>=aの形にする
	if(a > b)
	{
		i = a; a = b; b = i;
	}

	c = (b-a)/acc;
	
	for(i=a; i<b; i+=c)
		p += exp(-0.5 / var2 * pow((i-mean),2)) / (sqrt(2 * PI) * var) * c;
	return p;
}