function z=ex314(s,a,b,n)
    fx=inline(s);
    x=linspace(a,b,n);
    y=fx(x);
    z=trapz(x,y);