Diferencia entre revisiones de «Usuario:Carlos Castro»

Revisión del 11:20 31 ene 2013

y0 = f (t, y)=y(1-y);

y(t0) = y0:

y(0) = 1/10

y(n+1) = y(n) + h y(n)(1 - y(n))

t0=0; tN=4;

y0=1/10;

N=40; h=(tN-t0)/40;

yy=y0;

y(1)=yy;

for n=1:N-1

yy=yy+h*yy*(1-yy);

y(n+1)=yy;

end

x=t0:h:tN;

plot(x,y,'x');

@@ Línea 1: / Línea 1: @@
+===== Resolución de la Ecuación logística por el método de Euler =====
+* Logistic equation:
+:y0 = f (t, y)=y(1-y); <br />
+:y(t0) = y0:
+* Numerical scheme
+: y(0) = 1/10 <br />
+: y(n+1) = y(n) + h y(n)(1 - y(n))
+*  MATLAB code
+: t0=0; tN=4; <br />
+: y0=1/10; <br />
+: N=40; h=(tN-t0)/40; <br />
+:yy=y0; <br />
+: y(1)=yy; <br />
+: for n=1:N-1 <br />
+: yy=yy+h*yy*(1-yy); <br />
+: y(n+1)=yy; <br />
+: end <br />
+: x=t0:h:tN; <br />
+: plot(x,y,'x'); <br />
+<gallery>
+Archivo:solucion.jpg|Aproximación numérica
+Archivo:error.jpg|Error entre la solución exacta y la aproximación numérica
+</gallery>