用C#实现计算标准正态分布

最新推荐文章于 2024-09-30 08:30:00 发布

ling000000

最新推荐文章于 2024-09-30 08:30:00 发布

阅读量6.7k

点赞数

文章标签： c# button object

本文链接：https://blog.csdn.net/ling000000/article/details/1621342

版权

static double f(double x)
{
   double f=Math.Exp(-x*x/2)/Math.Sqrt(2*Math.PI);
   return f;
}
/// <summary>
/// 应用程序的主入口点。
/// </summary>
[STAThread]
static void Main()
{
   Application.Run(new Form1());
}
private void button1_Click(object sender, System.EventArgs e)
{
   int S=2;
   double Q=0;
   double b=Convert.ToDouble(textBoxU.Text);
   while(true)
   {

    double a=b-S;
    int M=1,N=1,k=1,m=1;
    double ep,I,h;
    ep=0.000000000001;
    h=b-a;
    I=h*(f(a)+f(b))/2;
    double[,] T=new double[5000,5000];
    T[1,1]=I;
    while(1>0)
    {
     N=(int)Math.Pow(2,m-1);
     if(N>5000)
     {
      MessageBox.Show("请缩小精度再计算!");
      break;
     }
     else
     {
      h=h/2;
      I=I/2;
      for(int i=1;i<=N;i++)
       I=I+h*f(a+(2*i-1)*h);
      T[m+1,1]=I;
      M=2*N;
      k=1;
      while(M>1)
      {
       T[m+1,k+1]=(Math.Pow(4,k)*T[m+1,k]-T[m,k])/(Math.Pow(4,k)-1);
       M=M/2;
       k=k+1;
      }
      if(Math.Abs(T[k,k]-T[k-1,k-1])<ep)
       break;
      m=m+1;
     }
    }
    I=T[k,k];
    Q=Q+I;
    if(Math.Abs(I)<ep)
     break;
    b=a;S=2*S;
    textBoxResult.Text=Convert.ToString(Math.Round(Q,10));
   }