function ydot=myfunc7(t,y)
 ydot=zeros(1,1);
 tau=1000 * 1.e-6;
 ydot(1)=(1-y(1))/tau;
endfunction
