Добрый вечер! Передо мной стоит задача минимизировать целевую функцию двумерного ур-я теплопроводности, чтобы к заданному моменту

распределение температуры было как можно ближе к распределению

.

Что я сделал:
1) Дискритизировал данный интеграл методом трапеций
2) Вычислил значение целевой функции в конечный момент

Пока не смог разобраться как в данном случае применить алгоритм глобальной оптимизации с помощью градиентного спуска. Например, как вычислить градиент целевой функции и так далее. Надеюсь, вы сможете меня направить в правильном направлении, спасибо!