I will present a computational method for the analysis and design of a temperature sensor based on optical fiber. I will combine symbolic and numerical computations using the following software: Maple, for symbolic computation; Ansys: and Quick-Field for the numerical-graph computation of temperature profiles; Opticfiber for the numerical-graph computation of the electromagnetic modes in the fiber optics. The design strategy is to convert the patterns of temperature in changes of the refractive index of the fiber and the detection of changes in the electro-optical normal modes in the fiber. The proposed method has many advantages for the design of optical fiber sensors nowadays, for temperature measurements as well for other physical variables. We will use many special functions of Mathematical Physics such as the error function, Bessel functions, Kummer functions, Heun functiions, Whittaker functions and Laguerre functions. We will use Maple to make very complex computations with such functions.