Function poly_impulse 

pub fn poly_impulse<T: Number + Float + Base<T> + FloatOps<T>>(
    k: T,
    x: T,
    n: T,
) -> T

returns a quadratic impulse (y position on a graph for x); n is the degree of the polynomial and k controls the stretching of the function