rho = 0.8
k = 0.262

sigma_F_over_sigma_S = (0.03/0.025) * (1/k)
print(sigma_F_over_sigma_S)

B2 = 1 - 2*k*rho*sigma_F_over_sigma_S + k^2 * sigma_F_over_sigma_S^2
B = sqrt(B2)
print(c(B2, B))