A novel reliability calculation algorithm is developed and applied to a robust optimization of electromagnetic device with uncertain design variables taken into account. The proposed algorithm, to calculate the reliability of a given design, incorporates sensitivity analysis based on finite element method with Monte Carlo simulation (MCS) method so that it has higher numerical efficiency than conventional methods: MCS and first-order reliability methods. Through numerical examples, its applicability to a reliability-based robust optimal design of electromagnetic device is demonstrated.