An efficient way to calculate 3-D stress distributions in silicon around embedded structures caused by different thermal expansion coefficients between silicon and inclusion is presented. The computational method is based on the solution for a parallelepiped, which is extended to arbitrarily shaped embedded structures.