pikoe, Proton-Induced KnockOut reaction calculation for Exclusive processes, is a Fortran 90 program that calculates triple- and quadruple-differential cross sections, vector analyzing powers, and momentum distributions of reaction residues, for proton-induced nucleon knockout reactions in normal and inverse kinematics. The distorted-wave impulse approximation with the factorization approximation to the nucleon-nucleon transition matrix is adopted, and the distorted waves of the incoming proton and the outgoing two nucleons are calculated quantum mechanically. Kinematics of the reaction particles are treated in a relativistic manner, which gives the proper asymptotics of the three-body scattering wave in the plane-wave limit. Program summaryProgram Title: pikoeCPC Library link to program files:https://doi.org/10.17632/m594h58kck.1Licensing provisions: MITProgramming language: Fortran 90Supplementary material: See https://www.rcnp.osaka-u.ac.jp/~kazuyuki/pikoe/index.phpNature of problem: Calculation of triple- and quadruple-differential cross sections, vector analyzing powers, and momentum distributions of reaction residues, for nucleon knockout reactions with proton beams/targets in normal/inverse kinematics.Solution method: Numerov method for solving single-channel second-order differential equations; multi dimensional integration with the Gauss-Legendre quadrature rule.Additional comments including restrictions and unusual features: The global optical potential by Koning and Deraloche [1] for nucleons at energies from 1 keV to 200 MeV, for (near-)spherical nuclei with mass number between 24 and 209, is implemented in the program. Tables of nucleon-nucleon (NN) elastic differential cross sections and NN transition amplitudes calculated with the NN effective interaction by Franey and Love [2] are provided as external files for pikoe (included in the package). Users are encouraged to prepare tables for the optical potentials and the cross sections of the elementary process according to the reaction system of interest by themselves. With external files thus prepared, pikoe can be applied to various kinds of knockout processes such as α or deuteron knockout by proton, α-induced knockout reactions, etc. Calculation of spin observables, the vector analyzing power in fact, however, is restricted to nucleon-induced nucleon knockout processes with the coplanar kinematics.