Матрица, SVD, сингулярное число - решение?

//Задание

2. Задать исходную матрицу А размерности (3 х 4).
3. Используя команду для сингулярного разложения MATLAB, получить представление для матрицы А через сингулярные числа, правые и левые сингулярные векторы в покомпонентной (в виде 2, 3) и в векторно-матричной формах (в виде 1).
4. Проверить условия ортогональности для вышеперечисленных форм представления.
5. Реализовать итеративную процедуру (4) вычисления максимального сингулярного числа и соответствующих ему правого и левого сингулярных векторов при произвольно заданных начальных значениях левого сингулярного вектора U0 и правого сингулярного вектора V0. Вычислить матричную компоненту, соответствующую найденным максимальному сингулярному числу и соответствующим ему правому и левому сингулярным векторам.
6. Используя процедуру метода исчерпывания, получить матричную невязку , для которой выполнить все перечисленные операции пункта 4.

//Мое видение решения

>> A=9*rand(3,4)

//9*rand - по указаниям было сделано, особой технической смазки общей картины не несет)

A =

6.7614 6.2917 4.9249 2.3176
2.2959 8.0181 1.2476 7.5665
4.5536 8.6336 1.3436 2.2885

>> [U S V]=svd(A)

U =

-0.5684 -0.6690 -0.4789
-0.6004 0.7352 -0.3145
-0.5625 -0.1088 0.8196


S =

17.4835 0 0 0
0 5.6734 0 0
0 0 2.7845 0


V =

-0.4452 -0.5871 -0.0820 0.6711
-0.7577 0.1316 0.5534 -0.3199
-0.2462 -0.4449 -0.5925 -0.6249
-0.4088 0.6634 -0.5796 0.2384

>> U*U'

ans =

1.0000 -0.0000 0.0000
-0.0000 1.0000 0.0000
0.0000 0.0000 1.0000

>> V*V'

ans =

1.0000 -0.0000 -0.0000 -0.0000
-0.0000 1.0000 0.0000 -0.0000
-0.0000 0.0000 1.0000 -0.0000
-0.0000 -0.0000 -0.0000 1.0000

>> A1=U(:,1)*S(1,1)*V(:,1)'

A1 =

4.4235 7.5291 2.4463 4.0625
4.6732 7.9540 2.5843 4.2918
4.3782 7.4519 2.4212 4.0209

>> A2=U(:,2)*S(2,2)*V(:,2)'

A2 =

2.2286 -0.4994 1.6885 -2.5179
-2.4491 0.5488 -1.8556 2.7671
0.3625 -0.0812 0.2746 -0.4095

>> A2=U(:,3)*S(3,3)*V(:,3)'

A2 =

0.1093 -0.7380 0.7901 0.7730
0.0718 -0.4846 0.5189 0.5076
-0.1870 1.2630 -1.3522 -1.3228

>> U0=ones(1,3)

U0 =

1 1 1

>> V0=ones(1,4)

V0 =

1 1 1 1

>> V1=UO*A
??? Undefined function or variable 'UO'.

>> V1=U0*A

V1 =

13.6109 22.9434 7.5162 12.1726

>> V1=V1/norm(V1)

V1 =

0.4496 0.7579 0.2483 0.4021

>> U1=A*V1'

U1 =

9.9637
10.4620
9.8451

>> U1=U1/norm(U1)

U1 =

0.5699
0.5984
0.5631

>> S0=U0*A*V1'

S0 =

30.2708

>> abs(S1-S0)

ans =

30.9322

Дело в том, что препод говорит о неверности данного решения, я же в умах разбился и не могу понять где ошибка, подскажите момент, который упустил или реализовал неверно, либо подскажите куда копать..

0