UNIVERSIDAD DE GUAYAQUIL

FACULTAD DE INGENIERIA INDUSTRIAL

MÉTODOS NUMÉRICOS II
TEMA:

APROXIMACION A UNA ECUACION DIFERENCIAL DE SEGUNDO ORDEN DE UN

CIRCUITO RLC POR DIFERENCIAS DIVIDIDAS FINITAS CENTRADAS.

ESTUDIANTE:
TITO MICHAEL LOOR GAVILANES

DOCENTE:
Ing. Edin Alex Garcés Coka

GUAYAQUIL- ECUADOR
Contenido
INTRODUCCIÓN
OBJETIVOS
OBJETIVO GENERAL
OBJETIVOS ESPECÍFICOS
PLANTEAMIENTO DEL PROBLEMA
JUSTIFICACION
MARCO TEÓRICO
El Método de Diferencias Finitas
Fórmula de Euler
CIRCUITO RLC
LEY DE KIRCHHOFF.
MATRIZ DE VANDERMONDE.
METODOLOGIA DE LA INVESTIGACION
RESULTADOS.
CONCLUSIÓN.
RECOMENDACIONES.
BIBLIOGRAFÍA.
ANEXOS.
INDICE DE TABLAS
Tabla 1 Tabla de datos
Tabla 2 Diferencia del i real y su aproximado
 INDICE DE FIGURAS.
Figura 1 diferencias centradas. Ing Jorge Mario Peña Consuegra.
Figura 2. Respuesta Natural de Circuito RLC. Elaboración propia.
Figura 3 Condiciones del Circuito Justo Antes que se Cierre el Interruptor. Elaboración propia.
Figura 4 Condiciones del Circuito Justo Después que se Cierra el Interruptor. Elaboración propia.
Figura 5 Ejercicio de la Respuesta Natural del Circuito RLC. Elaboración propia.
Figura 6 corriente 0; voltaje inicial 10v. Elaboración propia
Figura 7 t=0, i=0, vc=10v. Elaboración Propia.
Figura 8 Onda del Seno. Elaboración propia.
Figura 9. Comparación de la curva real y la aproximación. Elaboración propia
 INDICE DE ANEXOS
Anexo 1. Script Para Diferencias Dividivas Centradass
Anexo 2 Datos evaluados en Matlab.
Anexo 3. Polinomio Aproximado
 INTRODUCCIÓN

 En el siguiente trabajo se describe la construcción de un polinomio para aproximarlo a la
función verdadera y así analizar el comportamiento de la intensidad del circuito RLC en serie
con respecto al tiempo.
 Para lograr este propósito aplicamos el método de diferencias divididas centradas que nos
permitirá resolver analíticamente y de manera más sencilla la ecuación diferencial de
segundo orden ya que de la forma por el método exacto de análisis es muy difícil.
 En la actualidad cualquier fenómeno en la naturaleza puede ser descrito con las ayudas de la
física, en términos de ecuaciones algebraicas o de ecuaciones diferenciales.
 Los ingenieros y científicos estudian estos fenómenos de dos maneras:
 Formulación del proceso físico
 Análisis numérico del modelo trabajo.
 En este caso desarrollar el modelo matemático de un proceso es logrado a través de las
suposiciones de cómo trabaja el proceso. En una solución numérica utilizamos un método
numérico y una computadora para evaluar este modelo.
 OBJETIVOS

 OBJETIVO GENERAL
 Construir un polinomio que sustituya a la función original del circuito RLC y determinar la
corriente del mismo.

 OBJETIVOS ESPECÍFICOS

 Resolver un ejercicio de ecuación diferencial homogénea de segundo orden

analíticamente por el método exponencial.
 Aplicar la fórmula de diferencias divididas centradas a la fórmula original del circuito RLC.
 Permitir estudiar el comportamiento de las gráficas utilizando la herramienta de
programación MATLAB.
 PLANTEAMIENTO DEL PROBLEMA
 Vamos a resolver un ejercicio de circuito RLC por medio de una ecuación diferencial de
segundo orden y luego por el método de diferencias finitas centradas para calcular su
corriente (i) y luego comparar las curvas que nos dan ambos métodos por medio del software
Matlab.
El capacitor tiene un voltaje inicial de 10v.
 No hay corriente que fluya en el inductor al momento que se cierra el interruptor.

 JUSTIFICACION
 El circuito RLC ha sido base de nuestro estudio el cual tiene por objetivo encontrar un
