The computational investigation of the correlation between the achievable reaction yield and discharge energy for a plasma focus device (PFD) is presented. Radioisotope production in PFDs with applicable activities is highly dependent on establishing the related scaling law. Carbon target is bombarded by high energy deuterons and short-lived radioisotope of 13N is produced through 12C(d,n)13N in which the threshold energy is not very high. Both computed and measured ion energy spectra are used to estimate and optimize the scaling law. It is shown that the number of ions emitted from the pinch region for a device operating under optimized conditions is linearly proportional to the discharge energy of the PFD.