def prde_spde(a, b, Q, n, DE):
 (R, Z) = list(zip(*[gcdex_diophantine(b, a, qi) for qi in Q]))
 A = a
 B = (b + derivation(a, DE))
 Qq = [(zi - derivation(ri, DE)) for (ri, zi) in zip(R, Z)]
 R = list(R)
 n1 = (n - a.degree(DE.t))
 return (A, B, Qq, R, n1)

def prde_spde ( a , b , Q , n , DE ) : ( R , Z ) = list ( zip ( * [ gcdex_diophantine ( b , a , qi ) for qi in Q ] ) ) A = a B = ( b + derivation ( a , DE ) ) Qq = [ ( zi - derivation ( ri , DE ) ) for ( ri , zi ) in zip ( R , Z ) ] R = list ( R ) n1 = ( n - a . degree ( DE . t ) ) return ( A , B , Qq , R , n1 )
special polynomial differential equation algorithm: parametric version .