algorithm works for squared symmetric matrices - Other operators implemented but not tested.
initial commit for tast 3, coppied from task 3