def METHOD_NAME(L,M,timesL, timesM,root,scale=100.,s_species=1, maxSP=1000):
	ts=list()
	te=list()
	l_t=L[0]
	m_t=M[0]
	if timesL==0: timesL=np.array([-root,0])
	if timesM==0: timesM=np.array([-root,0])
	L,M,root=L/scale,M/scale,int(root*scale)
	for i in range(s_species): 
		ts.append(root)
		te.append(0)
	for t in range(root,0): # time
		for j in range(len(timesL)-1):
			if -t/scale<=timesL[j] and -t/scale>timesL[j+1]: l=L[j]
		for j in range(len(timesM)-1):
			if -t/scale<=timesM[j] and -t/scale>timesM[j+1]: m=M[j]
		TE=len(te)
		for j in range(TE): # extant lineages
			if te[j]==0:
				ran=random.random()
				if ran<l: 
					te.append(0) # add species
					ts.append(t) # sp time
				elif ran>l and ran < (l+m): # extinction
					te[j]=t
	te=array(te)
	return -np.array(ts)/scale, -np.array(te)/scale