function cumnord=cumnord(x);
   if (x>=0); cumnord=0.5+quad('phi',0,x); end;
   if (x<0); cumnord=0.5-quad('phi',x,0); end;
end;