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


syms Xi real;


L= ( (X(end)-X(1))^2 + ( Y(end)-Y(1))^2 )^0.5;

kee=bar2noded_1Dk([0 L],b,h,E);
%kee=b*h*E/L*[1 -1;-1 1]
c=(X(end)-X(1))/L;
s=(Y(end)-Y(1))/L;
k=[c s 0 0;0 0 c s]' *kee* [c s 0 0;0 0 c s];








   