Diferencia entre revisiones de «Usuario:Carlos Castro»
De MateWiki
| Línea 7: | Línea 7: | ||
: y(n+1) = y(n) + h y(n)(1 - y(n)) | : y(n+1) = y(n) + h y(n)(1 - y(n)) | ||
* MATLAB code | * MATLAB code | ||
| − | + | {{matlab|codigo= | |
| − | + | 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'); | ||
| + | }} | ||
<gallery> | <gallery> | ||
Revisión del 11:26 31 ene 2013
Resolución de la Ecuación logística por el método de Euler
- Logistic equation:
- y0 = f (t, y)=y(1-y);
- y(t0) = y0:
- Numerical scheme
- y(0) = 1/10
- y(n+1) = y(n) + h y(n)(1 - y(n))
- MATLAB code
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');
Aproximación numérica
Error entre la solución exacta y la aproximación numérica
