**TI89**mainRtestnorm£¥ZI() Prgm ClrIO Toolbar Title "testnorm" Item "Inpt mean,dev",lbl1 Item "Inpt mean",lbl2 Item "Inpt st dev",lbl3 Item "Inpt neither",lbl4 Item "quit",lbl10 EndTBar Lbl lbl1 Disp "mean" Input x Disp "standard deviation" Input s 0y Goto lbl9 Lbl lbl2 Disp "mean" Input x OneVar c1 ¨(sum(seq((dist[1][i]-x)^2,i,1,nStat))/nStat)s:1y Goto lbl9 Lbl lbl3 Disp "standard deviation" Input s OneVar c1 šx:1y Goto lbl9 Lbl lbl4 OneVar c1 šx:¨((nStat-1)/nStat)*Sxs:2y Goto lbl9 Lbl Lbl10 Stop Lbl lbl9 Disp "number of bins" Input m m-1-yl dist[1]list1 SortA list1 NewData dist,list1 OneVar c1 Define y2(x)=­3-6/(m-2)+6x/(m-2) 0j 0k Lbl lbl5 If ((dist[1])[j+k+1]-x)/sœy2(1) Then k+1k Goto lbl5 EndIf {k}list2 kj 2i Lbl lbl6 If m>i Then 0k Lbl lbl7 If j+k+1œnStat Then If ((dist[1])[j+k+1]-x)/s>y2(i-1) and ((dist[1])[j+k+1]-x)/sœy2(i) Then k+1k Goto lbl7 Else augment(list2,{k})list2 j+kj i+1i Goto lbl6 EndIf Else augment(list2,{k})list2 j+kj i+1i Goto lbl6 EndIf Else augment(list2,{nStat-j})list2 EndIf {nStat*(.5-s\nrmcdff(0,3,0,1))}list3 1j While jœm-2 augment(list3,{nStat*s\nrmcdff(y2(j),y2(j+1),0,1)})list3 j+1j EndWhile augment(list3,{list3[1]})list3 sum(seq((list2[i]-list3[i])^2/list3[i],i,1,m))z 1-s\Chi2CdfF(0,z,l)b NewData dist,list1,list2,list3 ClrIO Disp "test statistic" Disp z Disp "p value" Disp b EndPrgm äåܵ§