function[x,k]=cg(A,b,x0,maxit,tol)
% CG zum Loesen der Gleichung Ax+b=0
% Algorithmus gemaess Skript
x=...%%%
r0=...%%%
r=...%%%
p=...%%%
for k=1:maxit
    ...%%% Tipp: speichere das Skalarprodukt vom Residuum mit sich
        % selbst in einer Hilfsvariable
    if k>=2
        ...%%%
    end
    ...%%% Tipp: Speichere Matrix*Vektor in Hilfsvariable    
    if ...%%%
        return
    end
end

if k==maxit
    error(['keine Konvergenz nach ',int2str(k),' Schritte'])
end
