function [fevd_record,fevd_estimates]=panel5fevd(N,n,struct_irf_record,gamma_record,It,Bu,IRFperiods,FEVDband)






% run the Gibbs sampler to obtain FEVD draws
[fevd_record]=bear.fevd(struct_irf_record,gamma_record,It,Bu,N*n,IRFperiods,FEVDband);
% obtain point estimates and credibility intervals
[fevd_estimates]=bear.fevdestimates(fevd_record,N*n,IRFperiods,FEVDband);




















