for i=1:9
s1=r(i)^2*4*pi^2*10.4^2*1.02*10^(-16)*(1-exp(-8*x))/(8*2.82*10^(-5)*0.347*0.347*0.0037^2*1.36*10^(-12))-86400*8*x-1000*86400;
s2=-y+0.085+86400*2.82*10^(-5)*0.347^2*0.0037^2*1.36*10^(-12)*exp(8*x)/(4*pi^2*r(i)^2*10.4^2*1.02*10^(-16))-x;
% % % % [x,y]=solve('s1=0','s2=0');
s=solve(s1,s2,x,y);
m(i)=s.y
end
m
