f=@(x,y)1./(sqrt(x+y).*(1+x+y).^2);
ymax=@(x)1-x; %y ülemine raja

vastus=integral2(f,0,1,0,ymax)