--- >>> import Linear.Matrix ( Col4, Col7, frobenius_norm, fromList )
--- >>> import Naturals ( N3, N4, N7 )
---
--- >>> let big_A = const (1::Double)
--- >>> let c x = sin x
--- >>> let f x = x*(sin x)
--- >>> let bdy = Left (Dirichlet (0,1::Double))
--- >>> let pde = PDE big_A c f bdy
---
--- >>> let i1 = (0.0,1/3)
--- >>> let i2 = (1/3,2/3)
--- >>> let i3 = (2/3,4/5)
--- >>> let i4 = (4/5,1.0)
--- >>> let mesh = fromList [[i1], [i2], [i3], [i4]] :: Col4 (Double,Double)
--- >>> let pvec = fromList [[2],[3],[2],[1]] :: Col4 Int
--- >>> let params = Params mesh pvec :: Params N4 N3 N7 Double
+-- >>> import Linear.Matrix ( Col7, frobenius_norm )
+-- >>> import FEM.R1.Example1 ( pde', params' )