function I = F(V)
global a b c

%I = a*V.*(V-b).*(V-c)+d;

I=a*V.*(V.^2+b*V+c);