% Function Kinet(VX, VY, VZ, N) 
% Calculates the Kinetic Energy. It accepts as arguments the atomic velocity vectors
% for the X, Y, Z directions and N = number of atoms, calculates and returns the KE value.
function KE = Kinet(VX, VY, VZ, N);
KE = 0.0;
for i = 1:N
  KE = KE + VX(i,1)^2 + VY(i,1)^2 + VZ(i,1)^2;
end;
KE = KE * 0.5;
