UNIVERSIDAD NACIONAL DE SAN AGUSTIN
FACULTAD DE INGENIERÍA DE PRODUCCIÓN Y SERVICIOS
ESCUELA PROFESIONAL DE INGENIERÍA DE SISTEMA
Formato: Guía de Práctica de Laboratorio / Talleres / Centros de Simulación
Aprobación: 2022/03/01 Código: GUIA-PRLMN-001 Página: 1
INFORME DE LABORATORIO 11
(formato estudiante)
INFORMACIÓN BÁSICA
ASIGNATURA: Métodos Numéricos
TITULO DE LA
INTEGRACIÓN NUMÉRICA
PRÁCTICA:
NÚMERO DE NRO.
11 AÑO LECTIVO: 2024B 6
PRÁCTICA: SEMESTRE:
FECHA INICIO 25/11/2024 HORA FIN 30/09/2024
INTEGRANTE (s)
-Callo Ccagiavilca, Kevin Joel
-Hanccoccallo Achircana, Frank Duks
NOTA (0-20)
-Hidalgo Chinchay, Paulo Andre
-Chancuaña Alvis, Klismann
-Perez Huamani, Jeremy Joshua
DOCENTE(s):
Mg. José Rafael Delgado Bastidas
Mg. María Vilma Escobar Castillo
RESULTADOS
I. EJERCICIOS PROPUESTAS:
UNIVERSIDAD NACIONAL DE SAN AGUSTIN
FACULTAD DE INGENIERÍA DE PRODUCCIÓN Y SERVICIOS
ESCUELA PROFESIONAL DE INGENIERÍA DE SISTEMA
Formato: Guía de Práctica de Laboratorio / Talleres / Centros de Simulación
Aprobación: 2022/03/01 Código: GUIA-PRLMN-001 Página: 2
Link:
https://la.mathworks.com/matlabcentral/fileexchange/111490-numerical-methods-with-applicatio
ns
UNIVERSIDAD NACIONAL DE SAN AGUSTIN
FACULTAD DE INGENIERÍA DE PRODUCCIÓN Y SERVICIOS
ESCUELA PROFESIONAL DE INGENIERÍA DE SISTEMA
Formato: Guía de Práctica de Laboratorio / Talleres / Centros de Simulación
Aprobación: 2022/03/01 Código: GUIA-PRLMN-001 Página: 3
Método de Euler
Este código define una función simbólica f(x) como x^3 - 2x^2 + 1 y establece los parámetros para
aplicar el método de Euler: el intervalo [a, b] = [2, 30], la condición inicial y(a)=1, el número de
intervalos (7) y un tiempo de pausa entre pasos (0.5 segundos). Luego, utiliza la función
ShowEulersMethod para visualizar el método de Euler, pasando la función anónima derivada de f(x) y
los parámetros definidos.
Exercise 1:
UNIVERSIDAD NACIONAL DE SAN AGUSTIN
FACULTAD DE INGENIERÍA DE PRODUCCIÓN Y SERVICIOS
ESCUELA PROFESIONAL DE INGENIERÍA DE SISTEMA
Formato: Guía de Práctica de Laboratorio / Talleres / Centros de Simulación
Aprobación: 2022/03/01 Código: GUIA-PRLMN-001 Página: 4
CASO 1.
En este primer caso, aproximamos la solución a la ecuación diferencial dy/dx = 3x^2 - 1 con la
condición inicial y(0)=1 usando el método de Euler. Probamos diferentes valores para el tamaño de
paso hhh para encontrar aquel que genere un error menor a 0.01 al comparar con la solución
exacta y(x) = x^3 - x + 1.
UNIVERSIDAD NACIONAL DE SAN AGUSTIN
FACULTAD DE INGENIERÍA DE PRODUCCIÓN Y SERVICIOS
ESCUELA PROFESIONAL DE INGENIERÍA DE SISTEMA
Formato: Guía de Práctica de Laboratorio / Talleres / Centros de Simulación
Aprobación: 2022/03/01 Código: GUIA-PRLMN-001 Página: 5
CASO 2.
En este caso, resolvemos la ecuación diferencial dx/dy=sin(x) con la condición inicial y(0)=0, usando
el método de Euler para aproximar valores de y en el intervalo [0,π]. Este ejemplo ilustra cómo el
método se ajusta a una función oscilatoria como el seno.
UNIVERSIDAD NACIONAL DE SAN AGUSTIN
FACULTAD DE INGENIERÍA DE PRODUCCIÓN Y SERVICIOS
ESCUELA PROFESIONAL DE INGENIERÍA DE SISTEMA
Formato: Guía de Práctica de Laboratorio / Talleres / Centros de Simulación
Aprobación: 2022/03/01 Código: GUIA-PRLMN-001 Página: 6
CASO 3.
Aquí, usamos el método de Euler para aproximar la solución a dy/dx = e^x con la condición inicial
y(0)=1. Este caso demuestra cómo el método maneja funciones exponenciales, que crecen
rápidamente, verificando la precisión al comparar con la solución exacta y(x) = e^x.
UNIVERSIDAD NACIONAL DE SAN AGUSTIN
FACULTAD DE INGENIERÍA DE PRODUCCIÓN Y SERVICIOS
ESCUELA PROFESIONAL DE INGENIERÍA DE SISTEMA
Formato: Guía de Práctica de Laboratorio / Talleres / Centros de Simulación
Aprobación: 2022/03/01 Código: GUIA-PRLMN-001 Página: 7
CASO 4.
En este caso, aplicamos el método de Euler para resolver dy/dx=(1/(x^2+1)) en el intervalo [0,2] con
la condición inicial y(0)=0. Este ejemplo muestra cómo el método se adapta a una función que
decrece suavemente, comprobando la aproximación con la solución exacta y(x)=arctan(x).
UNIVERSIDAD NACIONAL DE SAN AGUSTIN
FACULTAD DE INGENIERÍA DE PRODUCCIÓN Y SERVICIOS
ESCUELA PROFESIONAL DE INGENIERÍA DE SISTEMA
Formato: Guía de Práctica de Laboratorio / Talleres / Centros de Simulación
Aprobación: 2022/03/01 Código: GUIA-PRLMN-001 Página: 8
Análisis
1. Errores y hhh: Un paso más pequeño (hhh) reduce el error, pero aumenta el costo
computacional.
2. Funciones diferentes: La precisión y el comportamiento del método dependen de la función
derivada.
3. Condiciones iniciales: Cambiar x0x_0x0o y0y_0y0cambia significativamente la solución
aproximada, ya que el método depende de estos valores para iniciar el cálculo.
Usando Series de Taylor para estimar el error
La estimación del error usando series de Taylor se basa en expandir la solución exacta de una
ecuación diferencial y(x) alrededor de un punto inicial X0. En el método de Euler, la solución
aproximada Yaproxutiliza solo los primeros términos de esta expansión, ignorando los términos de
orden superior. El error de un solo paso se deriva de la diferencia entre la solución exacta y la
UNIVERSIDAD NACIONAL DE SAN AGUSTIN
FACULTAD DE INGENIERÍA DE PRODUCCIÓN Y SERVICIOS
ESCUELA PROFESIONAL DE INGENIERÍA DE SISTEMA
Formato: Guía de Práctica de Laboratorio / Talleres / Centros de Simulación
Aprobación: 2022/03/01 Código: GUIA-PRLMN-001 Página: 9
aproximada, lo que resulta ser proporcional a h^2, donde h es el tamaño del paso. Al acumular este
error en varios pasos para alcanzar un punto final, el error total es aproximadamente proporcional a
h, lo que refleja que el método de Euler es de orden 1 en precisión.
Si comenzamos con una función y la integramos, obtenemos:
Al utilizar el desarrollo en Serie de Taylor para la función y(x), podemos escribirla como:
Método de Euler y error de un solo paso
El método de Euler aproxima la solución y(x) como:
donde h es el tamaño del paso.
Sustituyendo desde la serie de Taylor:
Testing Euler's Method
Ejemplos:
UNIVERSIDAD NACIONAL DE SAN AGUSTIN
FACULTAD DE INGENIERÍA DE PRODUCCIÓN Y SERVICIOS
ESCUELA PROFESIONAL DE INGENIERÍA DE SISTEMA
Formato: Guía de Práctica de Laboratorio / Talleres / Centros de Simulación
Aprobación: 2022/03/01 Código: GUIA-PRLMN-001 Página: 10
Resultados:
UNIVERSIDAD NACIONAL DE SAN AGUSTIN
FACULTAD DE INGENIERÍA DE PRODUCCIÓN Y SERVICIOS
ESCUELA PROFESIONAL DE INGENIERÍA DE SISTEMA
Formato: Guía de Práctica de Laboratorio / Talleres / Centros de Simulación
Aprobación: 2022/03/01 Código: GUIA-PRLMN-001 Página: 11
Example 3.
UNIVERSIDAD NACIONAL DE SAN AGUSTIN
FACULTAD DE INGENIERÍA DE PRODUCCIÓN Y SERVICIOS
ESCUELA PROFESIONAL DE INGENIERÍA DE SISTEMA
Formato: Guía de Práctica de Laboratorio / Talleres / Centros de Simulación
Aprobación: 2022/03/01 Código: GUIA-PRLMN-001 Página: 12
Try different values.
UNIVERSIDAD NACIONAL DE SAN AGUSTIN
FACULTAD DE INGENIERÍA DE PRODUCCIÓN Y SERVICIOS
ESCUELA PROFESIONAL DE INGENIERÍA DE SISTEMA
Formato: Guía de Práctica de Laboratorio / Talleres / Centros de Simulación
Aprobación: 2022/03/01 Código: GUIA-PRLMN-001 Página: 13
UNIVERSIDAD NACIONAL DE SAN AGUSTIN
FACULTAD DE INGENIERÍA DE PRODUCCIÓN Y SERVICIOS
ESCUELA PROFESIONAL DE INGENIERÍA DE SISTEMA
Formato: Guía de Práctica de Laboratorio / Talleres / Centros de Simulación
Aprobación: 2022/03/01 Código: GUIA-PRLMN-001 Página: 14
Ejecución:
Simpson´s rule
Visualizing Simpson´s rule
UNIVERSIDAD NACIONAL DE SAN AGUSTIN
FACULTAD DE INGENIERÍA DE PRODUCCIÓN Y SERVICIOS
ESCUELA PROFESIONAL DE INGENIERÍA DE SISTEMA
Formato: Guía de Práctica de Laboratorio / Talleres / Centros de Simulación
Aprobación: 2022/03/01 Código: GUIA-PRLMN-001 Página: 15
UNIVERSIDAD NACIONAL DE SAN AGUSTIN
FACULTAD DE INGENIERÍA DE PRODUCCIÓN Y SERVICIOS
ESCUELA PROFESIONAL DE INGENIERÍA DE SISTEMA
Formato: Guía de Práctica de Laboratorio / Talleres / Centros de Simulación
Aprobación: 2022/03/01 Código: GUIA-PRLMN-001 Página: 16
UNIVERSIDAD NACIONAL DE SAN AGUSTIN
FACULTAD DE INGENIERÍA DE PRODUCCIÓN Y SERVICIOS
ESCUELA PROFESIONAL DE INGENIERÍA DE SISTEMA
Formato: Guía de Práctica de Laboratorio / Talleres / Centros de Simulación
Aprobación: 2022/03/01 Código: GUIA-PRLMN-001 Página: 17
Ejecución:
4
UNIVERSIDAD NACIONAL DE SAN AGUSTIN
FACULTAD DE INGENIERÍA DE PRODUCCIÓN Y SERVICIOS
ESCUELA PROFESIONAL DE INGENIERÍA DE SISTEMA
Formato: Guía de Práctica de Laboratorio / Talleres / Centros de Simulación
Aprobación: 2022/03/01 Código: GUIA-PRLMN-001 Página: 18
Se construyó una tabla para comprobar el orden del error de una integral sobre límites finitos con
diferentes números de subintervalos.
UNIVERSIDAD NACIONAL DE SAN AGUSTIN
FACULTAD DE INGENIERÍA DE PRODUCCIÓN Y SERVICIOS
ESCUELA PROFESIONAL DE INGENIERÍA DE SISTEMA
Formato: Guía de Práctica de Laboratorio / Talleres / Centros de Simulación
Aprobación: 2022/03/01 Código: GUIA-PRLMN-001 Página: 19
1. La regla de Simpson con 𝑛=1 calcula (b), (c) y (d) exactamente, pero no (a).
2. Para la regla de Simpson, el error es 0 cuando la función que se está integrando es un polinomio
de grado ≤ 3. En este caso, el tamaño del paso ℎ no importa porque la aproximación es exacta
para cualquier ℎ. Sin embargo, para funciones de grado superior a 3, un valor ℎ más pequeño
reduce el error.
¿Qué sucede si se limita el tamaño del paso?
Mayor eficiencia: Menos evaluaciones de función en regiones suaves.
UNIVERSIDAD NACIONAL DE SAN AGUSTIN
FACULTAD DE INGENIERÍA DE PRODUCCIÓN Y SERVICIOS
ESCUELA PROFESIONAL DE INGENIERÍA DE SISTEMA
Formato: Guía de Práctica de Laboratorio / Talleres / Centros de Simulación
Aprobación: 2022/03/01 Código: GUIA-PRLMN-001 Página: 20
Mayor precisión: Mejor resolución en regiones donde la solución cambia rápidamente.
Control del error: Permite mantener el error dentro de un umbral tolerable.
5
Modificando la regla de Simpson
Actualizar la regla de Simpson requiere un poco más de esfuerzo porque debemos pensar detenidamente sobre
lo que significa encajar una cuadrática a través de los puntos que conocemos. Habrá diferentes casos para la
regla de Simpson dependiendo de si tenemos un número par o impar de subintervalos.
UNIVERSIDAD NACIONAL DE SAN AGUSTIN
FACULTAD DE INGENIERÍA DE PRODUCCIÓN Y SERVICIOS
ESCUELA PROFESIONAL DE INGENIERÍA DE SISTEMA
Formato: Guía de Práctica de Laboratorio / Talleres / Centros de Simulación
Aprobación: 2022/03/01 Código: GUIA-PRLMN-001 Página: 21
Número impar de subintervalos
Si tenemos un número impar de subintervalos, no podemos emparejarlos todos para crear curvas
cuadráticas apropiadas a través de cada triplete. Una forma de abordar este caso es eliminar el
subintervalo final para reducirnos al caso de subintervalos pares y luego usar los dos últimos
subintervalos para crear una curva cuadrática aproximada para usar en la integración numérica del
subintervalo final.
UNIVERSIDAD NACIONAL DE SAN AGUSTIN
FACULTAD DE INGENIERÍA DE PRODUCCIÓN Y SERVICIOS
ESCUELA PROFESIONAL DE INGENIERÍA DE SISTEMA
Formato: Guía de Práctica de Laboratorio / Talleres / Centros de Simulación
Aprobación: 2022/03/01 Código: GUIA-PRLMN-001 Página: 22
No se puede actualizar la cuadratura gaussiana
La cuadratura gaussiana se basa en conocer los valores en puntos particulares. La cuadratura
gaussiana no se puede usar si no se conocen esos valores de integrando.
UNIVERSIDAD NACIONAL DE SAN AGUSTIN
FACULTAD DE INGENIERÍA DE PRODUCCIÓN Y SERVICIOS
ESCUELA PROFESIONAL DE INGENIERÍA DE SISTEMA
Formato: Guía de Práctica de Laboratorio / Talleres / Centros de Simulación
Aprobación: 2022/03/01 Código: GUIA-PRLMN-001 Página: 23
Aplicación de la interpolación
Un método que se puede usar para regularizar datos irregulares para su uso en métodos numéricos
estándar es interpolar los valores necesarios en los puntos espaciados uniformemente deseados.
Test de funciones
1. Se crean carpetas con los nombres de las funciones
UNIVERSIDAD NACIONAL DE SAN AGUSTIN
FACULTAD DE INGENIERÍA DE PRODUCCIÓN Y SERVICIOS
ESCUELA PROFESIONAL DE INGENIERÍA DE SISTEMA
Formato: Guía de Práctica de Laboratorio / Talleres / Centros de Simulación
Aprobación: 2022/03/01 Código: GUIA-PRLMN-001 Página: 24
UNIVERSIDAD NACIONAL DE SAN AGUSTIN
FACULTAD DE INGENIERÍA DE PRODUCCIÓN Y SERVICIOS
ESCUELA PROFESIONAL DE INGENIERÍA DE SISTEMA
Formato: Guía de Práctica de Laboratorio / Talleres / Centros de Simulación
Aprobación: 2022/03/01 Código: GUIA-PRLMN-001 Página: 25
UNIVERSIDAD NACIONAL DE SAN AGUSTIN
FACULTAD DE INGENIERÍA DE PRODUCCIÓN Y SERVICIOS
ESCUELA PROFESIONAL DE INGENIERÍA DE SISTEMA
Formato: Guía de Práctica de Laboratorio / Talleres / Centros de Simulación
Aprobación: 2022/03/01 Código: GUIA-PRLMN-001 Página: 26
UNIVERSIDAD NACIONAL DE SAN AGUSTIN
FACULTAD DE INGENIERÍA DE PRODUCCIÓN Y SERVICIOS
ESCUELA PROFESIONAL DE INGENIERÍA DE SISTEMA
Formato: Guía de Práctica de Laboratorio / Talleres / Centros de Simulación
Aprobación: 2022/03/01 Código: GUIA-PRLMN-001 Página: 27
UNIVERSIDAD NACIONAL DE SAN AGUSTIN
FACULTAD DE INGENIERÍA DE PRODUCCIÓN Y SERVICIOS
ESCUELA PROFESIONAL DE INGENIERÍA DE SISTEMA
Formato: Guía de Práctica de Laboratorio / Talleres / Centros de Simulación
Aprobación: 2022/03/01 Código: GUIA-PRLMN-001 Página: 28
CÓDIGOS DE TESTEO
UNIVERSIDAD NACIONAL DE SAN AGUSTIN
FACULTAD DE INGENIERÍA DE PRODUCCIÓN Y SERVICIOS
ESCUELA PROFESIONAL DE INGENIERÍA DE SISTEMA
Formato: Guía de Práctica de Laboratorio / Talleres / Centros de Simulación
Aprobación: 2022/03/01 Código: GUIA-PRLMN-001 Página: 29
UNIVERSIDAD NACIONAL DE SAN AGUSTIN
FACULTAD DE INGENIERÍA DE PRODUCCIÓN Y SERVICIOS
ESCUELA PROFESIONAL DE INGENIERÍA DE SISTEMA
Formato: Guía de Práctica de Laboratorio / Talleres / Centros de Simulación
Aprobación: 2022/03/01 Código: GUIA-PRLMN-001 Página: 30
UNIVERSIDAD NACIONAL DE SAN AGUSTIN
FACULTAD DE INGENIERÍA DE PRODUCCIÓN Y SERVICIOS
ESCUELA PROFESIONAL DE INGENIERÍA DE SISTEMA
Formato: Guía de Práctica de Laboratorio / Talleres / Centros de Simulación
Aprobación: 2022/03/01 Código: GUIA-PRLMN-001 Página: 31
UNIVERSIDAD NACIONAL DE SAN AGUSTIN
FACULTAD DE INGENIERÍA DE PRODUCCIÓN Y SERVICIOS
ESCUELA PROFESIONAL DE INGENIERÍA DE SISTEMA
Formato: Guía de Práctica de Laboratorio / Talleres / Centros de Simulación
Aprobación: 2022/03/01 Código: GUIA-PRLMN-001 Página: 32
CONCLUSIONES
● vimos como la regla de Simpson y la interpolación de Lagrange son aplicadas en matlab. Básicamente,
estuvimos explorando cómo aproximar valores y resolver integrales de manera más eficiente usando
herramientas matemáticas y gráficas en MATLAB. Aunque son métodos útiles, sólo representan una
pequeña sección del panorama completo del análisis numérico y su aplicación práctica. Lo interesante es
cómo podemos combinar teoría y código para visualizar resultados y entender mejor los errores o
limitaciones que pueden surgir.
UNIVERSIDAD NACIONAL DE SAN AGUSTIN
FACULTAD DE INGENIERÍA DE PRODUCCIÓN Y SERVICIOS
ESCUELA PROFESIONAL DE INGENIERÍA DE SISTEMA
Formato: Guía de Práctica de Laboratorio / Talleres / Centros de Simulación
Aprobación: 2022/03/01 Código: GUIA-PRLMN-001 Página: 33
METODOLOGÍA DE TRABAJO
REFERENCIAS Y BIBLIOGRAFÍA
[1] R. P. C. Steven C. Chapra, Métodos numéricos para ingenieros, Mexico: Mc Graw Hill Education, 2015.
[2] J. G. Hernández, Métodos numéricos con applicaciones en Matlab, Quito: Universitaria Abya - Yala, 2013.
[3] K. D. F. John H. Mathews, Métodos Numéricos con MATLAB, Madrid: PEARSON EDUCATION S.A., 2000.
