double bar(double v){
    return v+v*v;
}