function x=ausgleich(A,b)

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

x=zeros(3,1);

for i=1:3
   x=x+U(:,i)'*b*V(:,i)/S(i,i);
end

inv(U)*b