function pi = Pi_matlab( n )

h = 1.0/n;
pi = 0.0;

for i=1:n
  x = (i + 0.5)*h;
  pi = pi+ (4.0/(1.0 + x*x));
  end
 pi=pi*h;
 disp(pi)
end