`thomson_fitting` error with `lmfit==1.3.0` about plasmapy HOT 3 CLOSED

import numpy as np
from lmfit import Model, Parameters

# Define some data
x = np.linspace(0,1,num=10)
data = 10*x+2



# Parameters list with one parameter
params = Parameters()
params.add("b", value=5, vary=True, min=0.5, max=25)


# Dictionary of constants to be passed to model function through **kws
# of `Model`
settings = {'a':2}


# Define a model function that takes the parameters and the settings keyword
def _model(x, settings=None, **params):
    
    # TypeError is raised here because `settings` keyword is not being 
    # passed through to `_model`, so `settings=None`
    a = settings['a']
    b = params['b']
    
    return x*b + a


model = Model(
     _model,
     independent_vars=["x"],
     settings=settings,
 )


# Try to fit
# This raises an exception
result = model.fit(
    data,
    params=params,
    x=x,
    method="leastsq",
)

print(result.best_values)

_spectral_density_model_lambda = lambda wavelengths, **params: _model(wavelengths, settings=settings, **params)

 return Model(
        _spectral_density_model_lambda,
        independent_vars=["wavelengths"],
        nan_policy="omit",
        )