A computational program for three dimensional stress analysis with boundary element method was developed on the basic of solution of Kelvin problem, and its availability was verified with some examples.