This code fails to run, ("rep" and "expandm" are unknown), has no help, no useful comments, creates a useless plot, and is poorly coded. I'll argue that for now this deserves 1 star, but after proper revision I am happy to change it accordingly.
No help. No explanation of what the variables mean at all. No H1 line. No explanation of what the function does. No error checking. This requires two functions square.m, and expandm.m to do trivial things. square.m comes from the signal processing toolbox, and is trivial to replace. expandm.m comes from where? A search of the website did not find it.
Code that does not run and has absolutely no help at all, should probably get a poor rating. So I'm arguably being remarkably generous by saying this "needs improvement".
Others might say this is worth a better rating. They would be very wrong here. At the very best, you can only argue that this is worth a 2 rating.
