A function that computes the impulse responses would be a welcome enhancement. Will be implemented in the new version.