Найти производную скалярного поля u = 3x^2 + x^3y + 5x + ln(z + z^2) в точке M_0 (1;3;1) по направлению от точки M_0 к точке M_1 (-2;2;1).