A few issues I can see,
-
You are not actually using sigma in your function at all and is
rather using it as a global variable so changing what you pass to the
rc0bc0gc0 function won’t make any difference. -
You are using sigma (as well as other values) as scalar. As I’ve
already said, you should figure out if you want to keep it that way
(and apply rc0bc0gc0 on each of the scalar values individually) or if
you want to “vectorize” rc0bc0gc0.The two should be almost equivalent though one will give you an
array of 3-tuples (since that’s the return value of rc0bc0gc0) and the
other will give you three arrays.
After you fix the function parameter (i.e. 1), you can fix this by
adding.
to appropriate places (i.e. where you want to apply element-wise
operation on the array)