1D Quadratic Element Matlab Code

I have created a Matlab code for a 1d steady state heat conduction problem, using linear or 2 node per element, using 'n' elements. I now have to change this for 1D Steady State using quadratic or 3 nodes per element using 'n' elements. I know what is to be changed and what values are to be changed. What I cannot figure out is how do I change them. I need to change nNodes, element co-ordinate system, Gauss 1 pt integration to Gauss 2 pt integration, shape functions, derivatives of shape functions and the jacobian. The shape functions are to be coded in natural or normalized co-ordinate system i.e using Xi which ranges from -1 to 1, instead of x, which ranges from 0 to length.

[Maryam ZD]

I have the same problem. Can you help me?

please connect to me with this e-mail :[e-mail]