Input: The single-valued neutrosophic decision matrix ${\left(R_{ij}^l\right)}_{n\times m}\left(l=1,2,\cdots ,z\right)$ , the alternatives is $S=\left\{S_1,S_2,\cdots ,S_n\right\}$ , and the distinguish coefficient is $\varsigma =0.5$ .

Output: The score function and the accuracy function values of each alternative.

function the score function and the accuracy function values of each alternative $\left({\left(R_{ij}^l\right)}_{n\times m},\varsigma ,\rho \right)$

for each $\text{DM}=1\to z$ do

function the general decision matrix ${\left({\delta }_{ij}\right)}_{n\times m}$

calculate the aggregated decision matrix ${\left({\delta }_{ij}\right)}_{n\times m}$ in Equation (11)

end for

return ${\left({\delta }_{ij}\right)}_{n\times m}$

end function

function the weight of criteria $w_j$

for each criterion $j=1\to m$ do

calculate the arithmetic mean value of the alternative under criteria ${\overline{\delta }}_i$

end for

calculate the ${\overline{\delta }}_i$ in Equation (12)

return ${\overline{\delta }}_i$

for each alternative under criterion $i,j=1\left(i=1,2,\cdots ,n;j=1,2,\cdots ,m\right)$ , determine the grey relation coefficient $\tau \left({\overline{\delta }}_i,{\delta }_{ij}\right)$ do

Calculate the $\tau \left({\overline{\delta }}_i,{\delta }_{ij}\right)$ in Equation (13), calculate the grey correlation degree $\pi \left({\overline{\delta }}_i,{\delta }_{ij}\right)$ in Equation (14), calculate the weight of criteria $w_j$ in Equation (15)

end for

return $w_j$

end function

for each alternative $i\leftarrow zeros\left(1,n\right)$ do

calculated the aggregated overall values ${\delta }_i$ of each alternative in Equation (16), calculated the score function and accuracy function values of each alternative $p\left({\delta }_i\right)$ and $q\left({\delta }_i\right)$ in Equation (17)

end for

return $p\left({\delta }_i\right)$ and $q\left({\delta }_i\right)$

end function