using Symata; isymata()   # load the package and enter Symata mode

Integrate(f(x+I),[x,0,Infinity])  # input is interpreted as Symata code

OutputStyle(Plain), Out(2)

OutputStyle(Unicode), Out(2)

Assume(a, Positive), Integrate(E^(-x) * x^(a-1), [x,0,Infinity])

OutputStyle(Plain), Out(5)

OutputStyle(IJulia), Out(5)

OutputStyle(Plain),Expand((x+y)^3)

FullForm(Out(8))  # Internal form of the previous output

Plus(Power(x,3),Times(3,Power(x,2),y),Times(3,x,Power(y,2)),Power(y,3)) # This is also valid input

OutputStyle(IJulia);

Integrate(g(x), [x,0,Infinity])

( g(x_) := Exp(-x), Out(12) )

(OutputStyle(Plain), CompactOutput(False), Out(10))

CompactOutput(True), Out(10)

OutputStyle(IJulia), (1/2 +  a^b)/(x+y) 

Sum(g(i,j), [i,0,Infinity], [j,0,Infinity]) + Sum(h(i,j), [i,0,Infinity], [j,1,n])

Integrate(g(x,y), [x,0,1], [y,0,1])

a < b < c/d

OutputStyle(Plain);

[π + 𝕖 + 𝕚 + a, Pi + E + I + a]   # \pi + \Bbbe + \Bbbi

OutputStyle(Unicode);

[π + 𝕖 + 𝕚 + a, Pi + E + I + a]

OutputStyle(IJulia);

[π + 𝕖 + 𝕚 + a, Pi + E + I + a]

[Pi == π, E == 𝕖, I == 𝕚, EulerGamma == γ, Gamma == Γ]

Cos([Pi, π])

? OutputStyle