Diferencia entre revisiones de «Reacciones complejas - Grupo 1 A»
De MateWiki
(→Apartado 3) |
(→Apartado 3) |
||
| Línea 16: | Línea 16: | ||
%PVI por Euler con una h = 0.1 | %PVI por Euler con una h = 0.1 | ||
clear all | clear all | ||
| − | + | ||
%condiciones dadas | %condiciones dadas | ||
y0 = 0; | y0 = 0; | ||
| Línea 30: | Línea 30: | ||
ya = zeros(1,N+1); | ya = zeros(1,N+1); | ||
yb = zeros(1,N+1); | yb = zeros(1,N+1); | ||
| − | |||
| − | |||
| − | |||
| − | |||
| + | |||
| + | %solucion obtenida | ||
| + | y =(3 - exp(2.*t))./(1 - exp(2.*t)); | ||
| + | |||
| + | %Concentracion de los reactivos | ||
ya = 3.* t; | ya = 3.* t; | ||
yb = 1.* t; | yb = 1.* t; | ||
| + | |||
%solucion grafica | %solucion grafica | ||
Revisión del 13:52 12 feb 2015
| Trabajo realizado por estudiantes | |
|---|---|
| Título | Reacciones complejas. Grupo 1A |
| Asignatura | Ecuaciones Diferenciales |
| Curso | Curso 2014-15 |
| Autores | María Ramírez
Ignacio Posada Antonio López-Mateos Pablo Bueno |
| Este artículo ha sido escrito por estudiantes como parte de su evaluación en la asignatura | |
Contenido
1 Introducción
2 Apartado 1
3 Apartado 2
4 Apartado 3
%PVI por Euler con una h = 0.1
clear all
%condiciones dadas
y0 = 0;
t0 = 0;
tN = 2;
h = 0.1;
%calculamos N
N = (tN - t0)/h;
%definimos la variable t
t = linspace(t0,tN,N+1);
%funcion y
y = zeros(1,N+1);
ya = zeros(1,N+1);
yb = zeros(1,N+1);
%solucion obtenida
y =(3 - exp(2.*t))./(1 - exp(2.*t));
%Concentracion de los reactivos
ya = 3.* t;
yb = 1.* t;
%solucion grafica
hold on
plot(t,y,'linewidth',2)
plot(t,ya,'g','linewidth',2)
plot(t,yb,'r','linewidth',2)
legend('Concentracion de C','Concentracion de A','Concentracion de B','Location','best');
hold off
Gráfica de las dos concentraciones y el producto. Utilizando el método de Euler y con una h = 0.1
