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