On computational methods for Gramians of continuous dynamic systems