is there any correction..?

I am doing Singular value decomposition using this.

 Xi = [1 2;3 4]; Yi = [1 2;3 4]; Zi = [1 2;3 4];
 A=[Xi,Yi,Zi,ones(length(Xi),1)];
 [U S V] = svd(A)

Now, i want to do the same without function. so, i written little as below.

 Xi = [1 2;3 4]; Yi = [1 2;3 4]; Zi = [1 2;3 4];
 A=[Xi,Yi,Zi,ones(length(Xi),1)];
 f = eig(A'*A);
 for l = 1:length(f);
 p(l) = sqrt(f(l));
 end
 S = diag(p);
 [U,D] = eig(A*A');
 [V,D] = eig(A'*A);
 V = V';
 U
 S

But in answer i am getting only some of its elements same in matrix. And remaining are different. So, please suggest any modification if there or solution..?

Answer 1

[Walter Roberson]

In your second bit of code, your S is a square diagonal matrix, not a diagonal matrix the same size as A.

[Lalit Patil]

Then what change to do in S.?

And before S i want to solve problem of U and V.

What change to do in them.?