function [k]=Bar3Noded_1Dk(X,b,h, E, Nu)
%clear all;                                   


D=E;

syms Xi real;



XiC=[-1 0 1];

Ni=[ ((XiC(2)-Xi)*(XiC(3)-Xi)) /((XiC(2)-XiC(1))*(XiC(3)-XiC(1))), ...
     ((XiC(1)-Xi)*(XiC(3)-Xi)) /((XiC(1)-XiC(2))*(XiC(3)-XiC(2))), ...
     ((XiC(1)-Xi)*(XiC(2)-Xi)) /((XiC(1)-XiC(3))*(XiC(2)-XiC(3))) ];

dNibydXi=diff(Ni,Xi);


x=Ni*X';

J=[ diff(x,Xi)];
Jinv=inv(J);

dNibydx=simplify(Jinv(1,1)*dNibydXi);

B=[dNibydx(1)     dNibydx(2)   dNibydx(3)];    

BTDBJ1=subs(B'*D*B*det(J)*b*h,{Xi},{ -1/sqrt(3)});
BTDBJ2=subs(B'*D*B*det(J)*b*h,{Xi},{  1/sqrt(3)});
k=BTDBJ1+BTDBJ2





   