A Radiation Therapy Design Problem consists in minimize the total radiation dosage at the patient. In this work, the dosage values are represented by fuzzy numbers. The transition from healthy tissue to tumour cells is continuous and the surprise function is adopted to model it. We propose to solve this problem by a specially tailored Primal-Dual Interior Point Method and present numerical experiments with real world large-scale problems.