Diferencia entre revisiones de «Modelos de mezclas (Grupo 15)»

De MateWiki
Saltar a: navegación, buscar
(Página creada con «==Resolución de los sistemas: Método de Euler== Sistema (1) {{matlab|codigo= clear all close all format long t0=0;tN=6*100/4.5; N=100;h=(tN-t0)/N; t=t0:h:tN; x0=20;y0...»)
 
(Resolución de los sistemas: Método de Euler)
Línea 1: Línea 1:
==Resolución de los sistemas: Método de Euler==
 
  
Sistema (1)
 
 
 
{{matlab|codigo=
 
 
clear all
 
close all
 
format long
 
t0=0;tN=6*100/4.5;
 
N=100;h=(tN-t0)/N;
 
t=t0:h:tN;
 
x0=20;y0=0;
 
x(1)=x0;y(1)=y0;
 
for n=1:N
 
x(n+1)=x(n)+h*(-3/100)*x(n);
 
y(n+1)=y(n)+h*(3/100)*x(n)-h*(4.5/100)*y(n);
 
end
 
X=20*exp(-3*t/100); %exacta
 
Y=40*(exp(-3*t/100)-exp(-4.5*t/100)); %exacta
 
figure(100)
 
plot(t,x,'b',t,y,'r',t,X,'.-b',t,Y,'.-r')
 

Revisión del 14:16 3 mar 2013