clear
f=@(y)exp(-y.^2).*sin(y);
vastus1=quadv(f,0,1)
y=0:1e-3:1;
z=f(y);
vastus2=trapz(y,z)