
function R = myfun1 (t, y)
emit = 60*10^-6;
kave = 20.3135;
kb = 0.0015;
R = [y(2) ; -kave*y(1)+ kb/y(1)+ emit^2/y(1)^3];
end