docs: def formatLij(Lij0,Ne):
  	"""This function transforms a list of laser conections of the form
  	[i,j,[l1,l2,...]] between states i and j by lasers l1,l2,... into a
  	Ne x Ne matrix whose elements are the lasers connecting the corresponding
  	indices."""
  	#We create Lij as a matrix of lists of laser indices
  	global Lij
  	Lij=[]
  	for i in range(Ne):
  		fila=[]
  		for j in range(Ne):
  			band=False
  			for tri in Lij0:
  				if [i+1,j+1] == tri[:2]:
  					band=True
  					break
  				elif [j+1,i+1] == tri[:2]:
  					band=True
  					break
  			if band:
  				fila+=[tri[2]]
  			else:
  				fila+=[[]]
  		Lij+=[fila]
  	return Lij