function ydot=myfunc7(t,y)
 ydot=zeros(1,1); % —ÌˆæŠm•Û
 tau=1000 * 1.e-6;
 ydot(1)=(1-y(1))/tau;
% end of function
