function [efcns,wn]=dirichlet(L,n)
 
  $Id: dirichlet.m,v 1.14 2003-01-05 14:43:15-06 brinkman Exp $
  (c) 2002, Peter Brinkmann (brinkman@math.uiuc.edu)
 
  dirichlet: finds eigenfunctions and eigenvalues for
 	Dirichlet boundary conditions, for use with coeff.m
 
  Usage: [ef,wn]=dirichlet(pi,25)
 
  Parameters:
 	L: max value in interval [0,L]
 	n: number of eigenvalues to be computed
 
  Return values:
 	efcns: name of function that computes eigenfunctions
 	wn: column vector of square roots of eigenvalues