is performed by reflection through a plane which contains these axes.

Your problem statement is confusing by the fact that (in general) most people thinkk of cartesian axes as x, y, x in a 'sequence'. You however have labelled these axes as Z, Y, X - so pretty much reverse order. So when you say you want " the X and Y-axes reversed ", I'm not completely clear whether you want the first and second axes reversed (ie x and y) or the second and third axes reversed (ie the ones you have labelled X and Y)

However the principle remains the same: the attache shows your original plot, then shows it reflected through all axes planes - ie firstAxis/secondAxis, firtAxis/thirdAxis, secondAxis/thirdAxis. One of them has to be correct. You just have to pick which one

I tried uploading the following with all its graphical output, but (for reasons unknown) this continually failed. Hence the following is 'input code' only

Download reverseAxes.mw

1. solve() in the attached was taking too long to run so I eventually gave up on it. If yu are more patient than me, then it may be worth 'uncommenting' to see if anything meaningful happens
2. the basic fsolve(), with range contstraints,  produces one solution, with very respectable residual errors, and is very quick
3. the SolveEquations() command from the DIrectSearch package, with the option AllSolutions=true produces six solutions, in pretty acceptable time. The first four of these have very respectable residuals - and I'd take all of them as 'solutions'. Both of the solutions generated in Kitonum's earlier answer are in this set of four. The residuals on the final two solutions are rather high - so I'd treat these with care. However, the residuals in these final two solutions probably ought to be compared with the magintudes of some of the terms in your original equations (as in O(10^13), so a residual of ~8, isn't too shabby. I wouldn't discount these final two solutions without some further investigation.

