Sorry for my poor english, I hope you can understand this post

Given the internal joint angles (joint states) and velocities. I want to compute the acceleration on an end effector of a floating base robot (the acceleration of a human foot in its ground contacts points) using the formula:
J_dot * state_dot + J * state_dot_dot
I need to compute both jacobian and its derivative. But I have problems to compute J_dot cortrectly. Someone colud help me? Any snapshot of code will be very helpfull.
Thanks
juanal