I have a huge set of coupled nonlinear integro-partial differential equations. After a long while trying to simplify the equations and solve them at least semi-analytically I have come to conclude there has been left no way for me but an efficient numerical method. Finite element seems most amenable as is based on Galerkin method which gives a weak form solution, so a great hope that it might finally solve the equations. But at the same time I am so new to this field to write the codes all from the scratch. Then I found FEniCS but everywhere I only read about FEniCS solving PDEs. Now I am interested if FEniCS/Dolphin can also solve integro-differential equations?
Here are the equations in my simplest case, just as example:
wherein, $i=1,2,3$ and $j_1,j_2=1,2,3,4$, thus, the number of unknowns and equations is: $$\overbrace{4}^{\text{components of }U^0}+\overbrace{4^2}^{\text{components of }U^1}+\overbrace{4^3}^{\text{components of }U^2}=84$$