The DirectSearch package is not a part of Maple but is an add-on available from the Maple Applications Center. Please note that until/unless you have the DirectSearch package loaded, the final execution group in the attached worksheet will not execute:-(
 

Download numSols.mw

Consider the following

#
# When Maple 'parses' an expression, it does so
# in a two-step process, known as simplification
# and evaluation. So
#
  x+2*x;
#
# will always return 3*x, because this is the
# 'simplification' process. Using uneval quotes
# which only affects the 'evaluation' stage will
# not help. Hence
#
  'x+2*x';
#
# will also return 3*x. Now you can control the
# evaluation process, so
#
  y:=2*x:
  x+y;
#
# will return 3*x, but
#
  y:=2*x:
  'x+y';
#
# will not - because y is not evaluated in the
# expression 'x+y'
#
# In your case of latex(someExpression), your
# problem is that Maple will simplify and
# evaluate the argument 'someExpression', before
# passing it to the latex() command. You can
# prevent the evaluation stage by using uneval
# quotes, but you cannot prevent the simplification
# stage, so
#
  latex(x+2*x);
#
# will give the same answer as
#
  latex('x+2*x');
#
# The following will give the same answer, since
# the evaluation stage is not prevented
#
  y:=2*x;
  latex(x+y);
#
# However the evaluation stage could be prevented with
#
  y:=2*x;
  latex('x+y');
#
# The fundamental problem which you have is that
# for the expression you are supplying, Maple will
# perform its basic simplification operation, which
# you cannot prevent. So
#
  latex( 5*x^3+7*x^2+2*x^3);
#
# and
#
  latex( '5*x^3+7*x^2+2*x^3');
#
# will produce identical results. The only way I can
# think to get around this is to perform the latex
# export as a string - and then sutably 'parse' the
# resulting string in the target application, as in
#
  latex(" 5*x^3+7*x^2+2*x^3")

Check the attached

Download solveEQ.mw

Having defined a 3-D surface plot (or a 2-D curve), you can use the plottools:-project() command to 'project' the plot on to any specified plane (or a 2d plot on to any line).

You can then use plots:-display() to combine the original plot and its projections in a single figure. See the "toy" examples in the code below

  restart;
  with(plots):
  with(plottools):
  p1:= plot3d( [ x+y+4, 25-x^2-y^2],
                x=-2..2, y=x^2..2-x^2
             ):
  pa:= shadebetween( x+y+4, 25-x^2-y^2,x=0..1,
                     y=x^2..2-x^2,
                     showboundary = false,
                     changefill = [color = khaki, transparency = 0],
                     axes=normal,
                     lightmodel=light1
                   );
#
# Use plottools:-project(), to project on to any desired
# plane
#
  p2:= project(p1, [[0,0,0], [1,0,0], [0,1,0]]); #xy-plane
  p3:= project(p1, [[0,0,0], [1,0,0], [0,0,1]]); #xz-plane
  p4:= project(p1, [[0,0,0], [0,1,0], [0,0,1]]); #yz-plane
#
# Combine the original surface plot and the projected plots
# in the same figure
#
  display([p1,p2,p3,p4]);
#
# Combine the 'shadebetween' plot and the projected plots
# in the same figure
#
  display([pa,p2,p3,p4]);

If the calculation which you perform depends on a simple loop index, then you can use an initialiser function with the Vector definition. So for a "toy" example

a:=Vector(10, i->i^2);

will load a 10-element Vector with the 'squares' of the indices. Obviously much more complicated functions can be achieved. The one thing you can't do this way is to have a particular element depend on the result of a previous element calculation - but a 'threaded' calclulation would also fail in this circumstance. With some care you could use the above with threads.

I am wondering how we can edit a diagram plotted in maple. for example in editdiagram.mw I want name tips of diagram as "a", "b",.... and write them besides the tips on diagram.

If you could be bothered to provide the procedure/function/whatever which generated this diagram then I am pretty sure I could figure out assorted extremal points. I could then use textplot() to insert any text I want at these extremal points. Then use display() to overlay the original plot and this additional text. Such a pity that you did not upload any way of generating your plot so that I could demonstrate!

############################################################################

Another question is how we can make color different parts of a single diagram by various colors? for example in above diagram we want make blue vertical line, make green horizontal line and make black the curve line?

Could be tricky!!  depends very much on how the procedure/function/whatever which determines the plot is defined - which, of course, you are not telling anyone. Basically if a single expression defines the complete curve, then you will not be able to colour parts of it separately. If it is defined using some kind of 'piecewise' process, then individual colouring of sections may be possible

##########################################################################

I suggest you provide sufficient information for your problem to be addressed

This looks very much like an addition to

https://www.mapleprimes.com/questions/222602-Numeric-Evaluation-With-Units#comment241089

It is considered bad practice here to post a 'new' problem, when your issue is actually an extension of a previously posted problem.

However the attached file should resolve your issue. Regrettably (again), I have used Maple's 1-D input mode, which you didn't like last time (and won't like this time). I use it, because it works. You, of course, are free to use 2-D input mode??

Download unitProb3.mw

If I simplify everything by using 1-D input (no smart plots, palette entry etct, etc), so that I just type exactly what I want to be calculated, then everything seems to work as expected, as in the attached

Download unitProb.mw

Just working from the code you have posted and ignoring the irrelevant(?) comments such as

I use some commands like the evalf[10]

But they don't work for the definite double integral in attached file.

It leads to "Error, (in type/algfun) too many levels of recursion"

You are using the function Sin() in a couple of places. I assume that this should be sin() - note the case. If I change your Sin() declarations to sin(), then the attached works (or at least provides a value with no error - only you know whether this value is vaguely plausible!)

 

Download doInt.mw

Since you don't provide expressions for the actual plots which you want, I can't guarantee to fix them :-(

If you want a quick lesson on how to plot more or less tickmarklabel you want, at more or less any tickvalue you want, in more or less any font you want, then examine the following carefully

   restart;
#
# Produce a plot where every tickmark (on x and y)
# is specified explicitly
#
   plot( x^2,
            x        = 0..4,
            axesfont  = [times, bold, 20],
            axis[1]     = [tickmarks=[0.5="0.500", 1.5="1.50", 3="3", 3.5="A"]],
            axis[2]     = [tickmarks=[3="3.00", 7.5="B", 15="15.000"]]
         );
#
# Now do the above in a more-or-less programmatic way
# where the 'format' of the x-axis tickmarks is being
# controlled. (Obviously I could do something similar
# with the y-axis tickmarks)
#
   xvals:= j$j=0.5..3.5:
   xreps:= [seq( xvals[j]=sprintf( "%4.3f", xvals[j]), j=1..4)]:
   plot( x^2,
            x        = 0..4,
            axesfont = [times, bold, 20],
            axis[1]    = [tickmarks=[xreps]],
            axis[2]    = [tickmarks=[3="3.00", 7.5="B", 15="15.000"]]
         );

If the above does not answer you question, then please provide equations for the 'curves' which you want (in a plottable form, so everything defined), together with a lists of x-axis and y-axis tickmarks which you want, in the formats you want

I duplicated your Error message - couldn't see any obvious reason, apart from the use of 'gamma', which, in Maple, is Euler's constant.

Unless you really, really know what you are doing, you should never, ever 'unprotect' Maple's known constants, so I decided just to change the name of this function to h(). Naturally one cannot use find-and-replace in 2D Math input, so I converted you worksheet to 1-D, input and substituted h() for gamma() everywhere.

Hey presto everything now works - See the attached

odeProb.mw

If you skip the second argument, as in

singular( 1/((x-b)*(x-a)));

then Maple correctly(!) returns the two singular points defined by

{a = a, b = x, x = x}, {a = x, b = b, x = x}

Why this doesn't happen when you define 'x' to be the variable is a mystery to me :-(

Checking the 'help' for surfdata(), you are using the 'second' calling sequence, ie

surfdata(g, r1, r2, opts)

and paraphrasing the 'help' for this calling sequence.

In this case, the x-values and y-values are assumed to correspond to those of a uniform grid over the region defined by ranges r1 and r2.

Now this means that your array of data can be arbitrarily scaled/translated/stretched, simply by appropriate selection of the ranges r1, r2.

Note that the help also clearly states

If r1 or r2 is not provided, then the range 0..1 is assumed.

Once you understand the above, you can 'combine' two surfdata() plots (more or less as Acer has suggested) using

p1:=(f,  r1, r2):
p2:=(g, r3 ,r4):
plots:-display([p1,p2]);

max[index]()

as in the 'toy' example

with(LinearAlgebra):
p:=RandomMatrix(10,10);
mval:=[max[index](p)];