add F function
add am and qr functions
plot initial values of coefficients
breaking for day. Need to fix recursion loop.