options crt;  ? 7e1
in ukelec;
smpl 1,42;
print mc4 mc6 mc8 gas6 gas8;
msd mc4 mc6 mc8 gas6 gas8;  ? (a)
shelec = expen/inc;
print inc kwh expen cap cust shelec;  ? (b)
avgpe = 240*expen/kwh;
f = expen - mc8*kwh/240;
print avgpe f mc8;  ? (c)
corr kwh inc mc6 gas6 cap;
lnx = log(kwh);
lnm = log(inc);
lnpe2 = log(mc6);
lnpg2 = log(gas6);
lnh = log(cap);
corr lnx lnm lnpe2 lnpg2 lnh;  ? (d)