Решить численно двумерную краевую задач

Есть такая задача, двумерная, нужно решить ее в какой-нибудь программе (mathlab,wolfram mathematica), численно, и построить график (можно двумерный с градиентами либо трехмерный).
функции А и В можно взять любые, краевые условия можно взять любые числа (функции), с начальным условием все несколько сложнее, в подобных одномерных задачах обычно берется u_init =1/2 * (u_1(y,0)-u_0(y,0)) tanh((y-0.25)/epsilon) +1/2 * (u_1(y,0)+u_0(y,0)) , в двумерном что-то подобное, только с иксом еще где-нибудь, видимо. но на этом можно сильно не останавливаться, мне важен рабочий код, как решать задачи такого вида. чтобы я мог подставлять свои функции и смотреть результат.