A computer code based on Method of Characteristics (MOC) is developed to solve neutron transport equation for mainly assembly level lattice calculation with reflective and periodic boundary conditions and to some extent core level calculation with vacuum boundary condition. The code is able to simulate square, circular and hexagonal geometries and their combinations. Delaunay triangulation together with the Bower–Watson algorithm is used to divide the problem geometry into triangular meshes. Ray tracing technique is developed to draw characteristics lines along different directions over the geometry and the transport equation is solved over these lines to obtain neutron flux distribution and multiplication factor for the geometry. A number of benchmark problems available in literature are analyzed to demonstrate the capability and validity of the code.