% Example 10.3b (Physics 330) f=inline('exp(-x).*cos(x)','x') quadl(f,0,2,1e-8)