c语言泊松分布的计算

【问题描述】 泊松分布是一种常用的离散型概率分布，数学期望为m的泊松分布的分布函数定义如下： P(m, k) = mk * e-m/k! (k = 0, 1, 2, 3, …) 对于给定的m和k (0<m<2000, 0<= k < 2500)，计算其概率，以科学格式输出，保留小数点后6位有效数字。 可以使用数学库函数，误差不超过0.000001。 【输入形式】 输入文件为当前目录下的poisson.in。文件中包含两个数字，分别为m,k的值。 【输出形式】 输出文件为当前目录下的poisson.out。文件中输出泊松分布的值，值以科学格式输

#include
#include
#include
#define MAX_INTER 1E280
#define MIN_INTER 1E-280
#define M_E 2.7182818284590452354
double poisson（double mint k）
{
	double x=1.0;
	int emmknk;
	em=mk=nk=0;
	while（em		while（mk			x*=m;
			mk++;
		}
		if（x			return 0.0;
		while（nkMIN_INTER）{
			nk++;
			x/=nk;
		}
		while（emMI