A new variant of the pseudospectral method for an activity-dependent dendritic branch model is proposed. This algorithm incorporates the Neumann boundary conditions in a more efficient way than in the algorithms proposed before for similar problems. Numerical experiments indicate that the new algorithm is more efficient than the previous algorithms discussed in the literature on the subject. (c) 2007 Elsevier B.V. All rights reserved.