I (only) use the classical interface (*.mws) and without too much pain
(which only needs some guessing which 'paragraph' stands for what)
you can create a style sheet through the menu bar 'styles' and then
save it as default.

Here you can choose colors, sizes and fonts as you like them. And
I use ~ dark red and ~ dark blue (for better contrast) and kick out
all the italics for better ergonomy.

It works similar like a normal.dot for Win Word and is very usefull.

First: even I often swear at a system - for example Maple - I meanwhile hesitate to use 'bug', simply because too often I finally recognized improper usage.

In your case you actually say " (r=0 or rr=0) and a*p*d-rr = 0; give me d "

And Maple refuses an answer - in a very unfriendly way.
Difficult to understand and to accept, but quite often the system assumes variables to be not zero, for example if asking to solve a*x=b for x.

p1 := 1/2*(p[L]+p[R]);
#  ^


You missed the double point to define p1

PS: where I would use 14/10 instead of 1.4 for 'gam' and then
f[L] = 0, df[L] = sqrt(5/7) exactly

?unapply

  [solve(a*x^2+b*x+c=0, x)]; # a list
  %[1];                      # take the first entry 
  f:=unapply(%, x);          # and make it a function of x

Note that this silently assumes a <> 0.

Without being an admin: that sound like a problem of the printer driver (may be related to Java version), since you say it works with an older model

try ?degree or ?lcoeff depending on your needs

As far as I see: 'out' does not have a value and thus the for-loop does not work.

What do you want to process? The elements of your S?

PS: Please do not duplicate threads, there is not need to post the same question twice,
http://www.mapleprimes.com/questions/100067-Code-Generation-Query

xi=ln(x); 
newCoord:=isolate(%, x);

                       newCoord := x = exp(xi)

ln(x); eval(%, newCoord); 
plot(%, xi=1 ..3);

Also see ?plots[logplot]

Sorry, I was quite absentminded.

It happens as well in a concurrent test version. And there:

J := 2*sin(Pi*k)/k;

J assuming k::integer;
                                  0

limit(J, k=0);
                                 2 Pi

Which is: no care is taken for the in-definte case.

Moreover:

sol:=int(exp(I*k*x),x);

                                -I exp(k x I)
                         sol := -------------
                                      k
discont(sol, x);
                                  {}

int(exp(I*k*x),x) assuming k::integer;

                            -I exp(k x I)
                            -------------
                                  k

So Maple cares only for the 'generic' case.

Not sure, whether that is the (?) reason.

Hint: what is the root of (x - 7) ?

Do you need to care for leading coefficients, looking at 2*(x-7)?

Do you need to care for degrees, looking at x*(x-7)?

 ?                     

temp:=phi(Lambda(tau, 0), tau, 0)+
((D[1](phi))(Lambda(tau, 0), tau, 0)*(D[2](Lambda))(tau, 0)+
(D[3](phi))(Lambda(tau, 0), tau, 0))*q;

eval(temp, Lambda = '(x,y) -> c[y](x)');

  phi(c[0](tau), tau, 0) + (
        D[1](phi)(c[0](tau), tau, 0) D[2]((x, y) -> c[y](x))(tau, 0)
         + D[3](phi)(c[0](tau), tau, 0)) q

eval(temp, Lambda = '(x,y) -> c[0](x)');

       phi(c[0](tau), tau, 0) + D[3](phi)(c[0](tau), tau, 0) q

As you see, that is 'sensitive' to the interpretation of input
and you may want to think about using indexed variables
(which often are used for 'pretty appearance' only).

# your sheet gives the following integrand E
E:=exp(-B*((1/8)*f/Pi^2+(1/2)*f*x^2));

# you want Int(Int(E, f),x) -  interchange variables
Int(Int(E, x),f); value(%); subs(int=Int, %);

# change variables, 0< B, 0 < f assumed for that
Change(%, B*f/4/Pi^2 = y^2, y);

# beautify it
combine(%) assuming 0 < y;
expand(%); lprint(%);

2/B*Pi^(3/2)*2^(1/2)*Int(erf(1/2*y*Pi*x*8^(1/2))*exp(-1/2*y^2),y);

               3/2  1/2   /             1/2          2
           2 Pi    2     |      y Pi x 8            y
           ------------  |  erf(-----------) exp(- ----) dy
                B        |           2              2
                        /

Now 1 +erf ~ cumulative normal distribution cdfN and exp(-1/2*y^2)
is like its derivative, pdfN.

So this is Int( (cdfN - 1) * pdfN ) = Int( cdfN * pdfN ) - Int(pdfN).

The last on the right hand side is cdfN. And for the first one has
to know, that this is a bivariate normal (see Abramowitz & Stegun?):

Int( pdfN(xi)* cdfN((y-rho*xi)/sqrt(1-rho^2)),xi=-infinity..x) is a
way to write the bivariate in x and y with correlation rho.

What remains is to determine the 'settings' to pass from your task
to the bivariate.

Hm ... at least any ranges for your variables and parameters B, x, f?

For me M 9.52 on Windows 7 worked in a test (not very deep and classic worksheets), essentially.

The only problem: exporting things (what kind ever) crashed the application