错误使用 plot 数据必须为可转换为双精度值的数值、日期时间、持续时间或数组

这是我的代码，错误处在plot，我想研究m随时间t的变化情况，观察在第几年时m降到m0以下。谢谢万能的网友们的帮助。
m=dsolve('Dm=(m0/6)-m/3','m(0)=5m0')
year=1:20
y=m
plot(year,y)
set(gca,'xtick',(2000:1:2020))
ylim([0 10m0])

无法用你的代码复现，建议查看一下求解出的m的格式，如果m是函数的话建议使用fplot(m)