Расчёты

Берём первую сф. умножаем на вторую, пропускаем через итерации для стационарной задачи, ставим в качестве начальной функции для Нестационарной в декартовых координатах.

Делаем расчёт с растяжением (hx *= 2) и сжатием (hx *= 0.5).

В первом случае в центре образуется область с малой температурой - значение в центре растёт ОЧЧЕНЬ медленно по сравнению с максимумом решения:

time = 0 center = 1.225234 sol_max = 1.859148

time = 0.5417 center = 2.689374 sol_max = 14.994066

time = 0.5810 center = 3.058009 sol_max = 120.199303

time = 0.5836 center = 3.088310 sol_max = 485.715073

При сжатии функция очень быстро вырождается в первую радиально-симметричную сф.

Неискажённая функция долго держит свою форму (досчитал до sol_max = 6000), хотя значение в центре всё равно сильно отстаёт.