function X=stepp(iv,T)
X = 1:(T+1);
X(1) = iv;
U = rand(size(X));
for t=2:(T+1)
    X(t) = updatep(X(t-1),U(t));
end