@acer I was rewriting my answer while you replied.

@Dark Energy Works for me...

@amrramadaneg Works fine...

with(plots):
V := view = [0..2, 0..2];
display(arrow([1,0]), arrow([0,1]), V);
display( plottools[transform]((x,y) -> [x+1, y+1])(%), V);

This is an incomplete specification of a problem. No one can help you with it until you post f & g.

Post the code text or the original file please.

@ajfriedlan See the difference in :
(1,2,3)[1], (1)[1];
[1,2,3][1], [1][1];

Your code uses the first version, it should be using the second. That is:
 Rootz := [solve](epsilon__k(q, J__a[l])-epsilon__k(kappa, J__a[l])+.16*e = 0);

You can always compose plots in the end, so I suggest you work on them sparately first.

@WJRAGG Simply evaluate the following list: [K, P, Q, ap]; and print the output here.

The parameters : {K, P, Q, ap} are missing from your code.

@Carl Love

g := subs(F = unapply(f,a,b,x), (a,b) -> (x-> F(a, b, x)));

@Carl Love This g does not do what was intended.

Also I have a Notification "Carl Love mentioned you in a Answer". Why does it say that?

@666 jvbasha or plots[pointplot3d]([seq(seq([bix[i],prx[j],X1[i,j]],i=1..NN),j=1..NN1)]);

@666 jvbasha try: print(eval(X1));

@Carl Love I dont think using indets will ever match subsindets code. I expect subsindets recurses using type & hastype.

Another difference:

e:= f(g(y), f(f(f(x), g(y)), f(x)), g(y)):
i:= 0:  SubsIndets(e, function, h);
i:= 0:  subsindets(e, function, h);

@ecterrab That was my mistake. I forgot type testing was done to the initial value not to the resulting value. "You find that which you seek".