/*===============================================================================================
probable.c
公算誤差確率計算
計算精度：小数点以下５桁まで保証
１９９９．９．１２
ｂｙ Ｔ．Ｋｕｂｏｔａ
================================================================================================*/
#include <stdio.h>
#include <math.h>
#define PI 3.141592653589793238462
#define EPS (1e-6)
double prob_inte(double); /* 確率積分関 */
double fact(long); /* 階乗関数 */
int main(void)
{
double a0,b0; /* 積分区間 */
double a,b; /* 正規化積分区間 */
double pa,pb,pab; /* 確率 */
double r; /* 公算誤差 */
printf("********** 公算誤差確率計算 **********\n\n");
printf(" 公算誤差を入力してください。--->");
scanf("%lf",&r);
printf("\n誤差区間 ａ〜ｂ（０＜ａ＜ｂ）を入力して下さい．--->);
scanf("%lf %lf",&a0,&b0);
a=a0*(0.476936/r);
b=b0*(0.476936/r);
if(a0==0){
pa=0;
}
else{
pa=prob_inte(a);
}
pb=prob_inte(b);
pab=pb-pa;
printf("誤差区間 %g 〜 %g の確率は %7.5lf です．\n",a0,b0,pab);
return 0;
}
/*
*/
/*