我一直在努力实现自己IGRF-13地磁模型C代码。我按照下面的方程从网站,我得到正确的结果B_phi和B_r磁场强度组件。但是现在我完全困在B_theta组件需要计算偏导数施密特的归一化相关的勒让德多项式对θ3 b根据方程。[![在这里输入图像描述][1]][1]现在,我该怎么做?递归方程19 a、b和c不产生正确的结果,和没有信息在互联网上。我会很感激如果有人可以提供一个示例代码,计算“相关的勒让德多项式的偏导数”或者提供将允许我计算方程。参考:https://hanspeterschaub.info/Papers/UnderGradStudents/MagneticField.pdf [1]: https://i.stack.imgur.com/AcGwh.png