This is not yet implemented. It is cited from [diffmat.m]() and should look like [poisson2D.m]() and [poisson3D.m]().