файлы матлаб / f123

clear, clc
rng(08), A1= rand(4);
Alpha = 1;
f1 = [1 2 3 4];
f1 = f1';
A=A1'*A1;
f=A1'*f1;
B = eye(4);
y0 = [1; 1; 1; 1];
r0=A*y0-f;
w0=B\r0;
C1=(A*w0)';
tau1 = (C1*r0)/(C1*w0)
y1=B\(B-tau1*A)*y0+B\(tau1*f);
n=0;
while norm (r0)>0.001
r1=A*y1-f;
w1=B\r1;
C=(A*w1)';
tau2 = (C*r1)/(C*w1)
Alpha = inv((1-((tau2/tau1)*((C*r1)/(C1*r0))*(1/Alpha))))
y2=B\((Alpha*B-Alpha*tau2*A)*y1)+B\((1-Alpha)*B*y0)+B\(Alpha*tau1*f)
tau1=tau2;
r0=r1;
w0=w1;
y0=y1;
y1=y2;
C1=C;
n=n+1
end