Two phase problems for linear elliptic operators with variable coefficients: Lipschitz free boundaries are C^{1,gamma}