z偏x=-sin3xy*3y+1/(x+y+1)
z偏y=-sin3xy*3x+1/(x+y+1)
dz=[-sin3xy*3y+1/(x+y+1)]dx+[sin3xy*3x+1/(x+y+1)]dy