polinomio que se aproxime a la ecuación original para encontrar una solución eficiente al
problema expresado matemáticamente.
 Para esto es necesario cumplir con ciertos objetivos específicos:
 Resolver la ecuación homogénea de segundo orden por el método exponencial para luego
compáralo con el método de diferencias finitas.
 Utilizar Matlab para estudiar el comportamiento de las gráficas.
 Esta investigación es importante porque nos ayuda a resolver de manera más sencilla la
ecuación diferencial.
 MARCO TEÓRICO
 Un problema muy común en el campo de las ciencias, la ingeniería o algún fenómeno
de la naturaleza puede ser representado por una ecuación. En esta investigación se
expondrá un método numérico para encontrar la solución a esa ecuación, pero antes
de eso analizaremos ciertos temas de interés que nos ayudaran a entender el proceso
que llevaremos a cabo.
 ¿Qué es un método numérico?
 Un método numérico es un procedimiento durante el cual, se obtiene casis siempre de
manera aproximada, la solución de ciertos problemas realizando cálculos puramente
aritméticos y lógicos. Un tal procedimiento consiste de una lista finita de instrucciones
precisas que especifican una secuencia de operaciones algebraicas y
lógicas(algoritmos), que producen o bien una aproximación de la solución del problema
(solución numérica) o bien un mensaje. La eficiencia en el cálculo de dicha
aproximación depende, en parte, de la facilidad de implementación del algoritmo y de
las características especiales y limitaciones de los instrumentos de cálculo (los
computadores).
 METODOLOGIA DE LA INVESTIGACION
 MÉTODO
 Método teórico
 Modelación matemática y computacional.
 VARIABLES
 Independiente: resistencia, inductancia, capacitancia, voltaje y corriente.
 Dependiente: amplitud de oscilación de la corriente en el tiempo.

 Modelamos el circuito RCL con una ecuación diferencial de segundo orden, con la
corriente con la variable independiente.
 Figura 3 Condiciones del Circuito Justo Antes que
Figura 4 Condiciones del Circuito Justo Después
se Cierre el Interruptor. Elaboración propia.
que se Cierra el Interruptor. Elaboración propia.
 Escribimos las ecuaciones para cada elemento individual:

𝑑𝑖
 𝑉𝑅 = 𝐼 ∗ 𝑅 𝑑𝑉𝑅 = 𝑑𝑡 ∗ 𝑅
𝑑𝑖 𝑑2 𝑖
 𝑉𝐿 = 𝐿 ∗ 𝑑𝑉𝐿 = 𝐿 ∗
𝑑𝑡 𝑑𝑡 2
1 1
 𝑉𝑐 = ‫𝑡𝑑𝑖 𝑐 ׬‬ 𝑑𝑉𝐶 = 𝑐 𝑖

 𝑉 = 𝑘𝑡𝑒 𝑑𝑉 = 0
 Escribimos la ley de Kirchhoff del voltaje y sumamos los voltajes conforme recorramos la
malla en sentido de las manecillas del reloj:
 𝑣𝐿 − 𝑣𝑅 − 𝑣𝐶 = 0
 Reemplazamos los términos de 𝑣 con los correspondientes 𝑖:
𝑑𝑖 1
 𝐿 𝑑𝑡 + 𝑅𝑖 + 𝑐 ‫ = 𝑡𝑑𝑖 ׬‬0
 Quitamos toda la integral derivando toda la ecuación:
𝑑 𝑑𝑖 1
 𝑑𝑡
𝐿 𝑑𝑡 + 𝑅𝑖 + 𝑐 ‫ = 𝑡𝑑𝑖 ׬‬0

 Esto se llama ecuación de segundo orden homogénea, ahora nos disponemos a resolver nuestra ecuación:
𝑑𝑖 2 𝑑𝑖 1
 𝐿 𝑑𝑡 2 + 𝑅 𝑑𝑡 + 𝑐 𝑖 = 0

 Aplicamos la forma exponencial para resolver la ecuación:

 𝑖 𝑡 = 𝑘𝑒 𝑠𝑡
𝑑2 𝑑 1
 𝐿 𝑑𝑡 2 𝑘𝑒 𝑠𝑡 + 𝑅 𝑑𝑡 𝑘𝑒 𝑠𝑡 + 𝑐 𝑘𝑒 𝑠𝑡 = 0

 Ahora trabajamos con las derivadas:

 Derivamos el primer término dos veces:
