proc nlin  data=data_prg.reg6_7b;
  parms  b1=1 to 2 by 0.005  b2=0.1 to 0.2 by 0.0005;
  model  y=b1*x1**2+b2*x2*x3;
  der.b1=x1**2;
  der.b2=x2*x3;
  output out=ccc p=yhat r=resid sse=scrs l95m=lm u95m=um;
run;
quit;