to be removed. Do not use.
compute I, D, D_k, R_k