pyNeutronics
Solution of neutron diffusion problem for 1 energy group in a very simple 1-dimensional slab geometry with vacuum boundary conditions.
Finite difference box integration method is used for spatial discretization.
It doesn't get any simpler than this...