Obtain the shape function and its derivative at each element

Hello,

I would like to get the shape function (or characteristic function) together with its derivative at each element. Is there a quick way in freefem to do so? Many thanks!

Note: I understand that I can use varf to get the shape function of each node in an integral sense. However, what I need is the shape function of each element separately, without summing up the contribution from surrounding elements. For example, for the following mesh, I would like to get the shape function N1 (and its derivative) of element 0 on node 0, and the shape function N1 (and its derivative) of element 1 on node 0, instead of the sum of two N1 (and their derivatives) from two elements on node 0.