Во-первых, я надеюсь, что вы понимаете, что таким образом вы считаете кривизну не абы где, а именно на горизонте.
Если вы посмотрите внимательно на ваши формулы, вы поймете, что мнимая единица может появиться ровно в одном месте - в

. Значит нужно понять как так получилось, что

оказалась отрицательной.
Проблему найти легко. Вы работаете в системе единиц, в которой

является размерной величиной, а значит исходя из формулы


должно быть другой размерности в сравнении с

и

. В то же время в формуле для гравитационного радиуса

их размерности совпадают. Значит в зависимости от определения

, которое вам больше нравится, одна из этих формул является неверной. Т.е. если взять

, то на самом деле
