shockvector

It should be possible to get the final shock values for each firm. I.e. both the h_d(T) and h_u(T). The R code does this with the option track_h. We should be able to do the same. Biggest hurdle: It's an enormous matrix that is not easy to put into a text file. Other solutions might be necessary.