(*   Kvadraten koren Vavezdat se samo elementite po i pod gl. diagonal; ostanalite - празни   *)
A=(({{4, \[Placeholder], \[Placeholder], \[Placeholder]}, {2, 2, \[Placeholder], \[Placeholder]}, {0, 2, 5, \[Placeholder]}, {2, 1, 1, 3}}));b={4,-1,-6,3};
n=Length[b];
For[i=1,i≤n,i++,(Subscript[A, [[i,i]]])=(Sqrt[(Subscript[A, [[i,i]]])-(∑_(k=1)^(i-1))((Subscript[A, [[i,k]]])^2)]);
For[j=i+1,j≤n,j++,(Subscript[A, [[j,i]]])=(((Subscript[A, [[j,i]]])-(∑_(k=1)^(i-1))(Subscript[A, [[i,k]]])(Subscript[A, [[j,k]]]))/(Subscript[A, [[i,i]]]))]]
Print["    T = ",MatrixForm[A]]
For[i=1,i≤n,i++,(Subscript[b, [[i]]])=(((Subscript[b, [[i]]])-(∑_(k=1)^(i-1))(Subscript[A, [[i,k]]])(Subscript[b, [[k]]]))/(Subscript[A, [[i,i]]]))]
Print["   Prav hod T*y=b   y = ",MatrixForm[b]]
For[i=n,i≥1,i--,(Subscript[b, [[i]]])=(((Subscript[b, [[i]]])-(∑_(k=i+1)^n)(Subscript[A, [[k,i]]])(Subscript[b, [[k]]]))/(Subscript[A, [[i,i]]]))]
Print["   Obraten ho T'*x=y   x = ",MatrixForm[b]]
Print["   Det(A) = ",(∏_(i=1)^n)((Subscript[A, [[i,i]]])^2)]