I have questions about the code (found here Estimation of Risk-Neutral Densities Using Positive Convolution Approximation - Python).
The synthetic price in Bondarenko paper includes two terms before the dot product of ak with double integral. These correspond to the double integral left of the listed strike range. I cannot find them in this code.
Also there is likely a typo in formula Bondarenko gave for double integral he wrote
double integral=density - x * integral of density
I get
double integral = x * integral of density - integral of x* density
the code uses a third expression
double integral = x* integral of density + density
.If you differentiate this or Bondarenko's formula, result doesn't look right. so there may be a typo there as well.