Since the photon has no mass, longitudinal electromagnetic waves cannot exist, and ordinary electromagnetic waves are always transversal waves. However, by speculating a bit we could introduce in the lagrangian density of the electromagnetic field a mass term for the photon, I understand that this is precisely what the Proca action does:
$$\mathcal{L}=-\frac{1}{2}(\partial_\mu B_\nu^*-\partial_\nu B_\mu^*)(\partial^\mu B^\nu-\partial^\nu B^\mu)+\frac{m^2 c^2}{\hbar^2}B_\nu^* B^\nu +\mathcal{L}_m, \quad \text{with }B^\mu = \left (\frac{\phi}{c}, \mathbf{A} \right)$$
being $\phi$ the electric potential and $\mathbf{A}$ the magnetic vector potential. My qüestion is: What would be the form of Maxwell's equations of this "electromagnetic" field formed by photons with mass (including sources)?
May guess is that the homogeneous equations remain the same:
$$\boldsymbol{\nabla}\times \boldsymbol{E} = -\frac{\partial \boldsymbol{B}}{\partial t}, \qquad \boldsymbol{\nabla}\cdot\boldsymbol{B} = 0$$
but the equations containing the sources are:
$$\boldsymbol{\nabla}\cdot \boldsymbol{E} + \frac{m^2c^2}{\hbar^2}\boldsymbol{A} = \frac{\rho}{\epsilon_0}, \qquad \boldsymbol{\nabla}\times\boldsymbol{B} + \frac{m^2c^2}{\hbar^2}\boldsymbol{A} = \mu_0\left( \boldsymbol{j} + \frac{\partial \boldsymbol{E}}{\partial t}\right)$$