#include <stdio.h>
#include <math.h>

#define PI 3.1415926535897932384626433 

//プロトタイプ宣言
double PNormal(double a, double b, double mean, double var);

void main(void)
{
	double a=8, b=14, mean=10, var=5;
	printf("平均 %lf 標準偏差 %lf\n",mean, var);
	printf("%lf から %lf までの正規分布に従う確率は",a,b);
	//引数は左から、開始範囲、終了範囲、平均、標準偏差
	printf(" %lf です\n",PNormal(a, b, mean, var));
}

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;
}