function kp=kiinto0(x0,al,tol,nmax)
k=0;x=x0;virhe=100;
while ((k<nmax) && (virhe>tol)) 
    y=x+al*((x-1).^5-0.001);
    virhe=norm(y-x)/norm(y);
    k=k+1;
    x=y;
end     
kp=x;    
