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

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

Notes
*****

	model to fit the partial molar heat capacity for the
	unfolding of a dimeric 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, in monomer units)
	a6= Tmax (K)
	a7= conc (M, monomer units)

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

c = x * 0.008314;
d = x * x * a[8] + 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[8])) * ((x*x) - (a[6]*a[6])));
s = 2.0 * s;
g = 2.0 * (f + a[5]);
h = exp((x * 2.0 * ((a[5] / a[6]) + s + 0.008314 * log(1.66 * a[7])) - g) / c);
v = ((-1.0 * h + sqrt(h * h + 8.0 * a[7] * h)) / 4.0 * a[7]);
w = (v * (1.0 - v) / (2.0 - v)) * (2.0 * g * g / (c * x));
answer = d + w + (v * (e-d));