Continuous dependence of solutions to a fictitious domain formulation with respect to boundary variations is established. This property is used in shape optimization.