I am trying to estimate the cumulative distribution function from a probability distribution function
for this I am using the lambda function (code given below), but it gives some negative values at the end of series
Dj_lin = lambda x: 1 - sp.integrate.quad(mypdf_lin, 0, x)[0]
here mypdf_lin is my pdf function
is there any other way to estimate the cdf from pdf ?