function demoros(mu,sig)
% just a demo
if nargin<2,
		  sig=[2,-1.8;-1.8,2];
end
if nargin<1,
		  mu=[0,0];
end

N=4000;
rho=0.005;
alpha=.7;
epsilon=.001;
beta=.95;
q=.925;

 ceros(mu,sig,rho,alpha,N,epsilon,beta,q);

