var('x,y') right=integrate(integrate(sin(x^2),y,0,x),x,0,1) wrong=integrate(integrate(sin(x^2),x,y,1),y,0,1) eq = real(wrong)==right; eq
bool(eq)
eq = SR(2) == SR(2); eq
type(eq)
parent(eq)