Mercurial > fem-fenics-eugenio
view example/Heat.ufl @ 63:3231f888f25d
New examples for the fem-fenics pkg
* Laplace: solve the laplace problem described on dolfin website
* Biharmonic: solve the biharmonic equation described on dolfin
* Heat: solve the heat evolutionary equation
author | gedeone-octave <marco.vassallo@outlook.com> |
---|---|
date | Fri, 26 Jul 2013 11:29:29 +0200 |
parents | |
children |
line wrap: on
line source
element = FiniteElement("Lagrange", triangle, 1) u = TrialFunction(element) v = TestFunction(element) u0 = Coefficient(element) f = Coefficient(element) g = Coefficient(element) k = Constant(triangle) eq = (1/k)*(u-u0)*v*dx + inner(grad(u), grad(v))*dx -f*v*dx - g*v*ds a = rhs (eq) L = lhs (eq)