**TI86** L < DISCRETE< : Menu(1,"BINOMIAL",A,2,"GEOMETRIC",B,3,"HYPERGEOM",C,4,"NEGBINOM",D,5,"POISSON",E,6,"QUIT",F)ÖLbl AÖClLCDÖDisp "NUMBER OF TRIALS"ÖInput NÖDisp "PROBABILITY"ÖInput PÖDisp "LOWER BOUND"ÖInput JÖDisp "UPPER BOUND"ÖInput KÖIf J==0ÖThenÖbicdf(N,P,K)CÖElseÖbicdf(N,P,K)-bicdf(N,P,J-1)CÖEndÖint ((N+1)P)MÖbipdf(N,P,M)AÖIf M1ÖThenÖbipdf(N,P,M-1)BÖElseÖ0BÖEndÖDisp "COMPLETE DIST?"ÖInput ZÖIf Z==1ÖThenÖClrLstsÖFnOff Öseq(I,I,0,N)xStatÖbipdf(N,P)yStatÖbicdf(N,P)fStatÖ0xMinÖNxMaxÖ1xSclÖ0yMinÖAyMaxÖ1ySclÖPlOffÖPlOn 1ÖPlot1(1,xStat,yStat,1)ÖEndÖClLCDÖOutpt(2,1,"PROB")ÖOutpt(2,8,C)ÖOutpt(4,1,"MEAN")ÖOutpt(4,9,N*P)ÖOutpt(5,1,"ST.DEV")ÖOutpt(5,9,(N*P(1-P)))ÖIf round(A)==round(B)ÖThenÖOutpt(7,1,"BIMODAL")ÖOutpt(7,9,{M-1,M})ÖElseÖOutpt(7,1,"MODE")ÖOutpt(7,9,M)ÖEndÖStopÖLbl BÖClLCDÖDisp "PROBABILITY"ÖInput PÖDisp "LOWER BOUND"ÖInput JÖDisp "UPPER BOUND"ÖInput KÖ(1-P)^(J-1)-(1-P)^KCÖDisp "COMPLETE DIST?"ÖInput ZÖIf Z==1ÖThenÖClrLstsÖint (ln(.025)/ln(1-P))+1MÖseq(I,I,1,M)xStatÖseq(gepdf(P,I),I,1,M)yStatÖseq(gecdf(P,I),I,1,M)fStatÖEndÖClLCDÖOutpt(2,1,"PROB")ÖOutpt(2,8,C)ÖOutpt(4,1,"MEAN")ÖOutpt(4,8,1/P)ÖOutpt(6,1,"ST.DEV")ÖOutpt(6,8,(1-P)/P)ÖStopÖLbl CÖClLCDÖDisp "POP. SIZE"ÖInput NÖDisp "TYPE A SIZE"ÖInput MÖDisp "SAMPLE SIZE"ÖInput RÖDisp "LOWER BOUND"ÖInput JÖDisp "UPPER BOUND"ÖInput KÖsum seq(M nCr I*(N-M) nCr (R-I)/N nCr R,I,J,K,1)CÖDisp "COMPLETE DIST?"ÖInput ZÖIf Z==1ÖThenÖClrLstsÖmax(0,R-(N-M))LÖseq(I,I,L,min({M,R,998}),1)xStatÖseq(M nCr I*(N-M) nCr (R-I)/N nCr R,I,L,min({M,R,998}),1)yStatÖFor(I,1,dimL xStat,1)Ösum seq(yStat(Z),Z,1,I,1)fStat(I)ÖEndÖEndÖClLCDÖOutpt(2,1,"PROB")ÖOutpt(2,8,C)ÖOutpt(4,1,"MEAN")ÖOutpt(4,8,R*M/N)ÖOutpt(6,1,"ST.DEV")ÖOutpt(6,8,(R*M/N*(N-M)/N*(N-R)/(N-1)))ÖStopÖLbl DÖClLCDÖDisp "PROBABILITY"ÖInput PÖDisp "NO. OF OCCURRENCES"ÖInput RÖDisp "LOWER BOUND"ÖInput JÖDisp "UPPER BOUND"ÖInput KÖsum seq((M-1) nCr (R-1)*P^R*(1-P)^(M-R),M,J,K,1)CÖDisp "COMPLETE DIST?"ÖInput ZÖIf Z==1ÖThenÖClrLstsÖ2*int (R/P)LÖseq(I,I,R,L)xStatÖseq((M-1) nCr (R-1)*P^R*(1-P)^(M-R),M,R,L)yStatÖcSum(yStat)fStatÖEndÖClLCDÖOutpt(2,1,"PROB")ÖOutpt(2,8,C)ÖOutpt(4,1,"MEAN")ÖOutpt(4,8,R/P)ÖOutpt(6,1,"ST.DEV")ÖOutpt(6,8,(R*(1-P))/P)ÖStopÖLbl EÖClLCDÖDisp "AVERAGE"ÖInput PÖDisp "LOWER BOUND"ÖInput JÖDisp "UPPER BOUND"ÖInput KÖIf J==0ÖThenÖpscdf(P,K)CÖElseÖpscdf(P,K)-pscdf(P,J-1)CÖEndÖint (P)RÖpspdf(P,R)AÖIf R1ÖThenÖpspdf(P,R-1)BÖElseÖ0BÖEndÖDisp "COMPLETE DIST?"ÖInput ZÖIf Z==1ÖThenÖClrLstsÖFnOff ÖR+1MÖWhile pscdf(P,M)<.975ÖM+1MÖEndÖseq(I,I,0,M)xStatÖseq(pspdf(P,I),I,0,M)yStatÖseq(pscdf(P,I),I,0,M)fStatÖ0xMinÖMxMaxÖ1xSclÖ0yMinÖAyMaxÖ1ySclÖPlOffÖPlOn 1ÖPlot1(1,xStat,yStat,1)ÖEndÖClLCDÖOutpt(2,1,"PROB")ÖOutpt(2,8,C)ÖOutpt(4,1,"MEAN")ÖOutpt(4,9,P)ÖOutpt(5,1,"ST.DEV")ÖOutpt(5,9,(P))ÖIf round(A)==round(B)ÖThenÖOutpt(7,1,"BIMODAL")ÖOutpt(7,9,{R-1,R})ÖElseÖOutpt(7,1,"MODE")ÖOutpt(7,9,R)ÖEndÖStopÖLbl FÖStopÖÖÖç½