Hi there I'm sure this is an embarrassingly simple problem but...
I've been using the shootf.h routine to find valid solutions to the 1D radial Schrödinger equation. Since the wavefunctions will be normalised I only want to match the logarithmic derivates at the match point, but I've been unable to create a Score function that does this, without returning on the first iteration or looping indefinitely. Any suggestions would be greatly appreciated.