function y = der_f1(x)
    y = 10*cos(10*x) - 3* sin(3*x);
end