ejercicio3

e3-pyl-g1
mail@pastecode.io avatar
unknown
matlab
2 years ago
333 B
3
Indexable
Never
clc; clear;
a=0.1; %parámetro a
b=0.2; %parámetro b

%condiciones iniciales en el vector x0
x0=[100,0]; %100 es el valor inicial de x, 0 es valor inicial de y
tf=20; %tiempo final


fg=@(t,x) [-a*x(1);a*x(1)-b*x(2)]; 
[t,x]=ode45(fg,[0,tf],x0);
plot(t,x)
grid on
xlabel('t')
ylabel('x,y');
title('Serie radioactiva')