/*  program to print a table of the normal c.d.f. */

x     = seqa(0,0.04,5);       @ creates col vector 0, .04,  .. .16     @
y     = seqa(0,0.2,11);       @                    0, .2        2.0    @
mat   = y + x';               @         matrix  0.0 thro' 2.8          @
table = cdfn(mat);            @ calculate table of normal prob.s       @
format 8,4;                   @ field width, precision                 @
print "         " x';          @ blank field, x as row                  @
print y~table;                @ y horizontally concatenated with table @