MATLAB程序设计与应用刘卫国(第三版)课后实验答案——5

31
五月
2021

MATLAB程序设计与应用刘卫国(第三版)课后实验答案——5

y=0;
	n=input('n=');
	for i=1:n
  		y=y+(1/(2*i-1))*(1/(3.^(2*i-1)));
	end
y

在这里插入图片描述

n=input('n=');
	i=(1./(2.*(1:n)-1)).*(1./(3.^(2.*(1:n)-1)));
	s=sum(i);
s

在这里插入图片描述
在这里插入图片描述

y=0;
n=0;
while y<3
   n=n+1;
  y=y+(1./(2.*n+1));
end
if y>3
   n=n-1;
end
n
y
i=(1./(2.*(1:n)+1));
y=sum(i)

在这里插入图片描述

a=input('正常数a=');
b=input('正常数b=');
x=1.0;
n=0;
while abs(x-a./(b+x))>=0.00001
   if n<500
      x=a./(b+x);
      n=n+1;
   end
end
x
r(1)=(-b+sqrt(b.^2+4*a))/2
r(2)=(-b-sqrt(b.^2-4*a))/2
s=r-x

(1)
在这里插入图片描述
(2)
在这里插入图片描述
(3)
在这里插入图片描述

f(1)=1;
f(2)=0;
f(3)=1;
for n=4:100
   f(n)=f(n-1)-2*f(n-2)+f(n-3);
end
disp(['max(f)=',num2str(max(f))])
disp(['min(f)=',num2str(min(f))])
disp(['sum(f)=',num2str(sum(f))])
disp(['c1=',num2str(length(find(f>0)))])
disp(['c2=',num2str(length(find(f==0)))])
disp(['c3=',num2str(length(find(f<0)))])

在这里插入图片描述

j=0;
s=0;
for i=2:50
   p=i*(i+1)-1;
   for k=2:p  %此循环判断是否为质数
       if rem(p,k)~=0
           continue
       end
   break
   end
   k;
   if k==p  %除到最后一个数判断完毕
       j=j+1;%亲密数对数
       s=s+p;%亲密素数和
   end
end
j
s

在这里插入图片描述

TAG

网友评论

共有访客发表了评论
请登录后再发布评论,和谐社会,请文明发言,谢谢合作! 立即登录 注册会员