function neuron1(W, b, x)
	relu = r -> max(r, 0) # for real values
	z = W * x # complex for MRI
	return relu(abs(z) + b) .* sign(z)
end