#NEWLINE# def test_normal_denom( ) : #NEWLINE# #INDENT# DE = DifferentialExtension ( extension = { #NEWLINE# #INDENT# 'D' : [ Poly ( 1 , x ) ] , #NEWLINE# #UNINDENT# } ) #NEWLINE# raises ( NonElementaryIntegralException , ( lambda : normal_denom ( Poly ( 1 , x ) , Poly ( 1 , x ) , Poly ( 1 , x ) , Poly ( x , x ) , DE ) ) ) #NEWLINE# ( fa , fd ) = ( Poly ( ( ( t ** 2 ) + 1 ) , t ) , Poly ( 1 , t ) ) #NEWLINE# ( ga , gd ) = ( Poly ( 1 , t ) , Poly ( ( t ** 2 ) , t ) ) #NEWLINE# DE = DifferentialExtension ( extension = { #NEWLINE# #INDENT# 'D' : [ Poly ( 1 , x ) , Poly ( ( ( t ** 2 ) + 1 ) , t ) ] , #NEWLINE# #UNINDENT# } ) #NEWLINE# assert ( normal_denom ( fa , fd , ga , gd , DE ) == ( Poly ( t , t ) , ( Poly ( ( ( ( ( t ** 3 ) - ( t ** 2 ) ) + t ) - 1 ) , t ) , Poly ( 1 , t ) ) , ( Poly ( 1 , t ) , Poly ( 1 , t ) ) , Poly ( t , t ) ) )