mc1.java

MC算法。简单的Monte Carlo 算法的java源码实现。

import java.io.*;

public class MC1
{
	static double r=0;
	public static void main(String[] args) 
	{
		int i;
		int k,m=8,n=100;
		double []x=new double[8];
		double []a=new double[800];
		double av, s=0, dx=0, si, t;

		for(i=0; i<n; i++)
		{
			MC1.rand(m,x);
			for(k=0; k<m;k++)
			{
				a[i*8+k]=x[k];
				s+=x[k];
				//System.out.print("x("+k+")	"+x[k]);
			}
			//System.out.println();
		}
		av=s/800;
		for(i=0;i<8;i++)
		{
			t=x[i]-av;
			dx=dx+t*t;
		}
		dx=dx/8;
		si=Math.sqrt(dx);
		System.out.println("av="+av+"	"+"dx="+dx+"	"+"si="+si);
	}

	public static void rand(int m, double x[])
	{
		int i;
		double[] q=new double[8];
	
		q[0]=0.4142135623;	q[1]=0.7320508075;
		q[2]=0.2360679774;	q[3]=0.6457513110;
		q[4]=0.3166247904;	q[5]=0.6055512755;
		q[6]=0.1231056256;	q[7]=0.3588989435;

		for(r++, i=0;i<m;i++)
			x[i]=MC.fmod(r*q[i],1);
	}
}