fun = @(x) sqrt(1-power(x,2));
format long
res=integral(fun,0,1)*4