I ran 10000 samples from the provided sample file through MAIN.m in R2010a. I get two artifacts, at around 5700 and 9200 samples. Also, I think the overall noise in the output is too large. There must be a problem with the code (besides its extraordinary slowness!)--is there something different in the fixed point functions in R2010a from R14SP3?