这个文件夹包含了我为本科毕业论文编写的代码。'schnackenberg_final.edp'是一个freefem ++文件，包含了Schnackenberg反应扩散系统，使用分数步长法进行时间前进。随后，它利用时间推进到最终稳态，作为牛顿-拉夫森迭代的初始条件。解决方案迅速收敛，确保达到了稳态。我编写了一些MATLAB脚本，以防止对来自'schnackenberg_final.edp'的Newton-Raphson迭代的Jacobian矩阵进行对角线化。 J用于牛顿-拉夫逊迭代的雅可比矩阵。 JStar在对称适应基础上的雅可比行列式。 RMatrix建立一个矩阵R，使得JSTAR = R”ĴR。为了建立R，即转换为对称适应的基础，我们需要运用一些组表示理论。有关详细信息，请参阅完整文档。如果您对全局有兴趣，请阅读文件“在存在对称性的情况下的数值连续和分叉（2014）.pdf”。这些内容在班加罗尔TIFR-CAM的“2014年有限元会议计算PDE会议”上做了介绍。