options crt, limwarn=0; 
freq q; smpl 47:1 95:3;
read(file='ex66.xls', format=excel) gcq gydq;
consum=gcq;
yd=gydq;
ydsq=yd**2;

olsq consum c yd ydsq;

frml eq1 consum=b0+b1*(yd**b2);
mmake @start 1 1 1;
param b0 b1 b2;
ml(maxit=100) eq1;

end;stop;