Diferencia entre revisiones de «Desintegración Radiactiva.Grupo 5»
De MateWiki
| Línea 4: | Línea 4: | ||
[[Categoría:ED14/15]] | [[Categoría:ED14/15]] | ||
[[Categoría:Trabajos 2014-15]] | [[Categoría:Trabajos 2014-15]] | ||
| + | |||
| + | ==Apartado 1== | ||
| + | {{matlab|codigo=clear all | ||
| + | |||
| + | %% TRABAJO 3 APARTADO 1 | ||
| + | % M'(t)=-1,24*10^(-4)*M(t) | ||
| + | % M(0)=1 | ||
| + | |||
| + | % M(t) representa la cantidad de C14 en un instante dado, como no nos | ||
| + | % dan una cantidad inicial M0 supondremos que esta es uno. Nuestro objetivo | ||
| + | % es determinar en que intanste queda un 8% de la cantidad inicial M0.7 | ||
| + | % Como M0=1, es como si trabajaramos todo el rato en tanto por uno | ||
| + | % Después demostraremos que quedará un 8% de M0 pasado el tiempo que | ||
| + | % tratamos de determinar independientemente del M0 adoptado. | ||
| + | |||
| + | t0=0; | ||
| + | h=input('Inserte el valor del paso, porfavor: '); | ||
| + | M0=input('Inserte la cantidad inicial de carbono 14, porfavor: '); | ||
| + | t(1)=t0; | ||
| + | M(1)=M0; | ||
| + | |||
| + | % Euler explícito | ||
| + | i=1; | ||
| + | while M(i)>(0.08*M0) | ||
| + | M(i+1)=M(i)+h*((-1.24*10^(-4))*M(i)); | ||
| + | t(i+1)=t(i)+h; | ||
| + | i=i+1; | ||
| + | end | ||
| + | |||
| + | disp('Tiempo final:') | ||
| + | disp(t(end)) | ||
| + | plot(t,M) | ||
| + | xlabel('Cantidad de Carbono 14'); | ||
| + | ylabel('Tiempo (años)'); | ||
| + | legend('Euler explícito','Location','best'); | ||
| + | }} | ||
Revisión del 00:02 25 feb 2015
| Trabajo realizado por estudiantes | |
|---|---|
| Título | Desintegración Radiactiva. Grupo 5. |
| Asignatura | Ecuaciones Diferenciales |
| Curso | Curso 2014-15 |
| Autores | Álvaro Ramón López, Diego García Vaquero, Noemí Palomino Bustos, Mercedes Ruiz Barrajón, Teresa Quintana Romero, Araceli Martín Candilejo. |
| Este artículo ha sido escrito por estudiantes como parte de su evaluación en la asignatura | |
Apartado 1
clear all
%% TRABAJO 3 APARTADO 1
% M'(t)=-1,24*10^(-4)*M(t)
% M(0)=1
% M(t) representa la cantidad de C14 en un instante dado, como no nos
% dan una cantidad inicial M0 supondremos que esta es uno. Nuestro objetivo
% es determinar en que intanste queda un 8% de la cantidad inicial M0.7
% Como M0=1, es como si trabajaramos todo el rato en tanto por uno
% Después demostraremos que quedará un 8% de M0 pasado el tiempo que
% tratamos de determinar independientemente del M0 adoptado.
t0=0;
h=input('Inserte el valor del paso, porfavor: ');
M0=input('Inserte la cantidad inicial de carbono 14, porfavor: ');
t(1)=t0;
M(1)=M0;
% Euler explícito
i=1;
while M(i)>(0.08*M0)
M(i+1)=M(i)+h*((-1.24*10^(-4))*M(i));
t(i+1)=t(i)+h;
i=i+1;
end
disp('Tiempo final:')
disp(t(end))
plot(t,M)
xlabel('Cantidad de Carbono 14');
ylabel('Tiempo (años)');
legend('Euler explícito','Location','best');
