x = 0:0.05:2; y = diag(f(x).*f(x+2)); plot(x,y); xlabel('\bfx'); ylabel('\bfy');