a=[];b=[];

a=0; b=1;

for m=1:11 do
  long=0;
  n=2^m;
  a(n+1)=b;
  for p=1:n do
    long=long+(a(n+1)-a(1))/n*sqrt(1+(2*a(1)+(2*p+1)*(a(n+1)-a(1))/n)^2);
  end
  disp("n =")
  disp(long)
end