options crt;  ? 9e1
in klem;
freq a;
smpl 47,71;
print pk pl pe pm py;
print qk ql qe qm qy;  ? (a)

dot k l e m;
  p.p = p./py;
enddot;
print pkp plp pep pmp;  ? (b)

dot k l e m;
  .y = q./qy;
enddot;
print ky ly ey my;  ? (c)

msd(corr) pkp plp pep pmp ky ly ey my;  ? (d)

ntcost = pk*qk + pl*ql + pe*qe + pm*qm;
nucost = ntcost/qy;
rucost = nucost/py;
print ntcost nucost rucost;
corr rucost year;  ? (e)

qvres = qv - (qy - qe - qm);
qvpvr = qv*pv - (qy*py - qe*pe - qm*pm);
msd qvres qvpvr;  ? (f)