'Plot graphical display of FEVD
Parameters
periods : int, default None
Defaults to number originally specified. Can be at most that number'

def plot(self, periods=None, figsize=(10, 10), **plot_kwds):

    import matplotlib.pyplot as plt
    k = self.neqs
    periods = (periods or self.periods)
    (fig, axes) = plt.subplots(nrows=k, figsize=figsize)
    fig.suptitle('Forecast    error    variance    decomposition    (FEVD)')
    colors = [str(c) for c in (np.arange(k, dtype=float) / k)]
    ticks = np.arange(periods)
    limits = self.decomp.cumsum(2)
    for i in range(k):
        ax = axes[i]
        this_limits = limits[i].T
        handles = []
        for j in range(k):
            lower = (this_limits[(j - 1)] if (j > 0) else 0)
            upper = this_limits[j]
            handle = ax.bar(ticks, (upper - lower), bottom=lower, color=colors[j], label=self.names[j], **plot_kwds)
            handles.append(handle)
        ax.set_title(self.names[i])
    (handles, labels) = ax.get_legend_handles_labels()
    fig.legend(handles, labels, loc='upper    right')
    plotting.adjust_subplots(right=0.85)
