Como transformar um aerofólio NACA em um círculo próximo usando o inverso Karman-Trefftz?

Existe alguma maneira de encontrar o centro de curvatura de ponta de um aerofólio NACA?

Estou tentando transformar um aerofólio NACA em um círculo próximo usando a transformação inversa de Karman-Trefftz. a transformação tem 2 pontos de singularidade, s1 que está localizado no bordo de fuga, e s0 está localizado a meio caminho entre o nariz da borda de ataque e seu centro de curvatura .. agora para aerofólios simétricos é fácil encontrar s0, pois o centro de curvatura está ligado o eixo xe a fórmula do raio da borda principal é conhecida, mas para aerofólios curvos é um pouco complicado com o MATLAB, eu fiz todo o trabalho duro, e agora preso aqui.