Puro <- function(theta)
{
	conc <- c(0.02, 0.02, 0.06, 0.06, 0.11, 0.11, 0.22, 0.22, 
		0.5600000000000001, 0.5600000000000001, 1.1, 1.1)
	rate <- c(76, 47, 97, 107, 123, 139, 159, 152, 191, 201, 207, 200)
	denom <- theta[2] + conc
	res <- rate - ((theta[1] * conc)/denom)
	sum(res^2)
}