function T=MatGauss(x,p)
n=length(x);
t=zeros(n,1);
t(p+1:n)=x(p+1:n)./x(p);
I=eye(n);
T=I-t*I(:,p)';