function x = mychol_solver(a,b)
% mychol_solver  Solve linear system a*x=b using Cholesky decomposition
% when is s.p.d.

% Output variables:
% x = a^-1*b, the solution

[r, f] = mychol(a);


y = forbacksubs(r,b,0);
x = forbacksubs(r',y,1);