𝑑2
 𝐿 𝑑𝑡 2 𝑘𝑒 𝑠𝑡 = 𝑠 2 𝐿𝐾𝑒 𝑠𝑡

 Derivamos el segundo término:

𝑑
 𝑅 𝑘𝑒 𝑠𝑡 = 𝑠𝑅𝐾𝑒 𝑠𝑡
𝑑𝑡

 Sustituimos en la ecuación diferencial:

1
 𝑠 2 𝐿𝐾𝑒 𝑠𝑡 + 𝑠𝑅𝐾𝑒 𝑠𝑡 + 𝑐 𝑘𝑒 𝑠𝑡 = 0

 Factorizamos el termino común:

1
 𝐾𝑒 𝑠𝑡 𝑠 2 𝐿 + 𝑠𝑅 + =0
𝑐

 Hacemos que, la ecuación se haga verdadera si igualamos s=0

1
 𝑠 2 𝐿 + 𝑠𝑅 + = 0
𝑐
 A esto llamamos ecuación característica del circuito RLC.
 ENCONTRAR LAS RAICES DE LA ECUACIÓN CARACTERISTICA.
 Aplicamos la fórmula cuadrática para resolver la ecuación:
𝑏2 −4𝑎𝑐
 𝑥 = −𝑏 ± 2𝑎

𝑅2 −4𝐿/𝐶
 𝑠 = −𝑅 ± 2𝐿


 La fórmula cuadrática nos dio dos soluciones para s, las cuales llamaremos s1 y s2 las
incluimos en la solución propuesta y tenemos una combinación lineal de dos términos
exponenciales separados con cuatro términos ajustables.
 𝑖 = 𝐾1 𝑒 𝑠1𝑡 + 𝐾2 𝑒 𝑠2𝑡
 S1 y s2 son frecuencias naturales.
 K1 y K2 son términos de amplitud.
EJERCICIO

Figura 5 Ejercicio de la Respuesta Natural

del Circuito RLC. Elaboración propia.
 La ecuación diferencial del circuito es:
𝑑𝑖 2 𝑑𝑖 1
 𝐿 𝑑𝑡 2 + 𝑅 𝑑𝑡 + 𝑐 𝑖 = 0

 Con los valores reales se vuelve:

𝑑𝑖 2 𝑑𝑖 1
 1 𝑑𝑡 2 + 2 𝑑𝑡 + 5 𝑖 = 0

 Suponemos una solución de la forma:

 𝑖 𝑡 = 𝑘𝑒 𝑠𝑡
 Realizamos el análisis anterior y nos resulta la siguiente ecuación característica:

 𝑠 2 + 2𝑠 + 5 = 0
 Aplicamos la formula cuadrática para resolver la ecuación:

𝑅2 −4𝐿/𝐶
 𝑠 = −𝑅 ± 2𝐿

22 −4∗1∗5
 𝑠 = −2 ± 2∗1

 𝑠 = −1 ± 𝑗2
 Utilizamos la j para el numero imaginario.
 Utilizamos i para la corriente.
 La solución propuesta para (𝑖) está escrita como la superposición de dos términos exponenciales diferentes:

 𝑖 = 𝐾1 𝑒 (−1+𝑗2.2)𝑡 + 𝐾2 𝑒 (−1−𝑗2.2)𝑡
 Separamos la parte real de la imaginaria:

 𝑖 = 𝐾1 𝑒 −1𝑡 𝑒 𝑗2𝑡 + 𝐾2 𝑒 −1𝑡 𝑒 −𝑗2𝑡

 Factorizamos el término común:

 𝑖 = 𝑒 −1𝑡 (𝐾1 𝑒 𝑗2𝑡 + 𝐾2 𝑒 −𝑗2𝑡 )

 FORMULA DE EULER
 Para todo número real x. Aquí, e es la base del logaritmo natural, i es la unidad imaginaria y sin, cos son funciones
