optimization - 使用 octave 中的优化包求解非线性方程组的约束系统

Question

我已经安装了 Optim1.6.1 包。我想在某些节点处使用约束为 u=1 来求解非线性有限元分析中的方程组。u=0 在某些节点。我不清楚如何输入刚度矩阵、残差矩阵，因为这些矩阵是在单独的脚本中作为多个变量的函数计算的。通常我发现 kmat 和 res 如下

[kmat,res]=systemmatrices(u,load,boundarycondition,material,constraints)

我能找到一个简单的例子吗？