for (out = 0; out < network.size.output; out++) {
        network.error.output[out] 
            = (network.neuron.target[out] - network.neuron.output[out]) 
              * sigmoid_derivative(network.neuron.output[out]);
    }