function x=SteilsterAbstieg(A,b,c,epsi)
[n,~]=size(A);
x=rand(n,1);
g=-(A*x+b);
k=0;
while norm(g)>epsi
    t=g'*g/(g'*A*g);
    x=x+t*g;
    g=-(A*x+b);
    
    k=k+1;
    
end
k
end
