You can get fixed font by editing in source mode (press "Source" button) and enclose text between <pre> </pre> tags.

You mean perhaps an ascii source and correct it with an editor?

One posibility is using lists of pairs [Re,Im] for each root, like this in Maple 13:

with(plottools):with(plots):
ball := proc(x,y) plots[pointplot]([[x,y]],color=blue,
symbol=solidcircle,symbolsize=30) end proc:
[solve(z^3-z+A,z)]:evalc~(%):
s2:=[Re,Im]~(%) assuming A::real:
a1:=animate(ball, s2[1], A=-2..2, scaling=constrained, frames=20):
a2:=animate(ball, s2[2], A=-2..2, scaling=constrained, frames=20):
a3:=animate(ball, s2[3], A=-2..2, scaling=constrained, frames=20):
display([a1,a2,a3]);

My interpretation of the output from these pieces

true,...,[{p = 0, q = 0, r = 0}],[[[e1, e2] = e5-e6, 
[e1, e3] = e4-e6, [e2, e3] = -e5+e6]]

is that a 3-dimensional Lie algebra generated by {e1,e2,e3} is meant. But it not clear either for me.

E.g.:

sys:={diff(x(t),t)=x(t)^2+y(t)+z(t)+w(t), 
diff(y(t),t)=x(t)+y(t)^2+z(t)+w(t),
diff(z(t),t)=x(t)+y(t)+z(t)^2+w(t), 
diff(w(t),t)=x(t)+y(t)+z(t)+w(t)^2}:
IC:={x(0)=0,y(0)=0,z(0)=1,w(0)=0}:
sol:=dsolve(sys union IC,numeric):
plots:-odeplot(sol,[x(t),y(t)],t=0..0.5);

Look at ?plots,odeplot

It sounds something like:

s:=dsolve({diff(r(t),t$2)+r(t),r(0)=1,D(r)(0)=0},numeric,output=operator):
dr:=rhs(s[3]):
(evalf@Int)(dr(t)^2,t=0..2*Pi);

                             3.141594697

diff(y(x[1],x[2]),x[1]);
                              d                 
                            ------ y(x[1], x[2])
                             dx[1]              
   
diff(sin(x[1])*exp(x[2]),x[1]);

                             cos(x[1]) exp(x[2])

Some closing brackets are missing in the system written above. Check them. Also if 'D' is meant as a constant, better use something else, D1 say.

You can differentiate the equation, and then solve this differential consequence:

d:=diff(pdeTest5P1,`ϕ`):
dsolve(d);
           _FF3(ϕ) = _C1 sin(ϕ) + _C2 cos(ϕ)

Just check then that this solution satisfies the initial equation. Better post the code text in Maple syntax.

If v means the absolute value of the velocity, hence v>=0, and m>0, etc, you have to use these pieces of information:

solve({E = (1/2)*m*v^2, v>=0},{v}) assuming positive;
                               1/2      1/2
                              2    (m E)
                         {v = -------------}
                                    m

In Maple 13:

A+~100;

For observing the communication between interface and kernel you need to capture the loopback device.

h:=proc()
INTERFACE_DIAG(Hola);
end proc:

h();

Hola

One posibility is a series expansion of the solution:  

Order := 5:
dsolve(sys union ics, {x(t), y(t)},series);

           P   2   C P   3
  {x(t) = --- t  - ---- t  + 1/24 P (-3 Delta M - KMIN M
          2 M         2
                   6 M

                                                     2
         - M KMIN exp(20 Pi) - M Delta exp(20 Pi) + C

            2               /   3                    4      5
         + C  exp(20 Pi))  /  (M  (1 + exp(20 Pi))) t  + O(t ),
                          /

                       P Delta exp(20 Pi)   4      5
        y(t) =  - 5/6 -------------------- t  + O(t )}
                       2                 2
                      M  (1 + exp(20 Pi))