算法算法算法

算法一

clear
clear
i=1;j=1;
n=1000;
s=0;

for i=1:n
s=s+1/j;
j=i+j;
if (1/j)<(10^(-5))
s
i
break;
end
end
结果如图

算法二
计算椭圆所有点到椭圆焦点的平均距离
代码如下
clc
clear
syms theta;
a=5;
b=3;
c=sqrt(a^2-b2);
e=c/a;
p=abs(a^2/c-c)
result=int(ep/(1-ecos(theta)),0,2pi)/2/pi
theta_div=0:0.01:2pi;
r_div=ep./(1-ecos(theta_div));
plot(theta_div,r_div)
hold on;
plot(theta_div,result)
% hold off;
% ezplot(‘x^2/(532)+ y^2/(32) = 1’,[-53 53 -3 3])
结果如下