diff --git a/README.md b/README.md
index 5a45742..b438cba 100644
--- a/README.md
+++ b/README.md
@@ -131,10 +131,28 @@ double const *genann_run(genann const *ann, double const *inputs);
 Call `genann_run()` on a trained ANN to run a feed-forward pass on a given set of inputs. `genann_run()`
 will provide a pointer to the array of predicted outputs (of `ann->outputs` length).
 
+### Activation Functions
+
+Genann uses a sigmoid activation by default. Each network has
+`activation_hidden` and `activation_output` members which can be set to
+`genann_act_sigmoid_cached`, `genann_act_tanh`, `genann_act_relu`,
+`genann_act_linear`, or `genann_act_threshold`:
+
+```C
+genann *ann = genann_init(2, 1, 3, 2);
+ann->activation_hidden = genann_act_relu;
+```
+
+Backpropagation training knows the derivatives of the built-in activation
+functions only. If you substitute your own function, `genann_train()` will
+assume the sigmoid derivative; other training methods (see above) work
+with any activation.
 
 ## Hints
 
 - All functions start with `genann_`.
+- The default sigmoid activation expects outputs between 0 and 1. Scale
+  inputs to roughly ±1 for best results.
 - The code is simple. Dig in and change things.
 
 ## Extra Resources
diff --git a/genann.c b/genann.c
index b05fa4f..e55ef29 100644
--- a/genann.c
+++ b/genann.c
@@ -105,6 +105,14 @@ double genann_act_threshold(const struct genann *ann unused, double a) {
     return a > 0;
 }
 
+double genann_act_tanh(const struct genann *ann unused, double a) {
+    return tanh(a);
+}
+
+double genann_act_relu(const struct genann *ann unused, double a) {
+    return a > 0 ? a : 0;
+}
+
 genann *genann_init(int inputs, int hidden_layers, int hidden, int outputs) {
     if (hidden_layers < 0) return 0;
     if (inputs < 1) return 0;
@@ -275,6 +283,17 @@ double const *genann_run(genann const *ann, double const *inputs) {
 }
 
 
+/* Derivative of an activation function, in terms of its output value.
+ * Recognizes the built-in activations; any other function is assumed to
+ * have the sigmoid's derivative. */
+static double genann_act_derivative(genann_actfun act, double y) {
+    if (act == genann_act_tanh) return 1.0 - y * y;
+    if (act == genann_act_relu) return y > 0 ? 1.0 : 0.0;
+    if (act == genann_act_linear) return 1.0;
+    return y * (1.0 - y);
+}
+
+
 void genann_train(genann const *ann, double const *inputs, double const *desired_outputs, double learning_rate) {
     /* To begin with, we must run the network forward. */
     genann_run(ann, inputs);
@@ -296,7 +315,7 @@ void genann_train(genann const *ann, double const *inputs, double const *desired
             }
         } else {
             for (j = 0; j < ann->outputs; ++j) {
-                *d++ = (*t - *o) * *o * (1.0 - *o);
+                *d++ = (*t - *o) * genann_act_derivative(ann->activation_output, *o);
                 ++o; ++t;
             }
         }
@@ -328,7 +347,7 @@ void genann_train(genann const *ann, double const *inputs, double const *desired
                 delta += forward_delta * forward_weight;
             }
 
-            *d = *o * (1.0-*o) * delta;
+            *d = genann_act_derivative(ann->activation_hidden, *o) * delta;
             ++d; ++o;
         }
     }
diff --git a/genann.h b/genann.h
index e4b7383..bd7ed63 100644
--- a/genann.h
+++ b/genann.h
@@ -99,6 +99,8 @@ double genann_act_sigmoid(const genann *ann, double a);
 double genann_act_sigmoid_cached(const genann *ann, double a);
 double genann_act_threshold(const genann *ann, double a);
 double genann_act_linear(const genann *ann, double a);
+double genann_act_tanh(const genann *ann, double a);
+double genann_act_relu(const genann *ann, double a);
 
 
 #ifdef __cplusplus