Neural networks are commonly used for classification tasks, in fact from this post it seems like that's where they shine brightest.
However, when we want to classify using neural networks, we often have the lastoutput layer to take values in $[0,1]$; typically, by taking the last layer to be the sigmoid function $x \mapsto \frac{e^x}{e^x +1}$. Is this theoretically justified
Can neural networks with a sigmoid as the activation function of the output layer approximate continuous functions? (i.e., isIs there an analogue to the universal approximation theorem for this case)? Why can the output of the neural network be in the range $[0,1]$ if it performs classification?