Matlab报错：矩阵拼接维度不一致，求解带质量矩阵的ODE方程组遇阻

我知道这个问题在论坛里挺常见的，已经翻了好多相关帖子，但到自己这个具体问题上还是卡壳了。我现在正试着用质量矩阵配合ode函数同时求解4个方程——其中2个是代数方程，2个是微分方程，变量定义如下：

• y(1)：半径1
• y(2)：浓度Cc
• y(3)：半径2
• y(4)：浓度Cc2

我的代码片段如下：

function ade 
% Equation set 
% y(1)== radius 1, y(2)== concentration Cc 
% y(3)== radius 2 y(4)== concentration Cc2 
%========================================================================== 
clc 
close all 
M = [1 0 0 0 0 0 0 0 0 0 1 0 0 0 0 0]; 
options = odeset('Mass',M...

运行这段代码时，Matlab抛出了错误：

Dimensions of matrices being concatenated are not consistent.

内容的提问来源于stack exchange，提问作者Dan Thomas