function y = specfun(s)
% Est appele par spec.m via la fonction fzero.m

global a
global b

y = sin (pi*s*b/a) .* sin (pi*s*(b-1/2))...
	+ a * cos (pi*s*(b-1/2)) .* cos (pi*s*b/a) ;
