PyDrake: Extracting AutoDiff gradients from a robot (controlled via Differential Inverse Kinematics)

Question

I am using PyDrake do build a simple model of a Franka Emika Panda robot arm which picks up and places a brick.

I would like to observe how a change in the initial ch