MATLAB:

function net = NNforward(net,input,targets)
  net.in = input;                
  bias = ones(size(input,1),1);
  net.hid = NNsigmoid([input bias]*net.wh); 
  net.out = NNsigmoid([net.hid bias]*net.wo); 
  net.error = targets - net.out;

function net = NNupdatewts(net)
  net.do = net.error .* net.out .* (1-net.out);
  net.dh = net.do * net.wo(1:net.nh,:)' ...
                          .* net.hid .* (1-net.hid);
  net.wo = net.wo + net.ro * [net.hid 1]' * net.do;    
  net.wh = net.wh + net.rh * [net.in 1]' * net.dh;

function y = NNsigmoid(x)
  y = 1 ./ (1 + exp(-x));