This paper is a generalization and extension of earlier research in pa
rametric optimal power flow. Its principal features are: (i) the use o
f a full non-linear OPF model, (ii) a more general parameterization of
equalities, inequalities and objective function and (iii) an algorith
m that exactly tracks the OPF behaviour in terms of the continuation p
arameter. The parameterization allows us to distinguish between two ph
ases. One serves to find the OPF solution to a static problem for a fi
xed load and network starting from an arbitrary initial condition. The
second, finds the trajectories corresponding to varying loads. The op
timal trajectories, in both phases, offer an excellent visualization o
f the complex nature of the OPF solutions, that is, the highly non-lin
ear behaviour and the sensitivity of the solutions to parameter variat
ions.