Scribd
Cargar
36 vistas2 páginas

Integrales Impropias

El documento describe la implementación del método de Euler para resolver una ecuación diferencial ordinaria de primer orden con una condición inicial específica. Se utiliza MATLAB para calcular y graficar la solución en el intervalo [0,2] con un paso de h=0.1. Además, se presenta una verificación simbólica de la solución utilizando la función dsolve de MATLAB.

Cargado por

William Steven Triana Garcia
Derechos de autor
© © All Rights Reserved
Nos tomamos en serio los derechos de los contenidos. Si sospechas que se trata de tu contenido, reclámalo aquí.
Formatos disponibles
Descarga como PDF, TXT o lee en línea desde Scribd
36 vistas2 páginas

Integrales Impropias

El documento describe la implementación del método de Euler para resolver una ecuación diferencial ordinaria de primer orden con una condición inicial específica. Se utiliza MATLAB para calcular y graficar la solución en el intervalo [0,2] con un paso de h=0.1. Además, se presenta una verificación simbólica de la solución utilizando la función dsolve de MATLAB.

Cargado por

William Steven Triana Garcia
Derechos de autor
© © All Rights Reserved
Nos tomamos en serio los derechos de los contenidos. Si sospechas que se trata de tu contenido, reclámalo aquí.
Formatos disponibles
Descarga como PDF, TXT o lee en línea desde Scribd

%Andrés Castro 15 de noviembre 2021

clc;
clear all;
%Hallar la solución por el método de Euler, de la EDO de primer orden y'-y-x
+x^2-
%con condición inicial y(0)=1
%En el intervalo [0,2] con h=0.1
fxy='1+x-y'

a=0;
b=1;
h=0.1;
x0=0;
y0=0

euler(fxy,a,b,x0,y0,h)

%comprobación
syms y(x) %manera de ingresar la función simbólica indicando la dependencia
de y con x. En caso de que no procese la solución, cambiar a syms x%edo=
diff(y)==y+x-x^2+1; %se ingresa la ecuación diferencial de la forma
y'=f(x,y). Se debe cambiar para cada caso.
ini=y(0)==1; %se ingresa la condición inicial. Debe modificarse de acuerdo a
cada problema.
edo= diff(y)==1+x-y; %se ingresa la ecuación diferencial de la forma
y'=f(x,y). Debe modificarse de acuerdo a cada problema.
y=dsolve(edo,ini) %función de Matlab para resolver la EDO. En caso de que no
procese la solución, ingresarla de forma manual

syms x
y= x + exp(-x); %se escribe la solución. Debe modificarse de acuerdo a cada
problema.
ezplot(y,a,b), grid on, hold on %se grafica la solución de la ecuación

fxy =

'1+x-y'

y0 =

tabla =

0 0
0.1000 0.1000
0.2000 0.2000
0.3000 0.3000

1
0.4000 0.4000
0.5000 0.5000
0.6000 0.6000
0.7000 0.7000
0.8000 0.8000
0.9000 0.9000
1.0000 1.0000

y =

x + exp(-x)

Published with MATLAB® R2021b

También podría gustarte