function [muu] = mu(A,betai,betac)

muu = A*(betai + betac);