第1个回答 2015-12-06
x = input('Please enter x: ');
y = input('Please enter y: ');
z = input('Please enter z: ');
e = input('Please enter e: ');
a = input('Please enter a: ');
n = input('Please enter n: ');
B0 = atan(z*(1 + e^2)/sqrt(x^2 + y^2));
Bi = B0;
for i = 1:n
temp = atan(z/sqrt(x^2 + y^2)*(1 + a*e^2*sin(Bi)/...
(z*(sqrt(1 - e^2*sin(Bi)^2)))));
B(i) = temp;
Bi = temp;
end
fprintf('Bn is: %f',B(n))