1) Adicionei meu codigo de leitor de dados: leitor.c e leitor.h.

Este leitor permite ler os parametros atraves de um arquivo de entrada e tambem pela linha de comando.
O programa pode ser chamado das seguintes formas

./a.out input.txt
./a.out input.txt parA valA parB val B ...
./a.out parA valA parB valB ...

Parametros nao encontrados no arquivo informado ou na lista de parametros na linha de comandos
receberao valores padrao definidos no codigo.

2) Limpei o codigo, mudei nome de algumas variaveis.

3) Corrigi termo de Coulomb na interacao de contato. Nao eh necessario introduzir constante dieletrica.
A intensidade de Hartree eh dada pelo termo eta e eh arbitrario por enquanto.
Incluindo arquivos

input.txt -- exemplo de arquivo de parametros
leitor.c e .h -- implementacao do leitor de dados