Solución de ecuaciones diferenciales elípticas en Fortran

Buen día. Necesito un programa desarrollado en el lenguaje de Fortran que utiliza ecuaciones diferenciales elípticas, ya que usa 2 matrices de tamaño dinámico (definido por el usuario) en las cuales las condiciones de frontera es la primera fila en cada matriz (estos valores también son definidos por el usuario), los demás valores de las matrices son calculados por iteraciones, las fórmulas para las iteraciones depende de la posición del elemento dentro de las matrices, hay 16 zonas entre las dos matrices. Hay una zona en cada matriz en la cual la fórmula involucra elementos de la otra matriz.