f=inline('exp(-x.^2/2)/sqrt(2*pi)');
quadl(f,-1,1)
quadl(f,-2,2)
quadl(f,-3,3)
quadl(f,-4,4)
quadl('sin',0,pi)
help quadl
quit
