Re: WHO STOPS HIM?

THREE DIFFENT ATTEMPTS Steve Lawford, T. Thadewald & H. Buning, and me, using Monte Carlo simulations achieved practically the same RESULTS.
It is not an uneducated persons as John Smith, Jack Tomsky) that are sufficiently credible to STATE THE VALUES´s WORTHLESS:
I INVITED REPETIBLY ALL THE UNBELIEVERS TO TEST MY PROGRAM (JB) IN ORDER TO FIND
********************WHERE****************
********************WHY******************
IT IS WRONG:
I am absolutely AVAILABLE to discuss IN DETAIL the program (listing below). EVERYONE, THE READERS; are free to
___RUN IT in whatever computer
___TRANSLATE IT IN OTHER PROGRAMING LANGUAGE
___ASK AN EXPERT TO ANALYSE IT.

OK?.
____


Luis Amaral Afonso


REM "JB"
CLS
DEFDBL A-Z
INPUT " N = "; n
LOCATE 1, 15
PRINT " CHI2 , 2DF ";
PRINT " 4.605 5.991 9.210"
PRINT " 10% 5% 1% "
DIM x(n)
DIM JB(8004)
pi = 4 * ATN(1): ju = 1
DEF fnJB (S, K, n) = (n / 6) * (S ^ 2 + (K - 3) ^ 2 / 4)
all = 400000
ball = all / 10
FOR ji = 0 TO 8004: JB(ji) = 0: NEXT ji
FOR v = 1 TO all
LOCATE 2, 50
PRINT USING "##########"; all - v
RANDOMIZE TIMER
S1 = 0
uu = 0: uuu = 0: uuuu = 0
FOR i = 1 TO n
a = SQR(-2 * LOG(RND))
x(i) = 1 * a * COS(2 * pi * RND)
S1 = S1 + x(i)
NEXT i
FOR ii = 1 TO n: d = x(ii) - S1 / n
uu = uu + d * d / n
uuu = uuu + d * d * d / n
uuuu = uuuu + d * d * d * d / n
NEXT ii
S = uuu / (uu ^ (3 / 2))
K = uuuu / (uu * uu)
z = fnJB(S, K, n): z = INT(z * 100 + .5)
IF z > 8000 THEN z = 8000
IF z < 0 THEN z = 0
JB(z) = JB(z) + 1
IF INT(v / ball) <> v / ball THEN GOTO 40
vc(0) = .9 * v: vc(1) = .95 * v: vc(2) = .99 * v
FOR cc = 0 TO 2
S = 0: vv = vc(cc)
FOR jj = 0 TO 8000
S = S + JB(jj)
IF S > vv THEN GOTO 20
NEXT jj
20 REM
LOCATE ju + 2, 20 * cc + 1
aa(ju, cc) = jj / 100
PRINT USING " ##.## #.#### "; aa(ju, cc); S / v;
NEXT cc
PRINT USING "############"; v
ju = ju + 1
40 NEXT v: PRINT
END
.