Metodo de euler-matlab

Solo disponible en BuenasTareas
  • Páginas : 2 (282 palabras )
  • Descarga(s) : 0
  • Publicado : 26 de diciembre de 2011
Leer documento completo
Vista previa del texto
SOLUCIÓN ESCUACIONES DIFERENCIALES
Problema del valor inicial
1. Método de Euler
yi+1=yi+hF(xi,yi)

h=(xf-x0)n

i |xi | yi | F(xi,yi) | yi+1=yi+hF(xi,yi) |
i=0 | x0=0 | y0=1 | Fx0,y0=-y0+x0+1=-1+0+1=0 | y1=y0+hFx0,y0=1+0.2*0=1 |
i=1 |x1=0.2 | y1=1 | Fx1,y1=-y1+x1+1=-1+0.2+1=0.2 | y2=y1+hFx1,y1=1+0.2*0.2=1.04 |
i=2 | x2=0.4 | y2=1.04 |Fx2,y2=-y2+x2+1=-1.04+0.4+1=0.36 | y3=y2+hFx2,y2=1.04+0.2*0.36=1.12 |
i=3 | x3=0.6 | y3=1.112 | Fx3,y3=0.488 | y4=1.2096 |
i=4 | x4=0.8 | y4=1.2096| Fx4,y4=0.5904 | y5=1.3277 |
i=5 | x5=1 | | | |

Algoritmo
* Ingreso F(x,y), x0, yf, n, y0
* Ingresoh=(xf-x0)n
* x0: h :yf
* Evaluar F(xi,yi)
* y=y0+hF(xi,yi)
* Actual y0=y ir a paso y



programa
%eulerclear all
x0=input('ingrese x inicial:');
xf=input('ingrese x final:');
n=input('ingrese el numero de divisiones:');y0=input('ingrese el valor de y0: ');
f=input('ingrese f:');
h=(xf-x0)./n;
x=x0:h:xf;
disp('i xi yi F(xi,yi) yi+1')disp('----------------------------------------')
y(1)=y0
for i=1:n
F=feval(f,x(i),y0);
y=y0+h*F;%euler
fprintf('%2.0f %4.2f %6.4f %6.4f %6.4f\n',i,x0,y0,F,y)
y0=y;
end
Desarrollo
ingrese x inicial:0
ingrese x final:1
tracking img