function y=skewtententierperturbed(x,p,Q)
% x : initial condition
% p: control parameter
% N:precision
if x<=p
       y=ceil(Q*(x/p));
    else 
    y=floor(Q*((Q-x)/(Q-p)))+1;
end
