This study is to develop a program, with which the motion of the system can be simulated in view of kinematics and dynamics, and then evaluate functional capabilities of many robots in order to select the most suitable one for any given working conditions. According to the theory of Denavit & Hartenberg, rectangular coordinate is set at each joint in the robot of six joints, as this one is the major concern, and transformation matrices are derived between any two coordinates. The necessary displacement of each joint for executing a given a work is obtained from the position and posture transformation and the in-stantaneous velocity of the joint is got from given velocity distribution of the path. The necessary torque of each joint is finally calculated by Newton-Euler Method. This simulation program is applied to a robot with six rotational joints, Rotational angle and torque to time which has a given motion trajectory are figured and also torque to r.p.m.