function val = radio_burbujas(t,y)
% PVI para describir radio de burbujas de aire microsegundos despues de
% comenzado un tratamiento por LEC
% La ecuacion diferencial original es de 2do orden y se convierte 
% a un sistema de ecuaciones diferenciales de 1er orden
% la parte derecha f(t,y(t)) del sistema resultante se evalua en esta
% funcion

% introduciendo z=r', se tiene el siguiente sistema de EDOs de 1er orden
% r' = z
% z' = (r^(-3*gamma)-1-(3/2)*z^2)/r
% el vector de entrada y = (r,z)
gamma = 1.4;
val = [y(2);...
    (y(1)^(-3*gamma)-1-(3/2)*y(2)^2)/y(1)];
