'1 pset, 1 set 1D x, 1 set 1D y, Legendre 1D polynomial'

def test_linear_fitter_1dlegend(self):

    expected = np.array([[1925.5000000000011, 3444.7500000000005, 1883.2500000000014, 364.4999999999996]]).T
    leg1 = models.Legendre1D(3)
    leg1.parameters = [1, 2, 3, 4]
    y1 = leg1(self.x1)
    pfit = fitting.LinearLSQFitter()
    model = pfit(leg1, self.x1, y1)
    assert_allclose(model.param_sets, expected, atol=(10 ** (-12)))
