function yp = format_f(t,y)

k = length(y);
yp = zeros(k,1);

yp(1) =         ;       % Should be function of y(1), y(2), ...
yp(2) =         ;
yp(3) =         ;