<div dir="ltr"><div>Bonjour,</div><div><br></div><div>Thank you for you reply. </div><div><br></div><div><p class="MsoNormal"><span lang="EN-US">I calculate the jabobian matrix using   J=[(x_node(2)-x_node(1))/2 0;  0 (y_node(4)-y_node(1))/2]. This is a special case for rectangle element (not applicable for general quadrilateral element). </span></p><p class="MsoNormal"><span lang="EN-US"><u></u> <u></u></span></p><p class="MsoNormal"><span lang="EN-US">If I calculate like this, all the terms will be zeros when the ordering changes. I have talked with Jeremy and I will change this formula to adapt to more general cases. </span></p><p class="MsoNormal"><span lang="EN-US"><br></span></p><p class="MsoNormal"><span lang="EN-US">Li Nan</span></p><p class="MsoNormal"><span lang="EN-US"><u></u> <u></u></span></p></div></div>