The summation limit should be whatever is required to meet your required accuracy. It is usually best to find that out in other tools (even a spreadsheet!) before worrying about VHDL.
For sufficiently smooth curves, a quadratic interpolation using coefficients from a small lookup table may be accurate enough, (again : evaluate in MatLab, Octave or spreadsheet) and is fairly easy to implement.