function y = model1(c, x)
    y = c(1) * sin(x) + c(2) * x.^2 + c(3) * x.*sin(x);
end