Численное решение УЧП

[Gortaur]

Добрый день, первый раз в этом разделе форума, очень надеюсь на Вашу помощь.
У меня есть УЧП на $u(t,x)$ где $x$ имеет размерность от трех и выше. УЧП следующее:

$$ u_t = \sum\limits_{i=1}^n a_i(x)\frac{\partial u}{\partial x_i} + \frac{1}{2}\sum\limits_{i,j=1}^n b_ij(x)\frac{\partial^2 u}{\partial x_i\partial x_j} $$

например

$$ u_t = \Delta x. $$

У меня также есть начальные условия $u(0,x) = h(x)$. Граничные если необходимы могу подобрать тоже.
Какой софт лучше использоваться для численного решения такого УЧП на $t\in[0,T]$? Хочу попробовать MATLAB, там вроде есть подходящий toolbox, но разобраться не смог, да и говорят, что он не подходит для решения УЧП с большой размерностью пространственной переменной.