ejercicio3

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')