2019/11/28 Initialize the package. 
2024/03/12 Convert all Fortran codes to c.
2025/08/27 Adding a function "blinding" to find the changing point (used in effect estimation based on blinding data from clinic trials).
