请教周期激励（如正弦）微分方程系统的分岔图的计算程序实现问题

zearo

求jacobian矩阵没有问题了，matlab里面就有这个命令jacobian。
又遇到新的问题的，请帮我看看吧
公式(3.2)第四个式子，也就是phi对u和lamda的二阶倒数中的最后一项 按照解析的公式是没有错的，但怎么发现维数对不上呢？
用P表示偏导的意思，2表示2阶倒数，（3.2）中第四个式子的最后一项为
P2fPxPlamda*PphiPu，但是 P2fPxPlamda是2×4矩阵，而PphiPu是2×2矩阵，这个乘法无法计算啊，我在那里理解错了么？