Some code examples:
Compute Laplace coefficients (could be improved to better cover alpha to 1 limit)
laplace.c
Convert orbital elements to cartesian coordinates and vice versa. Compute f,g functions.
kepcart.cpp
kepcart.h