Este documento presenta una novedosa técnica numérica para resolver ecuaciones diferenciales estocásticas fraccionarias con retardos neutros. El método emplea un esquema de colocación por pasos con poli-fractonómicos de Jacobi para considerar procesos estocásticos desconocidos. Con este fin, las ecuaciones diferenciales de retraso se transforman en ecuaciones aumentadas sin retardos. Esta transformación permite utilizar un esquema de colocación mejorado con poli-fractonómicos de Jacobi para resolver las ecuaciones modificadas repetidamente. En cada iteración, se genera un sistema de ecuaciones no lineales. A continuación, se analizan rigurosamente las propiedades de convergencia del método propuesto. Posteriormente, la utilidad práctica de la técnica numérica propuesta se valida a través de una serie de ejemplos de prueba. Estos ejemplos ilustran la capacidad del método para producir soluciones precisas y eficientes.