data d;
  filename d 'domrep86.dat';
  infile d;
  length distname $ 30;
  input distname & mag prsdvotp prsdsp prdvotp prdsp pldvotp pldsp;
  prsds=round(prsdvotp*mag/100,1);
  prds=round(prdvotp*mag/100,1);
  plds=round(pldvotp*mag/100,1);

  ovotp=100-prsdvotp-prdvotp-pldvotp;
  if ovotp<0 then ovotp=0;
  enpv=(prsdvotp/100)**2 + (prdvotp/100)**2 + (pldvotp/100)**2 +
       (ovotp/100)**2;
  enpv=1/enpv;
run;