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

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

Notes
*****

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

	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])));
g = f + a[5];
h = exp(x*((a[5]/a[6]) + s + 0.008314 * log(2.09*a[7]*a[7])-g)/c);
i = sqrt(1.0 + ( h / (20.25 * a[7]*a[7])));
j = (h / (6.0 * a[7]*a[7])) * (1.0 + i);
k = (h / (6.0 * a[7]*a[7])) * (1.0 - i);
v = pow(j, 0.3333333) + pow(k, 0.3333333);
w = (v * (1.0 - v) / (3.0 - (2.0 * v))) * (3.0 * g * g / (c * x));
answer = d + w + (v * (e-d));