function [v,dv] = c_5_p_1_f_n_fp(t,x) 
global k; 
v = -k * x;
dv = -k * v;