@henrylyl I don't know about the units conversion in plots. The answer may be on the help page ?plot,units . There is a lot of information there, and I don't fully understand it yet. I need to do some experiments.

Regarding your other question about displaying the array of plots, it can be done like this:

display(< display~(G), display~(< pG[] >)^%T >);

Do you have a numeric value for Pr?

Do you have any initial or boundary conditions? You'll need eight.

@bastorer 

It should work the way that you describe. But if Array2:= Copy(Array1) does not give an error message, then you can't be using ArrayTools:-Copy. Can you duplicate the situation?

@Preben Alsholm That works!

@bastorer

Array2:= Copy(Array1);

is not valid syntax for ArrayTools:-Copy. You must have another package loaded that has a Copy command.

@henrylyl The file link works now.

@Preben Alsholm It doesn't work recursively. Example:

restart:
Array1:= Array([Array([1, 2]), Array([3, 4])]):
Array2:= copy(Array1):
Array1[1][1]:= 11:
Array1;
Array2;
              

That was asked and answered 10 days ago, probably by a classmate of yours. See here.

@Markiyan Hirnyk 

I see your point. I wonder why the presence of r with no assumptions makes Maple think that it knows the integral.

@henrylyl File didn't attach. Try uploading again please.

@Markiyan Hirnyk Just plain integrate it:

int(BesselJ(2, r*k)*BesselJ(1, 1500*k), k= 0..infinity);

                                                0

@Markiyan Hirnyk Thanks. I got all files from Christopher's attached zip file.

Please supply (as an attached worksheet) the complete code that you used to read in the data. And please also attach the data files.

The worksheet that you display in your Question is not exactly the same as the one attached to the Question. For example, the displayed one defines a function TS1, but in the attached it is TS[1]. Please attach the most recent copy of the worksheet.