random-simulation-1.cpp

来自「很好的随机模拟例题」· C++ 代码 · 共 27 行

// 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, n=3000;
  double a, b, c, sum=0;

  for(i=1; i<=n; i++) {
	  a=myexp(1);
	  b=myn(2,1);
	  c=myu(0,3);
	  sum=sum+a+b*b+c*c*c;
  }

  printf("L = %3.4f\n",sum/n);
  return 1;
}