random-simulation-2.cpp

来自「（随机模拟）Stochastic Simulation（已加注）」· C++ 代码 · 共 30 行

// Monte Carlo Simulation
// Written by Microsoft Visual C++
// Copyright by UTLab @ Tsinghua University
// http://orsc.edu.cn/UTLab

#include <stdio.h>
#include <stdlib.h>
#include <math.h>
#include "UTLab.h"

main()
{
  int i, k=0, n=3000;
  double a, b, c;

  for(i=1; i<=n; i++) 
  {
	  a=myexp(1);  // 指数分布
	  b=myn(2,1);  // 正泰分布(2,1)
	  c=myu(0,3);  // [0,3] 均匀分布
	  if((a+b*b+c*c*c)<=30) k++;   
 // 统计循环总次数n内， 满足设定要求要求的次数
  }

  printf(" L =  %3.4f \n",(double)k/n);
  printf(" K =  %d \n",k);
  return 1;
}