scquad.m

          % essai de quad;  sc_quad.m
          f = inline(' cos(x.*sin (u))./pi', 'u', 'x');
          x = linspace(0,10,100);
          for i = 1:100
            y(i) = quad(f, 0, pi, [], [], x(i));
          end;
          plot(x,y);