M. Euler (Ecuación Integro Diferencial)
Publicado por Juan José (60 intervenciones) el 17/03/2015 08:49:45
Buenas chicos,
Voy a pediros ayuda para la solución de un problema utilizando el Método de Euler en una ecuación donde mezclan el método en sí, con una integral que tratan de aproximar mediante Regla del Trapecio.
El enunciado del problema lo subo en fotografia.
El código que yo he preparado lo pego aqui: (Funciona pero aparecen muchos Infinitos en el Resultado)
He subido la EJECUCIÓN del Código en txt
Muchísimas Gracias por vuestro tiempo.
Saludos!!
Voy a pediros ayuda para la solución de un problema utilizando el Método de Euler en una ecuación donde mezclan el método en sí, con una integral que tratan de aproximar mediante Regla del Trapecio.
El enunciado del problema lo subo en fotografia.
El código que yo he preparado lo pego aqui: (Funciona pero aparecen muchos Infinitos en el Resultado)
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
function E=eulerEj10(f,a,b,ya,M)
%Input - f is the function entered as a string 'f'
% - a and b are the left and right endpoints
% - ya is the initial condition y(a)
% - M is the number of steps
%Output - E=[T' Y'] where T is the vector of abscissas and
% Y is the vector of ordinates
% NUMERICAL METHODS: MATLAB Programs
%(c) 1999 by John H. Mathews and Kurtis D. Fink
%To accompany the textbook:
%NUMERICAL METHODS Using MATLAB,
%by John H. Mathews and Kurtis D. Fink
%ISBN 0-13-270042-5, (c) 1999
%PRENTICE HALL, INC.
%Upper Saddle River, NJ 07458
h=(b-a)/M;
T=zeros(1,M+1);
Y=zeros(1,M+1);
T=a:h:b;
Y(1)=ya; % k=0;
Y(2)=Y(1)+h*(1.3*Y(1)-0.25*(Y(1))^2-0.0001*Y(1)*(0+(h/2)*(0+Y(1)))); % k=1;
I(1)=0; % k=0;
I(2)=I(1)+(h/2)*(0+Y(1)); % k=1;
for j=2:M
I(j)=I(j-1)+(h/2)*(Y(j-1)+Y(j));
Y(j+1)=Y(j)+h*(1.3*Y(j)-0.25*(Y(j))^2-0.0001*Y(j)*(I(j)));
end
E=[T' Y'];
He subido la EJECUCIÓN del Código en txt
Muchísimas Gracias por vuestro tiempo.
Saludos!!
- Ejercicio-Ejecutado-10.zip(402,0 B)
Valora esta pregunta
0
