lvalues=[10,1,0,-1,-10**3,-10**6,-10**9,-10**12]

for i in range(len(lvalues)):
	print "def f" + str(i) + "(y,x):"
	print "\treturn " + str(lvalues[i]) + "*(y - sin(x)) + cos(x)"
	#print "\t\treturn " + str(i) + "*(y - sin(x)) + cos(x)"

print "\n"
print "def solutions(lmbda,init,xvalues,tol):"
for i in range(len(lvalues)):
	print "\tif(lmbda == " + str(lvalues[i]) + "):"
	print "\t\treturn odeint(f" + str(i) + ",init,xvalues,rtol=tol,full_output=1)"