rng('default')
   b1=[1 1];
   b2=[1 2.6];
   n1=40;
   n2=50;
   s1=0.1;
   s2=0.1;
   X1=rand(n1,1);
   X2=rand(n2,1);
   y1=randn(n1,1)*s1+b1(1)+b1(2)*X1;
   y2=randn(n2,1)*s2+b2(1)+b2(2)*X2;
   hold('on')
   plot(X1,y1,'o');
   plot(X2,y2,'o');
   title('Two simulated regression lines')
   y=[y1;y2];
   X=[X1;X2];
   figure
   [out]=FSRmdrrs(y,X,'constr','','nsimul',1000,'init',10,'plots',1,'cleanpool',false,'msg',0)