Dear all
If you have the following two maple package, kindly share it with me.
wkptest: For the symbolic computation of Painleve test for nonlinear PDEs.
ONEOptimal:-  Maple Package for obtaining Optimal System

Thanking you in advance
Debendra

How can I make a legend for this plot?

plot([17.85*(2.65*t^2 + 1 + 3*t)^2/t^2, 17.85*(2.65*t^2 - 47 + 3*t)^2/t^2, 17.85*(2.65*t^2 - 97 + 3*t)^2/t^2, 17.85*(2.65*t^2 - 147 + 3*t)^2/t^2, 17.85*(2.65*t^2 - 197 + 3*t)^2/t^2], t = 0 .. 25, labels = ['t', 'x'], labelfont = [Times, 12])

 I got the following figure.

Now I want to put a legend box in this plot for each color line. How can I do it?
 

Suppose I have an expression like the following:

h(t, x) :=  (m*t^2 + 6*t - 2*x)^2/(36*g*t^2)

I want to calculate multiple values of t in a single expression. Say, for t:= 1,5,10,15,20. I want to evaluate h(t,x) in vector form. so that I can extract for each t. How to write it?

How do we define a matrix or a vector of the partial differential operator?
 

Suppose I have, a row vector  R = [R₁ , R_{2  ,} R₃ ] and a column vector C =  [C₁ , C_{2  ,} C₃ ]. I need a multiplication like as follows

RC = [ R₁ [C], R₂[C], R₃ [C]  ], so that RC will be a 3/3 matrix.