Add some comments

rel_qkd/key_rates_qber.jl

@@ -317,11 +317,21 @@ function HAgB(alpha, q; eta)
 end
 
 
+"""
+    key_rate(α, q; eta)
+
+The key rate for a given alpha, qber `q` and transmittance `eta`.
+"""
 function key_rate(α, q; eta)
     return compute_entropy(α, q; eta) - HAgB(α, q; eta)
 end
 
 
+"""
+    key_rate(q; eta)
+
+The key rate for a given qber `q` and transmittance `eta`, optimized over alpha.
+"""
 function key_rate(q; eta)
     res = Optim.optimize(alpha -> -key_rate(alpha, q; eta), 0.0, 0.8)
     @show res