clear
g=@(y)exp(-y^2)*sin(y);
vastus=quadv(g,0,1)