function mycirc,p, t=t, xval=xval, yval=yval, xerr=xerr, yerr=yerr

	; p = [x0, y0, r, omega, phi]
	
	modelx = p[0] + p[2]*sin(p[3]*t + p[4])
	modely = p[1] + p[2]*cos(p[3]*t + p[4])

	modeldist = sqrt((xval-modelx)^2. + (yval-modely)^2.)
	errs = sqrt(xerr^2. + yerr^2.)

	return,modeldist/errs

end