Розглянуто метод радіальних функцій для наближеного розв'язання інтегро-диференціального рівняння з частинними похідними. Для апроксимації точного розв’язку використано лінійну комбінацію радіальних базисних функцій у певних центральних точках та лінійну комбінацію поліноміальних базисних функцій. Запропоновано розташування центральних точок як для дво-, так і для тривимірних областей. Колокація в центральних вузлах дає напівдискретну схему з інтегральними коефіцієнтами, для обчислення яких використано квадратурні формули Гаусса–Лежандра та трапецій. Для визначення параметра форми застосовано генетичний алгоритм. Наведено чисельні розрахунки для дво- та тривимірних областей, що підтверджують ефективність запропонованого алгоритму.