def _derivate_diff_eq(listofpoly):
 sol = []
 a = (len(listofpoly) - 1)
 sol.append(DMFdiff(listofpoly[0]))
 for (i, j) in enumerate(listofpoly[1:]):
  sol.append((DMFdiff(j) + listofpoly[i]))
 sol.append(listofpoly[a])
 return sol

def _derivate_diff_eq ( listofpoly ) : sol = [ ] a = ( len ( listofpoly ) - 1 ) sol . append ( DMFdiff ( listofpoly [ 0 ] ) ) for ( i , j ) in enumerate ( listofpoly [ 1 : ] ) : sol . append ( ( DMFdiff ( j ) + listofpoly [ i ] ) ) sol . append ( listofpoly [ a ] ) return sol
let a differential equation a0(x)y(x) + a1(x)y(x) + .