can anybody give me a numerical code to calculate the
kp equation?