# T is the 3x+1 function
restart: n := 87: # Change n!
 




 T := proc(n:posint)
 if type(n,odd) then (3*n+1)/2 else n/2 fi
 end:

 
 print(`T(` .n. `)` = T(n));