function [th,p,be]=junk3(x)
z=[]; be=[]; th=[]; theta=-.26;
for i=2:100
    [p,beta,psi,ps]=stableexample(x,theta,1.5,1000);
    theta=theta+(.02/sqrt(i))*p
    z=[z;p]; th=[th theta];
    be=[be;beta];
end
p=mean(p);
be=mean(be);
plotecdf(z)