En esta tesis se presenta un nuevo código de cálculo numérico diseñado para obtener la evolución de la función de distribución de iones y electrones en dispositivos toroidales de fusión por confinamiento magnético. El trabajo viene motivado por la necesidad de estudiar aspectos inabordables para los códigos tradicionales de transporte. Con la presente tesis, se ha pretendido establecer una base sólida desde la cual se pueda acceder a los estudios que la han motivado: la difícil física de los plasmas del stellarator TJ-II. Como es lógico, la puerta quedaría igualmente abierta para estudiar cualquier otra máquina de geometría toroidal. Así pues, como mínimo se ha pretendido (i) mostrar claramente el fundamento teórico del código de cálculo; (ii) escoger y probar las herramientas numéricas necesarias; (iii) codificar consecuentemente las ecuaciones y el software de apoyo para el an álisis de sus resultados (p. ej. la interfaz con el usuario a través de herramientas avanzadas de visualización tridimensional); y (iv) contrastar los resultados de todo el ensamblaje con la teoría básica.