trigonométricas.

 𝑒 +𝑗𝑥 = 𝑐𝑜𝑠𝑥 + 𝑗𝑠𝑖𝑛

 𝑒 −𝑗𝑥 = 𝑐𝑜𝑠𝑥 − 𝑗𝑠𝑖𝑛
 Aplicamos a la ecuación:
 𝐾1 𝑒 𝑗2𝑡 + 𝐾2 𝑒 −𝑗2𝑡
 𝑘1 𝑐𝑜𝑠2𝑡 + 𝑗𝑠𝑖𝑛2𝑡 + 𝑘2 (𝑐𝑜𝑠2𝑡 − 𝑗𝑠𝑖𝑛2𝑡)
 Multiplicamos:
 𝑘1 𝑐𝑜𝑠2𝑡 + 𝑗𝑘1 𝑠𝑖𝑛2𝑡 + 𝑘2 𝑐𝑜𝑠2𝑡 − 𝑗𝑘2 𝑠𝑖𝑛2𝑡
 Agrupamos los términos del seno y coseno:
 𝑘1 + 𝑘2 𝑐𝑜𝑠2𝑡 + 𝑗 𝑘1 − 𝑘2 𝑠𝑖𝑛2𝑡
 Reemplazamos en la ecuación donde se encuentra K con otras incógnitas en términos de A y se vera de la
siguiente manera:
 𝐴1=𝐾1+𝐾2
 𝐴2=𝐽(𝐾1−𝐾2)

 𝐴1 𝑐𝑜𝑠2𝑡 + 𝐴2 𝑠𝑖𝑛2𝑡
 Volvemos a la fórmula propuesta:
 𝑖 = 𝑒 −𝑡 (𝐴1 𝑐𝑜𝑠2𝑡 + 𝐴2 𝑠𝑖𝑛2𝑡)
 ENCONTRAR LAS CONDICIONES INICIALES.
 Es lo mismo que el circuito LC. Solo tomamos en cuenta el resistor.
 Sabemos que:
 𝑡 = 0− momento antes de que se cierre el interruptor.
 La corriente es 0 y el voltaje inicial en el capacitor es Vc=10v.
 El voltaje inicial del capacitor esta expresado como 𝑣𝑐 0− = 𝑣0
 Según algunas propiedades de los inductores y capacitores nos indica que pasa cuando
se cierra el interruptor yendo de 𝑡 = 0− a 𝑡 = 0+
 La corriente en el inductor no cambia de manera instantánea de modo que 𝑖 0+ =
𝑖 0− = 0
 El voltaje en el capacitor no cambia de manera instantánea de modo que 𝑣𝑐 0+ =
𝑣𝑐 0− = 𝑣0
 Ahora conocemos la condición inicial 𝑖 0+ = 0 y algo del voltaje pero no conocemos
𝑑𝑖/𝑑𝑡(0+ )
 Para la segunda condición aplicamos la ley de Kirchhoff alrededor de la malla:
 +𝑉𝐿 − 𝑉𝑅 − 𝑉𝐶 = 0
 Como 𝑖 0+ = 0 eso significa que el voltaje a través del resistor tiene que ser 0. También
sabemos que el voltaje a través del capacitor es 𝑣𝑐= 𝑣0 .
 +𝑉𝐿 − 0 − 𝑉0 = 0
 𝑉𝐿 = 𝑉0
 Usamos esto para derivar di/dt al usar la ecuación i-v del inductor.
𝑑𝑖
 𝑉𝐿 (0+ ) = 𝐿 𝑑𝑡 (0+ )
𝑑𝑖
 10 = 1 𝑑𝑡 (0+ )
𝑑𝑖
 0+ = 10𝐴/𝑠
𝑑𝑡
 Si es evaluamos i=0, podemos encontrar una de las constantes A. sustituimos 𝑡 = 0 e 𝑖 = 0
en la solución propuesta:
 0 = 𝑒 −0 (𝐴1 𝑐𝑜𝑠2 ∗ 0 + 𝐴2 𝑠𝑖𝑛2 ∗ 0)
 0 = 1(𝐴1 ∗ 1 + 𝐴2 ∗ 0)
 𝐴1 = 0
 El termino en el coseno se elimina y nos queda la solución de la siguiente manera:
 𝑖 = 𝑒 −𝑡 (𝐴2 𝑠𝑖𝑛2 ∗ 𝑡)
 Ahora vamos a buscar A2 para esto tomamos la derivada de la solución propuesta:
𝑑𝑖 𝑑
 = (𝐴2 𝑠𝑖𝑛2 ∗ 𝑡)
𝑑𝑡 𝑑𝑡

 Utilizamos la regla del producto y la respuesta nos queda:

𝑑𝑖
 = 𝐴2 𝑒 −𝑡 (2𝑐𝑜𝑠2𝑡 − 𝑠𝑖𝑛2𝑡)
