original_docstring: 
      Initialize the value of a variable inside a 
      force constant matrix, assuming the variable name
      has the form 'Cnab' where n is an integer identifying
      a set of independent variables and a and b identify 
      the directions (1,2 or 3).
  
      Return a random number between 5 and 10 for the ab=33 component,
      a random number between 1 and 4 for the 11 and 22 components,
      and a random number between -1 and 1 for the off-diagonal components.
  
      # Note: these values are expressed in the default GUI units that are 10^19 N/m^3
      
docstring: Initialize the value of a variable inside a
  force constant matrix, assuming the variable name
  has the form 'Cnab' where n is an integer identifying
  a set of independent variables and a and b identify
  the directions (1,2 or 3).
  
  Return a random number between 5 and 10 for the ab=33 component,
  a random number between 1 and 4 for the 11 and 22 components,
  and a random number between -1 and 1 for the off-diagonal components.
  
  these values are expressed in the default GUI units that are 10^19 N/m^3
docstring_tokens: Initialize the value of a variable inside a force constant matrix assuming the variable name has the form ' Cnab ' where n is an integer identifying a set of independent variables and a and b identify the directions ( 1 2 or 3 ) . Return a random number between 5 and 10 for the ab = 33 component a random number between 1 and 4 for the 11 and 22 components and a random number between - 1 and 1 for the off - diagonal components . these values are expressed in the default GUI units that are 10^19 N / m^3
code: def matrix_initialization(variable):
      """
      Initialize the value of a variable inside a 
      force constant matrix, assuming the variable name
      has the form 'Cnab' where n is an integer identifying
      a set of independent variables and a and b identify 
      the directions (1,2 or 3).
  
      Return a random number between 5 and 10 for the ab=33 component,
      a random number between 1 and 4 for the 11 and 22 components,
      and a random number between -1 and 1 for the off-diagonal components.
  
      # Note: these values are expressed in the default GUI units that are 10^19 N/m^3
      """
      # Range of value for 11 and 22
      min_11 = 2.0
      max_11 = 4.5
      # Range of random values for 33
      min_33 = 8.0
      max_33 = 13.0
      # range of values for off-diagonal
      min_offdiag = -0.8
      max_offdiag = 0.8
  
      if variable[-2:] == "33":
          return np.round(min_33 + (max_33 - min_33) * np.random.rand(), 1)
      if variable[-1] == variable[-2]:
          return np.round(min_11 + (max_11 - min_11) * np.random.rand(), 1)
      return np.round(min_offdiag + (max_offdiag - min_offdiag) * np.random.rand(), 1)
code_tokens: def matrix_initialization ( variable ) : min_11 = 2.0 max_11 = 4.5 min_33 = 8.0 max_33 = 13.0 min_offdiag = - 0.8 max_offdiag = 0.8 if variable [ - 2 : ] == "33" : return np . round ( min_33 + ( max_33 - min_33 ) * np . random . rand ( ) , 1 ) if variable [ - 1 ] == variable [ - 2 ] : return np . round ( min_11 + ( max_11 - min_11 ) * np . random . rand ( ) , 1 ) return np . round ( min_offdiag + ( max_offdiag - min_offdiag ) * np . random . rand ( ) , 1 )
short_docstring: Initialize the value of a variable inside a
  force constant matrix, assuming the variable name
  has the form 'Cnab' where n is an integer identifying
  a set of independent variables and a and b identify
  the directions (1,2 or 3).
short_docstring_tokens: Initialize the value of a variable inside a force constant matrix assuming the variable name has the form ' Cnab ' where n is an integer identifying a set of independent variables and a and b identify the directions ( 1 2 or 3 ) .
comment: """
      Initialize the value of a variable inside a 
      force constant matrix, assuming the variable name
      has the form 'Cnab' where n is an integer identifying
      a set of independent variables and a and b identify 
      the directions (1,2 or 3).
  
      Return a random number between 5 and 10 for the ab=33 component,
      a random number between 1 and 4 for the 11 and 22 components,
      and a random number between -1 and 1 for the off-diagonal components.
  
      # Note: these values are expressed in the default GUI units that are 10^19 N/m^3
      """$SEP$# Range of value for 11 and 22$SEP$# Range of random values for 33$SEP$# range of values for off-diagonal