1000	!
	!	WRITTEN BY JEFFREY ROBBINS
	!	CLASS OF '76 NORTHFIELD MT. HERMON
	!
	!	KEYED IN BY TIMOTHY SHANE MCCOY 
	!	AT CENTRAL STATE UNIVERSITY
	!	ON PDP 11/45
	!	VERSION RSTS/E  V06C
	!
1100	PRINT "POPULATION GENETICS PROGRAM"
1200	PRINT
	\ INPUT "INSTRUCTIONS";Y$
	\ IF Y$ = "NO"
	  THEN 1400
	  ELSE PRINT "     *GENES*  IS A SCIENCE PROGRAM"
	  \ PRINT "    IT WILL PRINT THE PER CENTAGE AND AMOUNT OF PERSON CARRYING A"
	  \ PRINT "DOMINATE TRAIT, AND A RECESSIVE TRAIT."
	  \ PRINT
	  \ PRINT
1300	PRINT "*NOTE*  TYPE '0' (ZERO) WHEN THE COMPUTER PRINTS 'TOTAL POPULATION?'"
	  \PRINT "AND ANOTHER WHEN THE COMPUTER ASKS FOR THE TOTAL HOMOZYGOUS RECESSIVE"
	\ PRINT "TO END THE PROGRAM."
	\ PRINT
	\ PRINT
1400	PRINT
	\ INPUT "TOTAL POPULATION";T;"TOTAL HOMOZYGOUS RECESSIVE";R
	\ GOTO 1400 IF T < R
1500	GOTO 2400 IF T = 0
1600	P = T - R
	\ D = P/T
	\ D1 = R/T
1700	M(2),M(4) = SQR(D1) 
	\ M(1),M(3) = 1.0 - M(2)
1800	PRINT
	\ PRINT
	\ PRINT "GENOTYPE","PER CENT","TOTAL"
	\ PRINT "========","========","====="
1900	PRINT "TTL. DOMINANT",D * 100;"%",T * D
2000	PRINT "HOM. DOMINANT",(M(1) * M(3)) * 100;"%",T * (M(1) * M(3))
2100	PRINT "HETEROZYGOUS",(2 * M(2) * M(3)) * 100;"%",T * (2 * M(2) * M(3))
2200	PRINT "HOM. RECESSIVE";(M(2) * M(4)) * 100;"%",T * (M(2) * M(4))
2300	PRINT
	\ PRINT
	\ GOTO 1400
2400	END