Excelente tutorial.

Porém verifiquei um erro na fórmula:
c = 2 * atan2(sqrt(a), sqrt(1-a))

O correto seria:
c = 2 * atan2(sqrt(1-a), sqrt(a))

Comparei fórmula corrigida com a Lei dos Cossenos e valores bateram.