using Symata  # load the package and enter Symata mode

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

OutputStyle(InputForm)
Out(2)

OutputStyle(UnicodeForm)
Out(2)

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

OutputStyle(InputForm)
Out(5)

OutputStyle(JupyterForm)
Out(5)

OutputStyle(InputForm), 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(JupyterForm);

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

g(x_) := Exp(-x)

Out(11)

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

CompactOutput(True), Out(10)

OutputStyle(JupyterForm), (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

InputForm( a< b < c/d )

N(Pi)

N(Pi) * 10^10

FloatFormat(Short);

N(Pi)

N(Pi) * 10^10

FloatFormat(3);

N(Pi)

N(Pi) * 10^10

Precision(N(Pi))

N(Pi,50)

BigFloatFormat(Short);

N(Pi,50)

Precision(N(Pi,50))

FloatFormat(Full);

N(Pi)

BigFloatFormat(Full);

N(Pi,50)

OutputStyle(InputForm);

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

OutputStyle(UnicodeForm);

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

OutputStyle(JupyterForm);

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

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

Cos([Pi, π])

? OutputStyle

VersionInfo()

Now()