
trials = 15;

u = zeros(trials,1);

for k = 1:trials
    N = 2^k;
    ff = @(x) sin(N*pi*x).* sin(N*pi*x);
    k
    u(k,1) = quad(ff,0,1);
end