𝑑𝑡

 Evaluamos esta expresión en 𝑡 = 0+
 10 = 𝐴2 𝑒 −0 (2𝑐𝑜𝑠2 ∗ 0 − 𝑠𝑖𝑛2 ∗ 0)
 10 = 𝐴2 (2𝐴1 ∗ 1 − 0)
 10 = 𝐴2 2
 𝐴2 = 5
 SOLUCION PARA LA CORRIENTE.
 𝑖 = 5𝑒 −𝑡 𝑠𝑖𝑛2𝑡
 PROPUESTA PARA EL PROBLEMA.
 Tenemos los datos:
𝑑𝑖 𝑑𝑖 2 𝑑𝑖 1
 = 𝐿 2 +𝑅 + 𝑖
𝑑𝑡 𝑑𝑡 𝑑𝑡 𝑐

 Condiciones iniciales:
 𝑖 0 =0
 𝑖 ′ 0 = 10𝐴/𝑠
 Rango= [0;6]
 H=0.5
 A continuación, tenemos la fórmula para aproximación de diferencia dividas finitas centradas para 2
puntos que es la que utilizaremos:
 Para la primera derivada:

𝑓 𝑥0 +ℎ −𝑓(𝑥0 −ℎ)
 𝑓 ′ 𝑥0 =
2ℎ

 Para la segunda derivada:

𝑓 𝑥0 −ℎ −2𝑓 𝑥0 +𝑓(𝑥0 +ℎ)
 𝑓 ′′ 𝑥0 =
ℎ2


 Ajustamos la fórmula a nuestro parámetro que es la corriente (i) y la formula nos queda de la siguiente
manera:

𝑖𝑖+1 −𝑖𝑖−1
 𝑖′ =
2ℎ
𝑖𝑖−1 −2𝑖𝑖 +𝑖𝑖+1
 𝑖 ′′ =
ℎ2


 𝑖 0 =0
 𝑖 ′ 0 = 10𝐴/𝑠

𝑖𝑖+1 −𝑖𝑖−1
 10 =
2ℎ

 𝑖𝑖−1 = 𝑖𝑖+1 − 10 ∗ 2 ∗ 0.5

 𝑖𝑖−1 = 𝑖𝑖+1 − 10
𝑖𝑖−1 −2𝑖𝑖 +𝑖𝑖+1
 𝑖 ′′ =
ℎ2
𝑅 1
 (𝑖𝑖−1 −2𝑖𝑖 + 𝑖𝑖+1 ) + (𝑖𝑖+1 − 𝑖𝑖+1 + 10) + 𝑖
𝐿 𝐿𝐶 𝑖

 𝑖𝑖+1 = 2.5
 Remplazamos las fórmulas en nuestra ecuación diferencial:
𝑑𝑖 2 𝑑𝑖 1
 𝐿 𝑑𝑡 2 + 𝑅 𝑑𝑡 + 𝑐 𝑖 = 0
 Pero primero dividimos nuestra ecuación para L para hacerla más sencilla:
𝑑𝑖 2 𝑑𝑖 𝑅 1
 𝑑𝑡 2
+ 𝑑𝑡 𝐿 + 𝐿∗𝐶 𝑖 = 0
𝑖𝑖−1 −2𝑖𝑖 +𝑖𝑖+1 𝑖𝑖+1 −𝑖𝑖−1 𝑅 1
 + + 𝑖 =0
ℎ2 2ℎ 𝐿 𝐿∗𝐶 𝑖

 Multiplicamos la expresión para 2ℎ2 para simplificar:

𝑅 2ℎ2 𝑖𝑖
 2 𝑖𝑖−1 − 2𝑖1 + 𝑖𝑖+1 + ℎ 𝑖𝑖+1 − 𝑖𝑖 − 1 + =0
𝐿 𝐿∗𝐶
𝑅ℎ 𝑅ℎ 2ℎ2
 2𝑖𝑖−1 − 4𝑖𝑖 + 2𝑖𝑖+1 + 𝑖
𝐿 𝑖+1
− 𝑖
𝐿 𝑖−1
+ 𝐿𝐶 𝑖
𝑖 =0
𝑅ℎ 𝑅ℎ ℎ 2 𝑖𝑖
 𝑖𝑖+1 2 + = 4𝑖𝑖 + 𝑖 − 2𝑖𝑖−1 − 2 =0
𝐿 𝐿 𝑖−1 𝐿𝐶
𝑅ℎ ℎ2 𝑖
4𝑖𝑖 +𝑖𝑖−1 𝐿 −2 −2 𝐿𝐶𝑖
 𝑖𝑖+1 = 𝑅ℎ
2+ 𝐿

