function L = Ucomplap (nodes,Nnodes,elements,Nelements,coeff) % L = Ucomplap (nodes,Nnode,elements,Nelements,coeff) % Computes the P1 finite element approximation of the % differential operator - d ( coeff d (.)\dx)\dx ## This file is part of ## ## SECS1D - A 1-D Drift--Diffusion Semiconductor Device Simulator ## ------------------------------------------------------------------- ## Copyright (C) 2004-2007 Carlo de Falco ## ## ## ## SECS1D is free software; you can redistribute it and/or modify ## it under the terms of the GNU General Public License as published by ## the Free Software Foundation; either version 2 of the License, or ## (at your option) any later version. ## ## SECS1D is distributed in the hope that it will be useful, ## but WITHOUT ANY WARRANTY; without even the implied warranty of ## MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the ## GNU General Public License for more details. ## ## You should have received a copy of the GNU General Public License ## along with SECS1D; If not, see <http://www.gnu.org/licenses/>. h = nodes(2:end)-nodes(1:end-1); d0 = ; d1 = [1000; -coeff./h]; dm1 = ; L = spdiags([dm1, d0, d1],-1:1,Nnodes,Nnodes); % Last Revision: % $Author: adb014 $ % $Date: 2008-02-04 16:26:27 +0100 (man, 04 feb 2008) $

Generated by Doxygen 1.6.0 Back to index