format long
f=@fun;
res=integral(f,0,1)*4

function y=fun(x)
y=sqrt(1-x.*x);
return
end
