Goldschmidt Solution

import pylab
import numpy
import math
a=1.5089
k=1.2
a*a
k*k
y=2.27677921+1.44
p=k+(1/2)math.sinh(2*k)
s=math.sinh(2*k)
t=k+(1/2)*s
b=y/t
print(b)