数値データを用いて偏微分する方法

MATLABで「変形勾配テンソル」を算出する方法を教えてください．

%            
%              ∂x/∂X   ∂x/∂Y  ∂x/∂Z
%            
%       F = [  ∂y/∂X   ∂y/∂Y  ∂y/∂Z  ]
%            
%              ∂z/∂X   ∂z/∂Y  ∂z/∂Z
%       

物体上に存在する複数の物質点の三次元座標データがあります．

物体を変形させた際の物質点の三次元座標を経時的に記録し，それらのデータから各物質点における変形勾配テンソルを算出したいと考えています．

変形勾配テンソルＦは各物質点の「基準配置Ｘ,Y, Z」と「現在配置 x, y, z」を用いて上記のように算出すると思うのですが，

①MATLABで偏微分を行う方法

②数値データから偏微分を行う方法

が分かりません．

どなたか教えていただけますと幸いです．