double fun(double x) {
    if (x < threshold)
        return k*x^2;
    return 50;
}