Get rid of Tables

From: Luis A. Afonso (licas__at_hotmail.pt)
Date: 02/07/05 

Date: Mon, 7 Feb 2005 19:51:33 +0000 (UTC)

 Get rid of Tables (a quite presumptuous title)

Do use a program that evaluate de Cumulative Normal Distribution

Get exact 7 decimal places values: you can input z as rigorous you
want: z=1.335628
(not the 2 decimal places of the ordinary Tables of textbooks). ONLY
POSITIVE z ....
IF NEGATIVE DROP THE SIGNAL (and do think!)

        Provides de Integrals: Example: z=1.2345678
                        
           from 0 to z 0.3915043
           from -z to z 0.7830086
               < = z 0.8915043
            right tail 0.1084957
            two tails 0.2169914

                
        REM "z1"
        DEFDBL A-Z
17 REM
        PRINT "***** NORMAL LAW *****"
        pi = 4 * ATN(1): c = 1 / SQR(2 * pi)
        INPUT " rigt limit (z>0 but not >6)"; x
        IF x <= 0 OR x > 6 THEN GOTO 17
 DEF fng (x, j) = -x ^ 2 * (2 * j + 1) / ((j + 1) * (2 * j + 3)) * .5
        s = c * x: before = c * x
        FOR j = 0 TO 100
        xx = before * fng(x, j)
        s = s + xx
        before = xx
        IF ABS(xx) < 5E-10 THEN GOTO 100
        NEXT j
100 PRINT
PRINT " from 0 to z "; : PRINT USING " #.####### "; s
PRINT " from -z to z "; : PRINT USING " #.####### "; 2 * s
PRINT " < = z "; : PRINT USING " #.####### "; .5 + s
PRINT " right tail "; : PRINT USING " #.####### "; .5 - s
PRINT " two tails "; : PRINT USING " #.####### "; 1 - 2 * s
        PRINT
        PRINT " do you want to proceed ?"
        PRINT " if YES tape ENTER "
        PRINT " if NOT any other "
        INPUT a$
        IF a$ = "" THEN GOTO 17
        CLS
        END