Continuation of curves using Lagrange Multipliers

This software computes a curve defined implicitly by a function F from R^m to R^n, with n=m+1. That is, we have n equations and n+1 unknowns. Assuming no degeneracy, this deifnes a curves in R^n. Soon more details about the theoretical background.

You will need GSL libraries installed.

How to compile: g++ -o execname cont_LagrangMult.c rk78.c Fun.c lu.c memoria.c -lgsl -lgslcblas

Fun.c is a file that you need to create. There you define the function F and its Jacobian. See the folder "examples" for examples. Remember to properly tune the parameters n and m in the main code and provide an initial guess according to your function.

Note also that the "kernel" function may provide a vector in oposite direction to the desired one. Tune the conditional after the call to "kernel" according to your needs.