Ниже привожу код в Maple 2022, а далее описываю проблему.
restart;
with(plots);
a := 0.5;
j := BesselJ(0, sqrt(2*I)*sqrt(x^2 + y^2)*a)/BesselJ(0, sqrt(2*I)*a);
A := sqrt(Re(j)^2 + Im(j)^2);
densityplot(A, x = -1 .. 1, y = -sqrt(-x^2 + 1) .. sqrt(-x^2 + 1), grid = [68, 68], style = patchnogrid, scaletorange = 0 .. 1);
Таким образом, имею вещественную функцию двух переменных:

. В нуле она минимальна-(значение 0.996). На границе единичного круга максимальна со значением равным 1. По идее, scaletorange = 0 .. 1 должна дать картнку со всюду в круге практически идеально белым так как значения функции очень близки к 1. Однако что со scaletorange, что без него команда выдает картинку с существенно темным кружком в центре...