Fixed major error: Forgot that 'L[i+1] = L[i] + ...' etc. in differential equations. Write to file instead of keeping all arrays. Fixed a lambda which is supposed to have a temperature treshhold. Set better initial values. Print stuff for debugging.