求30句左右的matlab程序

function x=followup(A,b)
n = rank(A);
for(i=1:n)
if(A(i,i)==0)
disp('Error: 对角有元素为0！');
return;
end
end;

d = ones(n,1);
a = ones(n-1,1);
c = ones(n-1);

for(i=1:n-1)
a(i,1)=A(i+1,i);
c(i,1)=A(i,i+1);
d(i,1)=A(i,i);
end
d(n,1) = A(n,n);

for(i=2:n)
d(i,1)=d(i,1) - (a(i-1,1)/d(i-1,1))*c(i-1,1);
b(i,1)=b(i,1) - (a(i-1,1)/d(i-1,1))*b(i-1,1);
end
x(n,1) = b(n,1)/d(n,1);

for(i=(n-1):-1:1)
x(i,1) = (b(i,1)-c(i,1)*x(i+1,1))/d(i,1);
end

追赶法解线性方程组。。。

clear
clc
X=imread('1.tif'); %你电脑上要有这两个名称的对应格式的图片就好了
Y=imread('5.tif');
X=imresize(X,[256 256]);
Y=imresize(Y,[256 256]);
mse=0;
mse=mse+sum(sum(X.*Y));
mse=mse/(256*256);
psnr=10*log(255*255/mse)/log(10);
t(1)=0;
for i=1:256
t(i)=0;
end
for i