function y=fhngen(x,opt)


opt.a=x(1);
opt.b=x(2);
opt.c=x(3);

odeopt = odeset('AbsTol',1e-2,'RelTol',1e-2);
sol = ode45(@(t,x)(fhn(t,x,opt)),opt.tspan,opt.y0,odeopt);

y = deval(sol,opt.tint);
y = y';

y=y(:);
