A package for solving coupled systems of time varying partial differential equations in two-dimensional

We present a package for solving time varying coupled systems of partial differential equations on a rectangle. The numerical solution technique employs Galerkin's method in space, using B-splines and a variable order, variable time-step extrapolated backward difference procedure in time.