Function Description
********************

	Name: Z_CP_MONOMER "Partial Molar Heat capacity of monomer"
	Function: Cp(temp) for unfolding of monomer 
	Parameters:  a0 a1 a2 a3 a4 a5 a6 a7 a8

Notes
*****

	model to fit the partial molar heat capacity for the
	unfolding of a monomeric protein

	see Yu et al. 1999 J. Phys. Chem. 103, 2270-

	NOTE: The temperature must be in K

	a0= first order parameter for Cp,N (kJ/K^2 mol) 
	a1= zero order parameter for Cp,N (kJ/K mol)
	a8= second order parameter for Cp,N (kJ/K^3 mol)
	a2= second order parameter for Cp,U (kJ/K^3 mol)
	a3= first order parameter for Cp,U (kJ/K^2 mol)
	a4= zero order parameter for Cp,U (kJ/K mol)
	a5= DHu (kJ/mol)
	a6= Tmax (K)
	a7= conc (M)

C-code for the function
**********************

c = x * 0.008314;
d = x * x * a[7] + x * a[0] + a[1];
e = x * x * a[2] + x * a[3] + a[4];
f = ((a[4] - a[1])*(x - a[6])) + ((0.5 * (a[3] - a[0])) * ((x*x) - (a[6] *
	a[6]))) + ((0.3333 * (a[2] - a[8])) * (pow(x,3.0) - pow(a[6],3.0)));
s = ((a[4] - a[1])*(log(x) - log(a[6]))) + ((a[3] - a[0])*(x-a[6]))
	+ ((0.5 * (a[2] - a[7])) * ((x*x) - (a[6]*a[6])));
g = f + a[5];
h = exp((x * ((a[5] / a[6]) + s) - g) / c);
v = h / (1.0 + h);
w = v * (1.0 - v) * (g * g / (c * x));
answer = d + w + (v * f);