Maxima使うといいよ

taylor((y^2)*(a^x)+x*(a^(y^2)), [x,y], 0, 3);