function f=fa(t,z)
global c q m;

f(1)=c*z(1);
end