FaCE is a self contained program, with namelist input, that solves the three body Faddeev equations. It enables the inclusion of excitation of one of the three bodies, whilst the other two remain inert. It is particularly useful for obtaining the binding energies and bound state structure compositions of light exotic nuclei treated as three-body systems, given the three effective two body interactions. A large variety of forms for these interactions may be defined, and supersymmetric transformations of these potentials may be calculated whenever two body states need to be removed due to Pauli blocking. Program summary Title of program: FaCE (Faddeev with Core Excitation) Catalogue identifier: ADTW Program summary URL: http://cpc.cs.qub.ac.uk/summaries/ADTW Program obtainable from: CPC Program Library, Queen's University of Belfast, N. Ireland Computers: The code is designed to run on any Unix/Linux workstation or PC. Operating systems: Linux or UNIX Program language used: Fortran-90 Numerical libraries used: Source code for 6 routines from the NAG and BLAS libraries is included to enable independent compilation. Memory required to execute with typical data: 9 Mbytes of RAM memory and 12 MB of hard disk space. No. of bits in a word: 32 or 64 No. of bytes in distributed program, including test data, etc.: 116 514 No. of lines in distributed program, including test data, etc.: 15 574 Distribution format: tar gzip file Nature of physical problem: The program calculates eigenenergies and eigenstates for the three body problem by solving the Faddeev equations. Method of solution: Given the two body effective potentials it performs the supersymmetric transformation in case where there are forbidden states to be removed. The three body wavefunction is expanded in hyperspherical coordinates, the hyper-angular part is a series of Jacobi polynomials and the hyper-radial part is written in terms of a Laguerre basis. Within this basis the three body matrix elements are calculated and the full three body Hamiltonian matrix is completed. The diagonalization process is performed after various reductions (isospin, orthonormal and Feshbach) to determine the energies. Finally the three body wavefunction is reconstructed and other bound state observables are calculated. Typical running time: 6 s on a 1.7 GHz Intel P4-processor machine.