use quadratic polynomial with a=0.001