我在解这样一个方程组：
dy/dt=27.23*erfc(y/sqrt(4*(t-0.006)))-0.59/(x-y)
dx/dt=7.5/(x-y)-54.46*erfc(y/sqrt(4*(t-0.006)))+1.18/(x-y)
关于求解区间[t0,tf]，t0=0.006，目前我还没找到方法确定tf的值，只知道x(tf)=1   
初值y(0.006)=0  x(0.006)=0.022
我最终是想得到y随t变化的曲线，
有人说用龙格库塔方法也许可以解，我想咨询大家不确定tf的值是否可以求解呢？还是说我必须要给定一个tf值呢
本人matlab初学者，希望大家不吝赐教~