100 BEGIN COMMENT PLOT A NORMAL DISTRIBUTION CURVE;
110
120 REAL PROCEDURE NORMAL(X); REAL X;
122    BEGIN NORMAL := EXP(-(X^2/2))/SQRT(2*3.14159265) END;
130 
140 PROCEDURE PLOTIT (T); INTEGER T;
142    BEGIN INTEGER I, J, K;
143       J := ENTIER(T/3); K := T - 3*J;
145       FOR I := 1 STEP 1 UNTIL J DO PRINT ("", "   ", "");
150       IF K = 0 THEN PRINT ("", "*")
155          ELSE IF K = 1 THEN PRINT ("", " *")
160          ELSE PRINT ("", "  *");
170    END PLOTIT PROCEDURE;
180
190 REAL X, Y; INTEGER T;
200 
210 FOR X := -2 STEP .1 UNTIL 2 DO
220    BEGIN Y := NORMAL(X); T := ENTIER(100*Y);
225       PLOTIT(T);
230    END LOOP;
235
240 END PROGRAM;
