10 VAR A,B,LOGPROB,LOG_NCR,N,R,T,X,M,Z,P,Y,Q,U,S =REAL.DOUBLE
15 VAR C$ = STRING
20 INPUT "N = ";N
30 LPRINTER
40 PRINT "N = ";N
45 CONSOLE
50 INPUT "R = ";R
60 LPRINTER
70 PRINT "R = ";R
80 CONSOLE
90 PRINT "RUNTIME = "; 0.4*R;"SECONDS"
100 INPUT "THETA =";T
110 LPRINTER
120 PRINT "THETA =";T
130 CONSOLE
140 PRINT "CALCULATIONS IN PROGRESS"
150 FOR X = 1 TO R
160 Y = Y + LOG(X)/LOG(10)
170 NEXT X
180 M = N - R + 1
190 FOR Z = M TO N
200 P = P + LOG(Z)/LOG(10)   
210 NEXT Z
220 LOG_NCR = P - Y
230 S = 1 - T
240 Q = R*LOG(T)/LOG(10)
250 U = (N - R)*LOG(S)/LOG(10)
260 LOGPROB =  LOG_NCR + Q + U
270 PRINT
280 LPRINTER
290 PRINT
295 PRINT "LOGPROB =";LOGPROB
300 CONSOLE
310 PRINT
315 PRINT "LOGPROB=";LOGPROB
320 A = 10^(LOGPROB - INT(LOGPROB))
330 B = INT(LOGPROB)
340 C$= "PROBABILITY = #.## TEN +#####"
350 PRINT USING C$,A,B
360 LPRINTER
370 PRINT USING C$,A,B
380 PRINT "--ooO0Ooo--"
390 PRINT
400 PRINT
410 CONSOLE
420 PRINT "--ooO0Ooo--"
430 PRINT
440 PRINT
450 END