实验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;%不会写函数  enddisp(a);m(1)=0.00005;for n = 2:8  m(n) = n*m(n-1);  enddisp(m); r(1) = 0.632;for n = 2:8  r(n) = 1-n*r(n-1);  r(n) = round(r(n)*1000)/1000;  enddisp(r); for n = 1:8  rr(n) = r(n)- a(n);enddisp(rr); plot(rr);%hold on;%plot(m);