## SisejoudWGmxpunktis.m
function SF=SisejoudWGmxpunktis(l,X,a,AlgP,mx,GIt,EIw,suurused)

## The forces of beam, ElementNr, at x = X
#
##=========================================================================
## This Program is written by Andres Lahe,   2015-03-30
##                    e-mail: andres.lahe@ttu.ee
## LAST MODIFIED: Andres Lahe,    2015-05-05
## Copyright (c)  2015  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
##=========================================================================
#
#if nargin != 8
if ~(nargin==8)
error(' function SisejoudTalaPunktis have wrong number of input arguments!') 
endif
#%%%%%%%%%%%%%%%%%
# global AP
baasi0=1.0;
xx=X;
#a=X;
#%%%%%%%%%%
vvF=zeros(4,8);
#APt=AlgP;
#%%%%%%%%%%%
 vvF=ylWGfhlin(baasi0,l,xx,GIt,EIw);
Fzx=yzWGmx(baasi0,l,xx,a,mx,GIt,EIw);
Fvv(:,1)=vvF*AlgP+Fzx; #%+vvFz;
Fvv(5,1)=Fvv(2,1)+Fvv(4,1);
#%  %9.2f %9.2e
disp('  ')
#disp(sprintf('%15s %9s  %4.10f   ', '','  x =', X))
disp(sprintf('%8s      %4.10f', ' x =', X))
#
#
for i=1:5 
disp(sprintf('%8s %12.3e ',suurused(i,:), Fvv(i,1)))
endfor
 # 
disp('  ')  
#
#disp('------------------')
#
SF=Fvv(:,1);
#
endfunction
