
/* input for the "integrated calculator" that came across the 
   net a few years ago */

define movavg(cur,avg,factor) {
	return (cur*8+factor*avg)/(factor+8);
}
define sim(once,factor,reps) {
	auto avg;
	avg = movavg(once,0,factor);
	while (reps > 0) {
		reps--;
		avg = movavg(0,avg,factor);
	}
	return avg;
}
define run(len,edu,hap) {
	auto i;
	printf ("reps	education	happiness\n");
	newi = 0;
	while (newi < len) {
		printf("%d	%f	%f\n", newi, sim(edu*5000, 48*4, newi), 
			sim(hap*6000, 48, newi));
		newi++;
	}
	return 0;
}
