lh = [0:.029:1.47]';
h  = exp(lh);
k  = h*.36/.64;

y  = [k;h];

options = foptions;
options(1)  = 1;
options(13) = 2;

st = input('start? ');

if st==1;
  x = constr('rody',y,options);
end;
