请问哪位大神知道上面这个迭代公式用c怎么编写。

j是一种工况，你们可以不管。
双竖线是欧几里得范数。beta是个常量，u是二维空间点的坐标。函数g是个隐函数，请问编写这个必须要函数g的表达式吗？
这个就是1阶偏导，就是梯度。

rho是一个大于等于0小于等于1的变量，g的显示我还没推导出来，可以假设一个简单的g函数，我就是想知道偏导，范数，迭代用c怎么表达！麻烦大家指导指导！

首先谢谢你的关注，这个是一个拓扑优化响应计算中的一个式子，因为现有软件只能用c和fortran语言，老师让我用c，所以没得选择，你提到的gsl，我先去补下知识，不懂得再请教，感谢