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