function [Y,N] = qiumin
for n = 1:10000
y=f(n);
if y>3
N = n;
Y = y;
return;
end
end
function y = f(n)
y=0;
for i = 1:n
y = y+1/(2*i-1);
end