You are not saying many about your problem, here general solution. I hate to use 'real' and 'ideal' because any model is always ideal, but i will use standard notation
CV, real = CV, ideal + A
A = IinfiniteV T (d2Tp) dV
"I" is integral and (d2T) denotes second derivative of p respect to T.
Once known CV. real can compute CP. real. Do you need or either equation of state of your model (there are thousand of equations available in engineering literature, Geodome can explain this better than me) or an empirical fit with you system.