T =[535 550 565 575 585 595 605 615 625]'; H0C= -226000; H0B = -123000; H0A = -66600; CpC = 46; CpB = 18; CpA = 35; CpM = 19.5; TR = 528; Ti0 = 535; vA0 = 46.62; vB0 = 46.62; VM0 = 233.1; V = 40.1; FA0 =43.04; FM0 = 71.87;; FB0 = 802.8; A = 16.96e12; E = 32400; R = 1.987;