added new library to solver vertical structure equation with rho(z)