@Jacobegerup Your attachment doesn't contain input exactly like your image. It's not clear which output(s) in it are presenting you with a similar issue. You should explain where the problem is, explicitly.

@Jacobegerup Could you upload and attach your Document (.mw file), using the green up-arrow in the Mapleprimes response editor's toolbar?

I don't know what version of Maple you're using.

But what happens if you use the name mass instead of m, for the assignment of 142g ?

@salim-barzani Sure you could do that in several different ways. For example,

plots:-display(seq(plot(last[2], l = -1 .. 1,
                        legend = typeset(`#mn("a");`=a)),
                   a = 0.3 .. 1, 0.1),
               legendstyle = [location = right]);

or just,

plot([seq(last[2], a = [0.3, 0.5, 1.0])], l = -1 .. 1);

etc. You can, of course, change those to other values of a. Or you could use other colors.

If you'd like those curves as contours on the 3D surface then that's possible in several ways. You haven't specifically asked about that, though.

If those aren't exactly what you're after then you're out of luck with me since you've now posted five times in this Question thread without properly and exactly explaining what you're after.

@Carl Love The coding style is mmcdara's, cribbed and tried on a new example...

@salim-barzani You last attempt tries contourplot mistakenly with ranges for x and y. But your expression has parameters l and a, not x and y. So, naturally, that doesn't work.

I don't understand your wording. You wrote this as one long sentence, with no punctuation: "this is good but i need to complet the work i add some part in below which beside the shape of 3D and 2D which must the result be smaller and smaller in three different of four level of changing parameter the 2D and 3D both importan like the first maple file   and yours which is so clear  there is anyway for for do something like that". It's incomprehensible to me.

Here are a half-dozen or so variants on plotting your bivariate expression.

g1_ac-explor_ac2ooo.mw

You should explain clearly what it is you want. You have not come even close to doing so, yet. I will not repond further unless you explain your goal and requirements much better.

A few of the included plots:

@salim-barzani I find your description of what you want to be confusing and incomplete. Please try and explain exactly what you want, with more explcit nouns and fewer vague pronouns.

You can get a 3D plot from your second example,

plot3d(G(l, a, 1, 1, 1, 1, 5), l=-2..2, a=0..1.5,
       colorscheme="turbo",
       style=surface, lightmodel=Light1,
       orientation=[-110,75,0]);

Or you could Explore the 2D plot, wrt parameter `a`.

Or you could get 2D contourplot or densityplot. And so on.

Bgraph_acc.mw

In your middle example (that shows error messages) your procedure G is malformed.

The mistake is that you have subscripted names entered using two different 2D Input syntax mechanisms. The procedure's parameters are defined as names with double-underscore (eg. alpha__3) while the body of the procedure has indexed names (eg. alpha[3]).

If you convert the 2D Input to 1D Input then the synax muddle becomes more visually apparent. Yes, there is a line of code to do assignments of one form to another, but that's fragile in the face of the way you do operator construction.

Bgraph1_ac.mw

You had kind of problem a few times before. I think that you'd be better off using just one of those two forms, consistently, throughout your worksheet, than mixing them. Especially if you use 2D Input syntax for operator definitions, etc.

To use Explore you just have to combine steps 4 and 5. The separation into steps 1-5 is not all necessary.

I upvote.

Are you considering putting this version onto the Application Center and/or the Maple Cloud?

@dharr 

[edit] I have made a correction in my code below, and have edited this Reply.

Using Maple 2024.2, the univariate numeric integration method=_d01amc appears to be able to handle directly the semi-finite integral (with your evalc'd form) faster than does method=_d01akc for the oscillatory integrand transformed to a finite range.

The plot3d with grid=[40,40] can then be done in  less than 3.3sec on my machine, in contrast to 9.6sec with that _d01akc.

I have not looked to see whether there was anything faster still.

Download Fokas_acc.mw

@sand15 Here's another variant, for Maple 2015.

I showed two uses with value, just to illustrate that while only one of them might display the nicer way for output, both could be used to compute with in some ways.

Yes, I realize that some variants produce that particular error (evaluation of an rtable indexed by a name rather than integer value). That's indeed why I looked again for Maple2015-specific workarounds...

Download Product_error_ac_2015_2.mw

Do you mean like with a legend, or caption, or inlined textplot, or...?

How do you want to "link" them?

Did it work for you?

@Muhammad Usman 

Download Usman_contplot3.mw

@Muhammad Usman Why didn't you mention all your requirements when you first posted the Question?