# Javier M. Hernández: Docencia

## Física Computacional

### Lecturas:

• L1: Cap. 6, del Survey Computational Physics, Landau
• L2: Cap. 1 y A.1-A.3, del Primer Scientific Programming
• L3: Cap. 7, sección 7.4 del Primer Scientific Programming
• L4: Cap. 9, del Survey Computational Physics, Landau
• L5: Álgebra Lineal
• L6: Ven alguna estructura artificial en esta imagen que tomó la NASA en Marte?
• notas LaTeX

### Programas:

1. Manejo de listas en python
2. Ejemplos de estructuras de control en python
3. Manejo de listas (II) en python
4. Método Babilónico para hallar raíces cuadradas.
5. Método Babilónico para hallar raíces cuadradas. (II)
6. Método Forward Euler para ecs. diferenciales
7. Archivo base latex
8. Funciones
9. Pi con MC. 10. Para calcular integrales de varias variables, , (cuyo resultado exacto es 155/6) podemos usar inte-mc-2.py (que usa como auxiliar a mcint.py, deben ponerlo en el mismo directorio). ejemplo de la salida.
11. Histogramas 1

### Tareas:

• A1: Leer la página del curso. Enviarme un email a jmhlopez+FisC en gmail.com, con encabezado: FisC-su nombre. En el cuerpo del email confirmen que han leido y aceptan las reglas escritas en la página. Anexen su matrícula.
• T1: coseno2.py
• T2: caida-libre.py
• T3: decaimiento.py (resolver N(t+dt)=(1-k dt) N(t), poner N(0)=10^3, usar k=1,100)
• T4: Experimento numéricos con la fórmula trapezoidal de integración
(Implementar el Linsting 6.1 (Landau), trapezoidal.py (pag. 587-588 del Primer) y programa trapezoid, pag. 203 del Numerical Method in Engineering).
2a. parte: integrar la función $$f(t)= t+1, \ si \ -2\le t \le -1; 2, \ si \ -1 \le t \le 0; -2t+4, \ si \ 0 \le t \le 1;$$ en el intervalo $[-2,1]$.
• T5: Implementar Listing 6.2 (Landau). Comparar resultados de ∫ e^{-x} dx, 0≤x≤10 usando trap(10), trap(1000) y Gauss-Legendre(4 puntos).
• T6: Del temario del curso, en la página inicial, definir puntos 1.1 y 1.2, en un archivo latex.
• T7: Implementar listado 9.2 del Landau. Resolver con él los puntos 3 y 5 de la pag. 191. Resolver las preguntas de la sección 9.6 (non linear oscillations) del Landau.
Tendrán problemas para que corra el listado 9.2 debido a que utiliza las bibliotecas de visual para las gráficas (ejemplo). Sin embargo podemos eliminar dicha dependencia (versión m de rk45.py) y obtener las mismas gráficas en gnuplot. Noten que salve la información a un archivo para manipularlo en gnuplot. Uds. pueden mandar los datos directamente a un archivo desde rk45-m.py. EL código para hacer plot conjuntos en gnuplot es:
set multiplot;
set size 1,0.5;
set origin 0.0,0.5; plot sin(x);
set origin 0.0,0.0; plot cos(x);

• T7. Con el programa pi-mc.py, para N=10000, correrlo 100 veces. Obtener promedio y desviación estándar. Cambiar N a 1 millón. Obtener promedio y desviación estándar. Revisen el código de mcint.py para el cálculo de la media y la desviación estándar que necesitan. The principle of science, the definition almost, is the following: The test of all knowledge is experiment. Experiment is the sole judge of scientific truth. But what is the source of knowledge? Where do the laws that are to be tested come from? Experiment, itself, helps to produce these laws, in the sense that it gives us hints. But also needed is imagination to create from these hints the great generalizations to guess at the wonderful, simple, but very strange patterns beneath them all. and then to experiment to check again whether we made the right guess. Richard Feynman Última revisión: 09/05/2018, javierh en {fcfm.buap.mx} , Puebla, Pue. México