'A case of 1 set of x and multiple sets of y and parameters.'

def test_wrong_pset(self):

    expected = np.array([[1.0, 0], [1, 1], [1, 2], [1, 3], [1, 4], [1, 5]])
    p1 = models.Polynomial1D(5, n_models=2)
    params = {}
    for i in range(6):
        params[p1.param_names[i]] = [1, i]
    p1 = models.Polynomial1D(5, model_set_axis=0, **params)
    y1 = p1(self.x1, model_set_axis=False)
    pfit = fitting.LinearLSQFitter()
    model = pfit(p1, self.x1, y1)
    assert_allclose(model.param_sets, expected, atol=(10 ** (-7)))
