turn GSL error handler off for numerical integration to allow it to continue even if integration doesn't coverge (also allowed by 'stop.on.error=FALSE' in R version)
added anclik.R: wrapper function for calling C likelihood and preparing ANC data
Updated header to describe changes to algorithm, omitted not yet optimised code
added an example of likelihood calculation using Botswana Urban ANC data