En este trabajo se implemento un modelo para el estudio de estructura de bandas de energia. Este modelo se baso en el metodo de diferencias finitas, usando un arreglo 3D. Inicialmente, la celda unitaria mallada emplea un paso h, con el fin de dividir la celda unitaria en partes mas pequenas. Cada punto de la malla se represento por medio de una funcion de onda. Se obtuvieron las derivadas aproximadas de cada punto en sus direcciones (x, y y z), empleando la diferencia entre los puntos vecinos a lo largo de los ejes. Las funciones de onda descritas por las funciones de Bloch, el laplaciano y las derivadas son incluidas en la ecuacion de Schrodinger. Luego, se resuelve dicha ecuacion con el fin de obtener los valores propios (las energias), graficadas en funcion del vector de onda, con el fin de construir las estructuras de bandas. Para estudiar su desempeno, este modelo fue aplicado a diferentes celdas cubicas, como cubica simple, cubica centrada en el cuerpo (BCC) y cubica centrada en las caras (FCC), con dimensiones unitarias