Consider the $1$-parameter family of functional equations on the critical strip $$\zeta(1-s)\Gamma[\hat\Omega_r,1-s ]=\zeta(s)\Gamma[\Omega_r,s]$$
where the Gamma factor $$ \Gamma[\Omega_r,s]= 2\int_0^1e^{\frac{r}{\log(|x|)}}x^{s-1}dx=4\sqrt{\frac{r}{s}}\, K_1\left(2\, \sqrt{r s}\right)$$
satisfies the degenerate linear PDE
$$ r^2\frac{\partial^3}{\partial r^3}\Gamma[\Omega_r,s]=s^2\frac{\partial}{\partial s}\Gamma[\Omega_r,s]. $$
We can reduce this to second order through the Fourier transform