function v = exam7f(t,y,Z,state)
xb = -0.427;
a = 0.16;
epsilon = 0.02;
u = 0.5;
tau = 1;
Delta = Z - xb;
v = (-y + pi * (a + epsilon*state - u*sin(Delta)^2)) / tau;