function y1 = f1(x,y)
    y1 = -2*x^3 + 12*x^2 - 20*x + 8.5;
end