added pod_J and calculate of Newton pod_du directly