PC使えるとこに移動しました。

n_abc=(x2-x1)(x3-x1)+(y2-y1)(y3-y1)+(z2-x1)(z3-x1)

n_bcd=(x3-x2)(x4-x2)+(y3-y2)(y4-y2)+(z3-z2)(z4-z2)

|n_abc|=[(x2-x1)^2+(y2-y1)^2+(z2-z1)^2]^(1/2)+[(x3-x1)^2+(y3-y1)^2+(z3-z1)^2]^(1/2)

|n_bcd|=[(x3-x2)^2+(y3-y2)^2+(z3-z2)^2]^(1/2)+[(x4-x2)^2+(y4-y2)^2+(z4-z2)^2]^(1/2)

n_abc*n_bcd=|n_abc||n_bcd|cosA

A=cos^(-1)(n_abc*n_bcd)/[|n_abc||n_bcd|]

であったいるでしょうか？