𝑅ℎ 𝑖+1 ℎ2 𝑖
4𝑖𝑖+1 +𝑖𝑖 𝐿 −2 −2 𝐿𝐶
 𝑖𝑖+2 = 𝑅ℎ
2+ 𝐿
 𝑅ℎ ℎ2 𝑖𝑖+2
4𝑖𝑖+2 +𝑖𝑖+1 −2 −2
 𝑖𝑖+3 = 𝐿
𝑅ℎ
𝐿𝐶
2+ 𝐿

𝑅ℎ 𝑖+3 ℎ2 𝑖
4𝑖𝑖+3 +𝑖𝑖+2 𝐿 −2 −2 𝐿𝐶
 𝑖𝑖+4 = 𝑅ℎ
2+ 𝐿

𝑅ℎ 𝑖+4 ℎ2 𝑖
4𝑖𝑖+4 +𝑖𝑖+3 𝐿 −2 −2 𝐿𝐶
 𝑖𝑖+5 = 𝑅ℎ
2+ 𝐿

𝑅ℎ 𝑖+5 ℎ2 𝑖
4𝑖𝑖+5 +𝑖𝑖+4 𝐿 −2 −2 𝐿𝐶
 𝑖𝑖+6 = 𝑅ℎ
2+ 𝐿
 𝑡𝑖 𝑖𝑖

𝑡0 0 𝑖0 0

𝑡1 0.5 𝑖1 2.5

𝑡2 1 𝑖2 1.25

𝑡3 1.5 𝑖3 -0.2083

𝑡4 2 𝑖4 -0.5208

𝑡5 2.5 𝑖5 -0.1910

𝑡6 3 𝑖6 0.0781

𝑡7 3.5 𝑖7 0.1027

𝑡8 4 𝑖8 0.0253

𝑡9 4.5 𝑖9 -0.0216

𝑡10 5 𝑖10 -0.0192

𝑡11 5.5 𝑖11 -0.0024

𝑡12 6 𝑖12 0.0052

 Aplicamos la matriz de vandermonde para encontrar nuestro polinomio
deseado, nos ayudamos de Matlab para resolver la matriz:
 𝑝 𝑥 = 𝑎0 + 𝑎1 𝑥1 + 𝑎2 𝑥 2 + 𝑎3 𝑥 3 + 𝑎4 𝑥 4 + 𝑎5 𝑥 5 + 𝑎6 𝑥 6 + 𝑎7 𝑥 7 + 𝑎8 𝑥 8 + 𝑎9 𝑥 9 +
𝑎10 𝑥 10 + 𝑎11 𝑥 11 + 𝑎12 𝑥 12
 RESULTADOS.
 Polinomio aproximado:
 𝑃 𝑥 = 0.000011𝑥 12 + 0.0004𝑥 11 − 0.0064𝑥 10 + 0.0537𝑥 9 − 0.2564𝑥 8 +
0.5915𝑥 7 + 0.2253𝑥 6 − 4.9606𝑥 5 + 10.98𝑥 4 − 5.27𝑥 3 − 10.70𝑥 2 + 10.59𝑥
𝑡 𝑖 𝑟𝑒𝑎𝑙 𝑖 𝑎𝑝𝑟𝑜𝑥.

0 0 0

0.5 2.55 2.5

1 1.6726 1.25

1.5 0.1574 -0.2083

2 -0.5121 -0.5208

2.5 -0.3936 -0.1910

3 -0.0696 0.0781

3.5 0.0992 0.1027

4 0.0906 0.0253

4.5 0.0229 -0.0216

5 -0.0183 -0.0192

5.5 -0.0204 -0.0024

6 -0.0067 0.0052
 CONCLUSIÓN.
 Los métodos numéricos es una herramienta muy importante que nos servirá de mucho
para resolver problemas de la vida cotidiana que son modelados por funciones; y nos
permite aproximarnos lo más posible a ellas ya que resolverlas de la manera tradicional
puede parecer un poco más compleja. En este caso nos ayudó a construir un polinomio
para modelar la curva de la intensidad del circuito con respecto al tiempo.
 El desarrollo de la informática ha aportado mucho a los métodos numéricos ya que sin el
resultaría difícil lograr el nivel actual de utilización de las técnicas numéricas en diferentes
ámbitos.

 RECOMENDACIONES.
 Conocer bien la teoría de los métodos utilizados para después aplicarlos en Matlab.
 El uso de más puntos en la evaluación de la derivada producirá mayor exactitud.
GRACIAS