+++ /dev/null
-module Main
-where
-
-import FEM.R1
-import Linear.Matrix hiding ((*))
-import Naturals
-import Piecewise
-
-big_A' :: Double -> Double
-big_A' = const 1
-
-c' :: Double -> Double
-c' x = sin x
-
-f' :: Double -> Double
-f' x = x*(sin x)
-
-bdy' :: BoundaryConditions Double
-bdy' = Left (Dirichlet (0,1))
-
-
-pde' :: PDE Double
-pde' = PDE big_A' c' f' bdy'
-
-mesh' :: Col4 (Double,Double)
-mesh' = fromList [[(0.0,1/3)],
- [(1/3,2/3)],
- [(2/3,4/5)],
- [(4/5,1.0)]]
-
-pvec' :: Col4 Int
-pvec' = fromList [[2],[3],[2],[1]]
-
-params' :: Params N4 N3 N7 Double
-params' = Params mesh' pvec'
-
-
-soln :: Piecewise Double
-soln = solution pde' params'