回答

まず、球面をどのようにパラメタ化していますか?それがわからないと判断できません。
p(u1,u2)をどのように設定したかによって計算過程も計算結果も異なります。

R1221=h11h22-h12^2
これが正しい表式です。
初めの式のR1212は
R1212=-R1221
ですから、ここがおかしいです。

この回答にコメントする
疑問は解決しましたか?