The global name gamma has a special meaning in Maple, and refers to a particular real constant. (That's why you cannot assign a new value to it. The name is protected.)

The error message from the implicitplot command is due to your use of the global name of the constant gamma on the left hand side of R.

The simplest course of action in your Maple 11 is to use another name instead.

For your first question, it seems to me you are asking about the functionality provided by the unapply command.

And for your second question, a DataFrame might possibly serve.

See a link to the attachment below.

unapply_df.mw

The first step can be done using the rationalize command.

The second step can be done by rewriting the given substitution equation.

Download Resonant_xformer_ac.mw

Here's another variant, using InertForm with the inert operators for the divisions and additions, and with the active operator for multiplications. The active multiplications are implicitly displayed in output (ie. the multiplication symbol suppressed).

Download SérieAléaEqFractionnaires1_ac.mw

Your document had the same operations done multiple times. You don't have to load the InertForm package (utilizing the with command). You certainly don't have to call with multiple times.

You could turn the example generator into a procedure, to make it all cleaner.

You could also use Embedded Components, with each new example replaing the old one. (It could be as simple, or as fancy, as you choose.)

Your document also made me think of the Quiz package.

Download sqrt_evalc.mw

How do you feel about t=0 ?

The odeplot command makes it easy to specify what goes on which axis.

Download odeplot_v.mw

A small modification of your code gets you the following. Is it something like your goal?

You could easily make it a procedure.

Download adjoint.mw

Download rad_combine.mw

I have not utilized the cosine property cos(-z)=cos(z) however.

You can represent this operation with a call to %.  as an inert noncommutative multiplicaion.

The dot symbol renders in gray by default, but you can also use InertForm:-Display to render it in black.

The value command is one mechanism (of several) that turns it into an active computation.

Download inertmatmul.mw

[edit] Changed to an infix call, see comments below.

simplify(eval([y,z],aa[1])-eval([y,z],bb[2]))
   assuming alpha[1,8]>0;

                [0, 0]

radnormal(eval([y,z],aa[1])-eval([y,z],bb[2]))
   assuming alpha[1,8]>0;

                [0, 0]

I do not recommend in general your original approach of comparing a pair of simplify results using evalb, ie. testing whether they are the identical same structure. In general there may not be a definitive normal form, and so it is an unnecessarily weakened test -- despite the fact that it succeeds in this example involving radicals.

Under the assumption used above,

# I don't much like this first one, even if
# it works for your polynomial solutions.
evalb(simplify(aa[1,1])=simplify(bb[2,1]))
   assuming alpha[1,8]>0;

                true

# I'd prefer a check more like this.
is(simplify(aa[1,1]-bb[2,1]))
   assuming alpha[1,8]>0;

                true

# Alternatively,
is(simplify(eval(y,aa[1])-eval(y,bb[2]))=0)
   assuming alpha[1,8]>0;

                true

You didn't state what you wanted the titles to be.

>  restart; >  dim:=3: flat:=(j - 1 + dim*(i - 1))+1; pltopts:=color=blue,size=[250,250],tickmarks=[0,0]: s:=[seq(     [seq(plot(diff(x^9,[x$(flat)]),pltopts,               title=typeset(InertForm:-Display(%diff(x^9,[x$(flat)]),                                                inert=false)                             = diff(x^9,[x$(flat)])),               titlefont=["Times",10]        )     ,j=1..dim)], i=1..dim)]: A:=Array(1..dim,1..dim,s): >

Download QuestionTabulation_ac.mw

Put both in one more call to plots:-display.

You have mistmatched brackets in the expression assigned to wres.

You forgot the colon in the colon-minus ( := ) in the line that seems to define z3. As you have it, that line is just an equation, rather than as assignment.

You used lowercase pi when perhaps you intended the popular constant Pi which is capitalized in Maple.

From the Help page (Maple 18) for the Compiler,
   "Arrays can only appear as parameters of a procedure that is to be compiled."

In other words, you can construct your Array outside of the procedure, and pass it in as an argument. You shouldn't be able to assign an Array to a local variable within the compiled procedure. (The same should go for Vector and Matrix.)