求问of编写偏微分方程求解问题

[小泽同学]

在编写激光加工数值仿真的时候遇到了一个偏微分方程的求解，应该是很简单的，但是我一直求不出准确的解，所以来问问是不是编写错误，请各位大佬帮帮忙。
公式如下：

下面是我的代码：


引入一个Nc的三次方项就一直报错，但是用Nc直接乘三遍再乘系数就没事，想问问代码编写的对不对？感谢各位！

[李东岳]

你的方程，除了$n_e$，都是已知量么

[小泽同学]

谢谢大佬回复，是的，除了Nc都是已知的，然后NcCubed就是Nc的三次方，

[李东岳]

solve
(
    fvm::ddt(Nc)
  ==
    ...
  + beta*Nc
  - fvm::Sp(gamma*Nc*Nc, Nc)
);

[小泽同学]

@李东岳 谢谢大佬，我去试试感谢