<code>매트랩(ode) 도와주세요 ㅠㅠ

안녕하세요. 매트랩 코드 관련 질문입니다.
숫자들이 많지만 실제로 내용은 별로 없어요. 한번만 봐주시면 감사하겠습니다.
실제로 식은
S'=n(H-S)-aBS
I'=aBS-(r+m+n)S
B'=hB+bI
인 ODE 입니다.
S 를 f(1) 로 I 를 f(2) 로 B 를 f(3) 으로 했습니다.

function dfdt = col(f,a,b)
y= [2175,3075,3279,3806,4270,4583,5285,5876,6548,8137,11183,15983,19499,27431,37204,48647,56092,62607,69761,78140,86107];
z=[22,26,27,33,32,33,35,35,35,41,51,60,66,74,85,108,120,131,139,163,181];
H=44000000;
n=-0.2;
r=0.2;
h=-10^(-3);
m=sum(z,1)/sum(y,1);
dfdt=zeros(3,1);
dfdt(1)=n*(H-f(1))-a*f(3)*f(1);
dfdt(2)=a*f(3)*f(1)-(r+m+n)*f(2);
dfdt(3)=h*f(3)+b*f(2);
end

으로 col 이라는 함수를 저장하고
이를 이용하여

a=0.8;
b=10^(-6);
tspan=[2000 86107];
f1_0=44000000; -초기조건들-
f2_0=2175;
f3_0=10000;
f0=[f1_0,f2_0,f3_0];
syms t
[t,f]=ode23(@(t,f) col(f,a,b),tspan,f0);
plot(t,f(:,2),'-')

을 코딩하였는데
syms t에서 오류가 나고
또 밑에서 두번째 줄에서도 오류가 나는데 왜 이런건지 모르겠습니다.ㅠㅠ
고수님들 도와주세요ㅠㅠ