这是我的集成功能：
function dy_dt=myode(t,y)
dy_dt=(-2*y+3);
这是我的主要剧本：
time_period=[0 10];
IC=[21,12];
[t,y]=ode45(@myode,0:0.1:10,IC);
这个程序运行正常。我给出的第一个初始条件。
（注意：我有一组y值（yref）作为数组）
在第一次迭代结束时，会创建一个y值。一开始，y(1)=yref(1)。在下一次迭代中，我想检查y(2)和yref(2)之间的错误。同样，对于每次迭代，都必须这样做。我怎么能这样做？