Matlab - decaimiento radiactivo

hola profesor, buen dia, estoy haciendo un programa que me compara con nuemros aleatorios si un nucleo decae o no, utilizo el bucle for y sentencia if, ya lo he conseguido obtener los datos, pero ahora quiero graficarlos en función del tiempo,
aquí esta el scrip
disp('***********************************************')
no=input('numero de nucleos = ');
tf=input('time = ')
y=no;
for t=0:1:tf
b=0.01*t;
a=rand(1);
if a<b
y=y-1;
else
y=y;
end
x=y;
disp(x)
end