带参数的非线性方程组

Question

我正试图在matlab中解决非线性方程组。我正在尝试使用fsolve来解决相同的系统但使用不同的参数。

k1=0.01;
k2=0.05;
k3=10;
k4=100;
Ca0=1;
f2=@(C,tau) [Ca0-C(1)-tau*(k1*C(1)-k2*C(2)+k4*C(1).^2);C(2)-tau*(k1*C(1)-k2*C(2)-k3*C(2))];
t = linspace(0,1);
for i = 1:length(t)
 options=optimset('display','off');
 solution(i) = fsolve(f2,[2 0],options,t(i));
end

然后我需要绘制解决方案，但我不断得到Subscripted assignment dimension mismatch.