options crt;  ? 4e6
in chow;
smpl 1,137;
mem = words*binary*digits;
dot rent mult access add mem;
  ln. = log(.);
enddot;
dot(value=jyr) 55-65;   ? jyr = 55 to 65
  d. = (year = jyr);
  select d. .or. (year = (jyr-1));
  olsq(silent) lnrent c d. lnmult lnmem lnaccess;
   set b. = @coef(2);
  select 1;
enddot;
mmake bm b55-b65;
?
d54 = (year=54);
olsq(silent) lnrent d54-d65 lnmult lnmem lnaccess;
  unmake @coef c1-c15;
  mmake cm c1-c12;
?
freq a;
smpl 54,65;
unmake cm cs;
smpl 55,65;
dc = cs-cs(-1);
unmake bm bs;
bsum = bs;
smpl 56,65;
bsum = bsum(-1) + bs;
smpl 55,65;
print bs dc cs;  ? (a)
prind = exp(cs-c1);
chind = exp(bsum);
print prind chind;  ? (b)