function network = apply_gradients(network, gradients, learning_rate)
    network.Learnables = network.Learnables - learning_rate * gradients;
end