function [x] = thomas(alp,bet,gam,b)
n=length(b);
x=zeros(n,1);
% a=zeros(n);
% for k=1:n
%     a(k,k)=alp(k);
% end
% for k=1:n-1
%     a(k,k+1)=bet(k);
%     a(k+1,k)=gam(k);
% end
for k=1:n-1
    m=gam(k)/alp(k);
    alp(k+1)=alp(k+1)-m*bet(k);
    b(k+1)=b(k+1)-m*b(k);
    gam(k)=m;
end
x(n)=b(n)/alp(n);
for k=n-1:-1:1
    x(k)=(b(k)-bet(k)*x(k+1))/alp(k);
end
% l=eye(n);
% for k=1:n-1
% l(k+1,k)=gam(k);
% end
% u=zeros(n);
% for k=1:n
%     u(k,k)=alp(k);
% end
% for k=1:n-1
%     u(k,k+1)=bet(k);
% end
% norm(a-l*u)

end