function f = addNode(net)

ro = length(net.left(1).w);
co =  length(net.right(1).w);
n = net.regressors + 1;

left(n+1).w = (2*rand(1, ro) - 1);
right(n+1).w = (2*rand(co, 1) - 1); 
d_left(n+1,:) = zeros(1,ro);  %derivatives
d_right(n+1,:) = zeros(co,1);
dleft(n+1,:) = zeros(1,ro);   %deltas
dright(n+1,:) = zeros(co,1);
   

f = net;
f.regressors = n;

f.left = left;
f.right = right;
f.d_left = d_left;
f.d_right = d_right;
f.dleft = dleft;
f.dright = dright;

f.d_b = zeros(1,n);
f.weights = 2*rand(1,n+1) - 1;
f.dweights = zeros(1,n+1);
f.d = zeros(1,n+1);
f.bias = (2*rand(1,n) - 1);
f.dbias = zeros(1,n);
f.numparams = n+1 + n*(ro + co + 1);