实验0.2

吉林大学 计算方法实验课 实验0.2

代码用octave语言编写。

题目：

实验代码：

a = [1,2,3,4,5,6,7,8];%e矩阵
r = [1,2,3,4,5,6,7,8];%给的矩阵
m = [1,2,3,4,5,6,7,8];%误差矩阵
rr = [1,2,3,4,5,6,7,8];%自己算的误差矩阵

a(1) = 1- exp(-1);
for n = 2:8
  a(n) = 1- n*a(n-1);
  %a(n) = round(a(n)*1000)/1000;%不会写函数
  end
disp(a);

m(1)=0.00005;
for n = 2:8
  m(n) = n*m(n-1);
  end
disp(m);
 
r(1) = 0.632;
for n = 2:8
  r(n) = 1-n*r(n-1);
  r(n) = round(r(n)*1000)/1000;
  end
disp(r); 

for n = 1:8
  rr(n) = r(n)- a(n);
end
disp(rr);
 
plot(rr);
%hold on;
%plot(m);

实验现象：