The calculation of voltages and currents is quite basic, it is a simple mathematical calculation
.
Starting from the Shockley diode equation
$$I_d=I_s (e^{\frac{V_d}{nV_T}} -1)$$
which is more than sufficient for our calculations, we can start by removing the -1.
This value is only of interest for reverse biased diodes (Vd<0). For forward biased diode we have
$$I_d \approx I_s e^{\frac{V_d}{nV_T}}$$.
However, we do not know the value of n (quality factor), but we do know the characteristics of the diode: Is = 93.2pA and Vd = 72 for Id = 1A.
So can calculate n :
$$V_d \approx nV_T \ln{(\frac{I_d}{I_s})}$$
and therefore nVt = 3.117386 and n = 120.58
We can check this value with the first diodes: V = 63.527, which gives 66mA
Now let's use ohm law $$V_{d1}=R*I_{d2} + V_{d2} \;\;\;\;\; V_d = nV_T \ln{(\frac{I_d}{I_s})} \;\;\;\;\; I = I_{d1} + I_{d2}$$
therefore
$$nV_T \ln{(\frac{I_{d1}}{I_s})} = R*I_{d2} +nV_T \ln{(\frac{I_{d2}}{I_s})}$$
which gives
$$I= I_{d2} + I_{d2} e^{\frac{R*I_{d2}}{nV_T}}$$
By using tools such Scilab, Matlab, ... we easily obtain the solution
Example in scilab, for R = 440 and I = 78mA
deff ("[y] = f (x)", "y = x + x * exp ((440 * x) /3.117386) -0.078")
y = fsolve (0.040, ff)
gives 12.04mA and to get the voltage, just use previous equations
Now, you just have to test different values of R to get the desired current.