import?math
import?numpy
def?LUDecompose?（table）:
????#table?that?contains?our?data?#table?has?to?be?a?square?array?so?we?need?to?check?first
????rowscolumns=numpy.shape（table）
????L=numpy.zeros（（rowscolumns））
????U=numpy.zeros（（rowscolumns））
????if?rows!=columns:
????????return
????for?i?in?range?（columns）:
????????for?j?in?range（i-1）:
????????????sum=0
????????????for?k?in?ran