Я так понимаю, функция
у Вас — это гладкая "шапочка", которая "живёт" в районе нуля (в интервале
, вне этого интервала нуль). Если
возрастает в окрестности нуля от 0 до 1 (или от -1 до 1), то из неё можно сделать эту "шапочку", продифференцировав. Затем её надо нормировать (домножить на константу) так, чтобы интеграл от неё по вышеуказанному отрезку был равен 1.
Стандартная процедура для того, чтобы сделать из разрывной функции гладкую — сглаживание с помощью свёртки с такой вот гладкой "шапочкой". Пусть
— Ваша кусочно-постоянная функция концентраций. Тогда сглаженная
будет считаться так:
Дальше надо подставлять, интегрировать, упрощать. Может быть, проще для понимания окажется один раз сделать интегрирование для функции, равной единице в одном квадратике и нулю в остальных, получить формулу, а затем просуммировать такие формулы по всем квадратикам с коэффициентами-концентрациями в этих квадратиках. Контрольное отношение для понимания: в точках, отстоящих от границ более чем на
,
должна совпадать с
.