function [s x] = ARprocess (a,varu,varn,N)
s(1)=randn()*sqrt(varu/(1-a^2));
for k=2:N
    s(k) = a*s(k-1)+sqrt(varu)*randn();
end;
x = s + sqrt(varn)*randn(size(s));