here is a matrix

A := Matrix(2, 2, {(1, 1) = ((1/2)*I)*E[0]*d[ba]/`ℏ`, (1, 2) = -((1/2)*I)*delta, (2, 1) = ((1/2)*I)*delta, (2, 2) = ((1/2)*I)*E[0]*d[ab]/`ℏ`})

it is clear that there is a common coeff with each element by ((1/2)*I,i just want to show the matrix with traditional way.

#A:="((1/2)*I*Matrix(2, 2, {(1, 1) = -E[0]*d[ba]/`ℏ`, (1, 2) = delta, (2, 1) = -delta, (2, 2) = -E[0]*d[ab]/`ℏ`})"

A1:=((1/2)*I;

here is an equation

eq := I*`ℏ`*(sum((diff(c[n](t), t))*f[n](r)*exp(-I*omega[n]*t), n = l .. k)) = (1/2)*E[0](e_.r_)*e*(sum(c[n](t)*f[n](r)*omega[n]*(exp(I*t*(-omega[n]+Omega))+exp(-I*t*(omega[n]+Omega))), n = l .. k))

req := Int(conjugate(f[m](r))*rhs(eq), r)

leq := Int(conjugate(f[m](r))*lhs(eq), r)

sup:=Int(conjugate(f[m](r))*f[n](r), r) = delta[m, n]

ans:=applyrule(sup,leq=req)

how can get the result as

ans := I*`ℏ`*(diff(c[m...

I have got an equation of Sum.i want to split the eqution three parts of Sum,how can do?

ans:=sum(I*`ℏ`*(diff(c[n](t), t))*phi(r)*exp(-I*omega[n]*t)-(1/2)*E[0](e_.r_)*e*c[n](t)*phi(r)*omega[n]*exp(-I*omega[n]*t+I*Omega*t)-(1/2)*E[0](e_.r_)*e*c[n](t)*phi(r)*omega[n]*exp(-I*omega[n]*t-I*Omega*t), n = l .. k)

i want to  result like this

ans1:=sum(I*'ℏ'*(diff(c[n](t), t))*phi(r)*exp(-I*omega[n]*t), n = l .. k) = sum((1/2)*E[0...

restart; with(Physics); Setup(mathematicalnotation = true):

am := Annihilation(N); ap := Creation(N)

Setup(op = {b[1], b[2], b[3]})

alg := seq(%Commutator(ap, b[i]) = 0, i = 1 .. 3), seq(%Commutator(am, b[i]) = 0, i = 1 .. 3), %Commutator(b[1], b[2]) = b[3], %Commutator(b[1], b[3]) = -Physics[`*`](2, b[1]), %Commutator(b[2], b[3]) = Physics[`*`](2, b[2])

q := %Commutator(Physics[`.`](ap+am, b[2]), b[1])

p := expand(q)

ans := sort(p, [a, b])

p := (g+Delta*g*Zeta/omega-g*Zeta)*a[1]*b[2]+(g+Delta*g*Zeta/omega-g*Zeta)*a[2]*b[1]+(-Delta*g*Zeta/omega+g-g*Zeta)*a[2]*b[2]+(-Delta*g*Zeta/omega+g-g*Zeta)*a[1]*b[1]

q := factor(add(op(i, p), i = 1 .. 2))

q := g*(omega+Delta*Zeta-Zeta*omega)*(a[1]*b[2]+a[2]*b[1])/omega

here,throuth factor ,simplify the expression .yet,if i want the expression like this:

(g+Delta*g*Zeta/omega-g*Zeta))*(a[1]*b[2]+a[2]*b[1])

how do i?

thanks