MaplePrimes - Newest Questions

Local curve distance sought

Alfred_F — Mon, 10 Aug 2026 14:31:32 Z

For the purpose of practicing with Maple, I have worked on a definite integral in the attached file. I found it on the Internet. According to a side calculation, it can only be calculated using famous special functions or series expansions. That is why I am only interested in the numerical evaluation here. According to auxiliary calculations, it can only be computed using famous special functions or a series expansion. It is known that the series expansion of the antiderivative converges very rapidly. Therefore, for comparison purposes, I calculated the values for integration limits of 10 and 100 instead of infinity. In the first graph, the two result curves appear to coincide perfectly. Is there a way to zoom in locally - like using a magnifying glass - to make the difference in the ordinate values visible?

Download testint3.mw

CPU vs GPU matrix calculations

Gillee — Mon, 10 Aug 2026 13:26:27 Z

Hello,

I found this sample code that determine the time to do matrix multiplications using cpus and gpus. I was surprised that the speedup factor was less than one. I also did a similar speed test using python pytorch code on the same computer. The speedup factor was remarkable. Please find my mistakes, thanks.

MAPLE CODE:

kernelopts(version);
with(LinearAlgebra);
with(CUDA);
n := 4000;
M1 := RandomMatrix(n, n, datatype = float[8]);
M2 := RandomMatrix(n, n, datatype = float[8]);
printf("CUDA enabled? %a\n", IsEnabled());
tNoCUDA := time[real](M1 . M2);
printf("Time without CUDA: %.3f seconds\n", tNoCUDA);
prevSetting := Enable(true);
printf("Previous CUDA setting: %a\n", prevSetting);
printf("CUDA enabled now? %a\n", IsEnabled());
tCUDA := time[real](M1 . M2);
printf("Time with CUDA: %.3f seconds\n", tCUDA);
printf("Speedup factor: %.2f\n", evalf(tNoCUDA/tCUDA));
props := Properties();
printf("CUDA Device Properties:\n");
props;
Enable(false);
printf("CUDA enabled after disabling? %a\n", IsEnabled());

RESULTS:

Maple 2026.0, X86 64 WINDOWS, Mar 05 2026, Build ID 2001916

CUDA enabled? false
Time without CUDA: 0.288 seconds
Previous CUDA setting: false
CUDA enabled now? true
Time with CUDA: 0.439 seconds
Speedup factor: 0.66
CUDA Device Properties:

[TABLE(["Clock Rate" = 2692000, "Resisters Per Block" = 65536,

"ID" = 0, "Texture Alignment" = 512,

"Max Grid Size" = [2147483647, 65535, 65535],

"Memory Pitch" = 2147483647, "Major" = 12,

"MultiProcessor Count" = 36, "Shared Memory Per Block" = 49152,

"Minor" = 0, "Name" = "NVIDIA GeForce RTX 5060 Ti",

"Max Threads Dimensions" = [1024, 1024, 64], "Warp Size" = 32,

"Total Global Memory" = 4294967295,

"Max Threads Per Block" = 1024, "Device Overlap" = 1,

"Total Constant Memory" = 65536,

"Kernel Exec Timeout Enabled" = true])]

CUDA enabled after disabling? false

Ordering of symbols in textplot

Michael — Mon, 10 Aug 2026 04:46:46 Z

Hello:

For educational purposes, I want to define the zeta function on a plot using textplot, so I write

textplot([2,3,Zeta(sigma+i*t)=Sum(1/n^(sigma+i*t),n=1..infinity)]...

But the characters "sigma+it" show up as "it+sigma" rather than "sigma+it" which is how most people would naturally read it.

Is there any way that I can convince textplot to order the output the way I set it?

I've tried single quotes, i.e. 'sigma+it' but that doesn't work.

Thank you.

Dual numbers and automatic differentiation

one man — Sat, 08 Aug 2026 14:01:20 Z

I learned about automatic differentiation from an AI. Could you please tell me, what does Maple have for working with dual numbers in this regard?
(I am experiencing some minor internet issues, so my response to your reply may be significantly delayed.)

Type checking for a kind of set

Ronan — Fri, 07 Aug 2026 00:48:53 Z

I am trying to add this to my type list in a package, but cannot get it to work.

TypeTools:-AddType(_L3DP, set(satisfies(s -> type(s, [algebraic $ 3])),'Vector[column](3, algebraic)'));

when I test with this I get an error.

 type({[7,8,9],<1,2,3>},:-_L3DP);
Error, testing against an invalid type

or this

type({[7,8,9],<1,2,3>},set(satisfies(s -> type(s, [algebraic $ 3])),'Vector[column](3, algebraic)'));
Error, testing against an invalid type

How make beautiful figures?

Gabriel Barcellos — Wed, 05 Aug 2026 18:33:15 Z

How do I vectorize a Maple graph—that is, export it in PDF quality—and add X and Y axis labels along with annotations?

By "annotations," I mean something like "green curve = 39T" written directly on the figure.

I need a more complete comand. After the comand, eu can ajust how I liked.