Reacciones con autocatálisis, Grupo C26 - Historial de revisiones

Alejandropcr en 20:23 6 mar 2015

2015-03-06T20:23:00Z

2015-03-06T20:18:32Z

2015-03-06T20:16:31Z

2015-03-06T20:09:17Z

2015-03-06T20:06:15Z

2015-03-06T20:05:16Z

2015-03-06T20:00:08Z

2015-03-06T19:58:24Z

2015-03-06T19:56:13Z

2015-03-06T19:54:08Z

@@ Línea 264: / Línea 264: @@
 [[Archivo:cod2.png]]
@@ Línea 465: / Línea 468: @@
 [[Archivo:cdg5.png]]
 {{Matlab|codigo=%Introducimos los datos

@@ Línea 465: / Línea 465: @@
 [[Archivo:cdg5.png]]
 [[Archivo:cod7.png]]
 [[Categoría:Ecuaciones Diferenciales]]
 [[Categoría:ED14/15]]
 [[Categoría:Trabajos 2014-15]]

@@ Línea 466: / Línea 466: @@
 [[Archivo:cdg5.png]]
 [[Categoría:Ecuaciones Diferenciales]]
 [[Categoría:ED14/15]]
 [[Categoría:Trabajos 2014-15]]

← Revisión anterior		Revisión del 20:09 6 mar 2015
Línea 263:		Línea 263:

	[[Archivo:cod2.png]]		[[Archivo:cod2.png]]
		+
		+
		+

	'''5. REACCIÓN DE LOTKA. INTERPRETACIÓN'''		'''5. REACCIÓN DE LOTKA. INTERPRETACIÓN'''

← Revisión anterior		Revisión del 20:06 6 mar 2015
Línea 261:		Línea 261:
	legend('Concentración de B','Concentración de A','Location','Best');}}		legend('Concentración de B','Concentración de A','Location','Best');}}

		+
		+	[[Archivo:cod2.png]]

	'''5. REACCIÓN DE LOTKA. INTERPRETACIÓN'''		'''5. REACCIÓN DE LOTKA. INTERPRETACIÓN'''

← Revisión anterior		Revisión del 20:05 6 mar 2015
Línea 219:		Línea 219:

	[[Archivo:cooood1.png]]		[[Archivo:cooood1.png]]
		+
		+	'''4.2 METODO DE RUNGE-KUTTA'''
		+
		+	El método Runge-Kutta de orden 4 es de orden superior al de Euler, con lo cual nos ofrecerá una mayor aproximación a la solución real. El código del programa es el que se muestra a continuación.
		+
		+	{{Matlab\|codigo=%metemos los datos de k, h y creamos el vector t
		+	k=1;
		+	t0=0;
		+	tN=10;
		+	h=0.1;
		+	y0=0.01;
		+	x0=1;
		+	t=t0:h:tN;
		+	%creamos los vectores x e y.
		+	y=zeros(1,length(t));
		+	x=y;
		+	%Introducimos los valores iniciales de las concentraciones
		+	y(1)=y0;
		+	x(1)=x0;
		+	%Definimos las funciones
		+	F=inline('x*y','t','y','x');
		+	G=inline('-x*y','t','y','x');
		+	%Aplicamos el método Runge Kutta de orden 4
		+	for k=1:length(t)-1
		+	K1y=F(t(k),y(k),x(k));
		+	K1x=G(t(k),y(k),x(k));
		+	K2y=F(t(k)+(h/2),y(k)+hK1y/2,x(k)+hK1x/2);
		+	K2x=G(t(k)+(h/2),y(k)+hK1y/2,x(k)+hK1x/2);
		+	K3y=F(t(k)+(h/2),y(k)+hK2y/2,x(k)+hK2x/2);
		+	K3x=G(t(k)+(h/2),y(k)+hK2y/2,x(k)+hK2x/2);
		+	K4y=F(t(k)+h,y(k)+K3yh,x(k)+K3xh);
		+	K4x=G(t(k)+h,y(k)+K3yh,x(k)+K3xh);
		+	y(k+1)=y(k)+(h/6)(K1y+2K2y+2*K3y+K4y);
		+	x(k+1)=x(k)+(h/6)(K1x+2K2x+2*K3x+K4x);
		+	end
		+	%Dibujamos ambas gráficas.
		+	hold on
		+	plot(t,y);
		+	plot(t,x,'r');
		+	hold off
		+	legend('Concentración de B','Concentración de A','Location','Best');}}