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...