'This case covers:
N param sets , 1 set 1D x --> N 1D y data'

def test_non_linear_NYset(self):

    g1 = models.Gaussian1D([10, 10], [3, 3], [0.2, 0.2], n_models=2)
    y1 = g1(self.x1, model_set_axis=False)
    assert np.all(((y1[0, :] - y1[1, :]).nonzero() == np.array([])))
