load stack_x.dat

load stack_y.dat



[mb,s]=bay_reg(stack_y,stack_x);



[MB,S,LAM]=rreg(stack_y,stack_x,10);





