function spF=yspSTlfhlin(x)
# sparse fransfer matrix of the of the first order analysis 
# for a Gerber's beam (Sign Convention 2)
#  Zx=F*Z0+lq  
#  Z0'=[Q0 M0]
#
#
##=========================================================================
## This Program is written by Andres Lahe,   2009-02-14
##                    e-mail: andres.lahe@ttu.ee
## LAST MODIFIED: Andres Lahe,    2010-03-12
## Copyright (c)  2009  by Tallinn University of Technology 
##                Department of Mechanics
##                http://www.ttu.ee/
##
## This program is free software; you can redistribute it and/or
## modify it under the terms of the GNU General Public License
## as published by the Free Software Foundation; either version 2
## of the License, or (at your option) any later version.

## This program is distributed in the hope that it will be useful,
## but WITHOUT ANY WARRANTY; without even the implied warranty of
## MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the
## GNU General Public License for more details.

## You should have received a copy of the GNU General Public License
## along with this program; if not, write to the Free Software
## Foundation, Inc., 51 Franklin Street, Fifth Floor, Boston, MA  02110-1301, USA.
## http://www.gnu.org/licenses/old-licenses/gpl-2.0.html
##=========================================================================
#
#
#%  i0=baasi0;
# F=eye(6);
# F(1,4)=-i0*x/EA;
# F(2,3)=-x;
# F(2,5)=i0*x^3/(6*EJ);
# -i0*x/GAr;
# F(2,6)=i0*x^2/(2*EJ);
# F(3,5)=-i0*x^2/(2*EJ);
# F(3,6)=-i0*x/EJ;
# F(4,4)=-1;
# F(5,5)=-1;
# F(6,5)=-x;
# F(6,6)=-1;
#  Fsisse=[1  1  1;      #    Fsisse=[1  1  1;
#2  2  1;                #  2  2  1;
#3  3  1;                #  3  3  1;
#1  4  -i0*x/EA;         #  1  4  -i0*x/EA;
#2  3  -x;               #  2  3  -x;
#2  5  i0*x^3/(6*EJ);    #  2  5  i0*x^3/(6*EJ);
#2  6  i0*x^2/(2*EJ);    #  2  6  i0*x^2/(2*EJ);
#3  5  -i0*x^2/(2*EJ);   #  3  5  -i0*x^2/(2*EJ);
#3  6  -i0*x/EJ;         #  3  6  -i0*x/EJ;
#Fsisse=[1  1  -1;       #%  4  4  -1;
Fsisse=[1  1  -1;        #%  5  5  -1;
2  1  -x;                #%  6  5  -x;
2  2  -1];               #%  6  6  -1];
spF=spconvert(Fsisse);
#
endfunction

