function yp=fopc_f(t,s,flag,tu,name)
% Subroutine for FOPC - fwd integ. sdot=f(s,u); 2/2/98
%
u=interp1(tu(:,1),tu(:,2:end),t);  
yp=feval(name,u,s,t,1); 