function [integrand] = inteval(t);
global aglob
global bglob
a = aglob;  b = bglob;

dxdt = -a*sin(t);  dydt = b*cos(t);
integrand = sqrt( dxdt.^2 + dydt.^2 );
