Update

1. Modified the equation with scipy