/* A0320.java　３章２．解答
 *        幾何分布シミュレーション  
 *         (C) H.Ishikawa 2008 
 */

package simulation;
import java.applet.*;
import java.awt.*;
import java.text.*;


public class A0320 extends Applet{

	public void paint(Graphics g) {
		
		int M = 10000;					/*発生させる乱数の数*/
		int    j;                      /* forのカウンタ */
		int    k = 0;                  /* 失敗続くの回数 */
		int    n = 0;                  /* 成功の回数の合計 */
		int f[] = new int [100];      /* 度数 */
		double p = 0.5;                /* 成功の確率 */
		double g0;                     /* 幾何分布理論値 */
		DecimalFormat exFormat1 = new DecimalFormat(" #.000"); 
	
		/* メイン */
		for (j = 1; j <= M; j ++) {
			if (p <Math.random()) { 	/* 失敗 */
				k = k + 1;
			} else {
				n = n + 1;
				f[k] = f[k] + 1;
				k = 0;
			}
		}
		g.drawString("      k        simulation          theory", 10, 10);
		for (k = 0; k < 16; k ++) {
			g0 = p * Math.pow((1 - p) , (double)(k)); 	/* 幾何分布理論値 */
			g.drawString("     " + k, 10 , 20 + 10 * k);
			g.drawString("     " + exFormat1.format((double)f[k] / n), 40 , 20 + 10 * k);
			g.drawString("     " + exFormat1.format( g0), 120 , 20 + 10 * k);
		}
	}
}