In this one, you have inserted a wrong sign:
function dydt = fSystem(t, y, beta1, alpha2, beta2, p, psi, T)
f2 = -alpha2*y(2)+beta2+y(1)*y(2); %% ERR: 'beta2+' has to be 'beta2*'
g = beta1*y(2)+p*f2*y(1)^(-1)-psi*(T*y(1))^(-1);
dydt(1) = g*y(1)-beta1*y(1)*y(2);
dydt(2) = f2;
dydt = dydt';
end
