ex2.m

重点内容是计算中的误差、函数方程求根、插值呆逼近、数值积分和微分、线性代数方程组解法、常微分方程初值问题的数值解法、数学软件。

n=input('输入矩阵阶数：');
t=cputime;
a=zeros(n,1);b=zeros(n,1);c=zeros(n,1);d=zeros(n,1);
y=zeros(n,1);x=zeros(n,1);
for i=1:n;
    b(i)=2;a(i)=1;c(i)=1;d(i)=-5;d(1)=-7;
end
for j=2:n;
    q(1)=b(1);
    p(j)=a(j)/q(j-1);    
    q(j)=b(j)-p(j)*c(j-1);
end
for j=2:n;
    y(1)=d(1);    
    y(j)=d(j)-p(j)*y(j-1);
end
x(n)=y(n)/q(n);
for k=n-1:-1:1;
    x(k)=(y(k)-c(k)*x(k+1))/q(k);
end
e=cputime-t