我想用 Matlabpdepe来解决这个系统：

$$s_t =(sr)_x + s_{ xx } \\ r_t =(\frac{ A }{ B }r^2+s)_x + \frac{ A }{ -K } r_{ xx }$$

在哪里$A$,$B$和$K$是常数。我需要使用 MATLAB 解决这个问题pdepe。问题是我方程中的非线性项：我不知道如何处理非线性部分。边界和初始条件在该部分不是问题。

因此，我尝试以pdepe期望的形式编写我的 PDE：

$$\left( \begin{array}{c} 1 \\ 1 \\ \end{array}\right) \times \frac{ \partial }{ \partial t }\left( \begin{array}{c} s \\ r \\ \end{array}\right)=\left( \begin{array}{c} s\frac{ \partial r }{ \partial x } + r\frac{ \partial s }{ \partial x } \\ 2\frac{A}{B}r\frac{ \partial r }{ \partial x } +\frac{ \partial s }{ \partial x }\\ \end{array}\right)+\frac{ \partial }{ \partial x }\left( \begin{array}{c} \frac{ \partial s }{ \partial x } \\ \frac{A}{-K}\frac{ \partial r }{ \partial x } \\ \end{array}\right)+\left( \begin{array}{c} 0 \\ 0 \\ \end{array}\right)$$

在 MATLAB 中：

function [c,f,s] = eqtn(x,t,u,DuDx)

c = 1;

f = ????

s = 0;

如何对函数进行编码$f$?