[求助]有点头大的一阶微分方程的问题

[wolfnet]

方程 dy/dt+i=240*sin(2*pi*50*t)

其中i=(0.0006455*exp(0.4182*(10*y/0.001+21)) + (-3.971e+004)*exp(-0.4854*(10*y/0.001+21)))*0.1/100

我想用ODE45来得到这个微分方程的数值解,其中y(0)=0;
请高手帮忙.
我自己试了多次,总是出错.我不知道问题是不是出在这里.这个i和y的方程是我用CFtoolbox得到的.下面的y和微分方程中的y是同一个参数.在下面的程序里面,我已经设定了y的范围.也就是说如果y超出这个范围,i就不是我所想要的值了.而且y的值不能这样取-0.0015:0.0015,如果这样取的话,也是错误的.
如果问题是在这里的话,应该怎么解决呢?
设定如下:
y=-0.0015:0.00001:0.0015;
i=(0.0006455*exp(0.4182*(10*y/0.001+21)) + (-3.971e+004)*exp(-0.4854*(10*y/0.001+21)))*0.1/100;