function I = F(V)
global a b c d
I = a*V.*(V-b).*(V-c)+d;