`energy_power=xx` not working in `gammapy.astro.darkmatter`'s `table_model.plot` about gammapy HOT 6 CLOSED

>     for channel, label, linestyle, linewidth, color in zip(channels, labels, linestyles, linewidths, colors):
>         fluxes.channel = channel
>         fluxes.table_model.plot(
>             energy_bounds=[mDM / 100, mDM],
>             ax=ax,
>             label=label, 
>             linestyle=linestyle,
>             energy_power=2,
>             linewidth=linewidth,
>             color=color,
>             yunits=u.Unit("1/TeV"),

System:

	python_executable      : /Users/mdoro/opt/mambaforge/envs/gammapy-1.1/bin/python 
	python_version         : 3.9.18     
	machine                : x86_64     
	system                 : Darwin     


Gammapy package:

	version                : 1.1        
	path                   : /Users/mdoro/opt/mambaforge/envs/gammapy-1.1/lib/python3.9/site-packages/gammapy 


Other packages:

	numpy                  : 1.24.4     
	scipy                  : 1.10.1     
	astropy                : 5.3.4      
	regions                : 0.7        
	click                  : 8.1.3      
	yaml                   : 6.0.1      
	IPython                : 8.7.0      
	jupyterlab             : 3.5.3      
	matplotlib             : 3.7.3      
	pandas                 : 1.5.2      
	healpy                 : 1.16.6     
	iminuit                : 2.21.3     
	sherpa                 : 4.15.1     
	naima                  : 0.10.0     
	emcee                  : 3.1.3      
	corner                 : 2.2.1      
	ray                    : not installed 

import numpy as np
import matplotlib.pyplot as plt
%matplotlib inline

from gammapy.astro.darkmatter import (
    profiles,
    JFactory,
    PrimaryFlux,
    DarkMatterAnnihilationSpectralModel,
)

from gammapy.maps import WcsGeom, WcsNDMap
from astropy.coordinates import SkyCoord
from matplotlib.colors import LogNorm
from regions import CircleSkyRegion
import astropy.units as u

fig, axes = plt.subplots(1, 3, figsize=(13, 4))
mDMs = [0.5, 5, 50] * u.TeV

labels = [r'$b\bar{b}$', r'$\tau^+\tau^-$',r'$W^+W^-$',r'$\mu^+\mu^-$' ]
linestyles = ['-', '--', 'dotted','dashdot']
linewidths = [2,2,1,1]
colors = ['tab:blue','tab:orange','black','black']
channels = ["b", "tau", "W", "mu"]

# zip it for the loop
for mDM, ax in zip(mDMs, axes):
    fluxes.mDM = mDM
    ax.set_title(rf"m$_{{\mathrm{{DM}}}}$ = {mDM}")
    ax.set_yscale("log")

    ax.set_ylabel("dN/dE")
    
    for channel, label, linestyle, linewidth, color in zip(channels, labels, linestyles, linewidths, colors):
        fluxes.channel = channel
        fluxes.table_model.plot(
            energy_bounds=[mDM / 100, mDM],
            ax=ax,
            label=label, 
            linestyle=linestyle,
            energy_power=2,
            linewidth=linewidth,
            color=color,
            yunits=u.Unit("1/TeV"),
        )

axes[0].legend()
plt.tight_layout()
plt.subplots_adjust(wspace=0.3)
plt.subplots_adjust(hspace=0.3)