#yzThqz.m
function Zq0=yzThqz(baasi0,x,qz,EJ)
# 
#  computed a loading vector for uniformly distributed load for a beam 
#  qz -- uniformly distributed load along the z-axis
#  EJ -- bending stiffness of the element
#  baasi0 - scaling multiplier for the displacements (i = EJ/l)
#
##=========================================================================
## This Program is written by Andres Lahe,    1998-05-20
##                    e-mail: andres.lahe@ttu.ee
## LAST MODIFIED: Andres Lahe,    2009-02-13
## Copyright (c)  2004  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;
#
#Zq0(1,1)=-i0*qx*x^2/(2*EA);  - talal pikijõud puudub
#Zq0(2,1)=i0*qz*x^4/(24*EJ);
#Zq0(3,1)=-i0*qz*x^3/(6*EJ);
#Zq0(4,1)=-qx*x;  - talal pikijõud puudub
#Zq0(5,1)=-qz*x;
#Zq0(6,1)=-qz*x^2/2;
#%%%%%%%%%%%%% Tala %%%%%%%%%%%%%%
Zq0(1,1)=i0*qz*x^4/(24*EJ);
Zq0(2,1)=-i0*qz*x^3/(6*EJ);
Zq0(3,1)=-qz*x;
Zq0(4,1)=-qz*x^2/2;
#
endfunction
