+ int a,int b,int c,t_atom **atom);
+int destroy_lattice(t_atom *atom);
+int thermal_init(t_moldyn *moldyn,t_random *random,int count);
+int scale_velocity(t_moldyn *moldyn,int count);
+double get_e_kin(t_atom *atom,int count);
+double get_e_pot(t_moldyn *moldyn);
+double get_total_energy(t_moldyn *moldyn);
+t_3dvec get_total_p(t_atom *atom,int count);
+
+double estimate_time_step(t_moldyn *moldyn,double nn_dist,double t);
+
+int moldyn_integrate(t_moldyn *moldyn);
+int velocity_verlet(t_moldyn *moldyn);
+
+double potential_harmonic_oscillator(t_moldyn *moldyn);
+int force_harmonic_oscillator(t_moldyn *moldyn);
+double potential_lennard_jones(t_moldyn *moldyn);
+int force_lennard_jones(t_moldyn *moldyn);