I’m trying to run a simulation. I will number sentences to make response easier. (1) Here is my first equation:

`𝑗[LBar_, y_, x1_, σ_, X2M_, w_, X1M_]:=(LBar(𝜎−1)⁄𝜎 + 𝑦(𝜎−1)⁄𝜎−x1(𝜎−1)⁄𝜎)𝜎⁄(𝜎−1)−(𝑦−X2M+𝑤(LBar+X1M)−𝑤x1) `

(2) I insert some parameter values and Reduce:

`Reduce[𝜕x1(𝑗[100,𝑦,x1,13⁄,42,.23,80])==0 && x1>0, x1, Reals] `

(3) This produces output:

`x1 == 119.575 `

(4) When I try to call the output value (119.575) I run into problems. (5) For instance

`2 % `

results in:

`2 (x1 == 119.575) `

(6) and

`j[100, 80, %[[2]], 1/3, 42, .23, 80] `

produces this output:

`-79.4 + 1/Sqrt[41/160000 - 1/239.149[[2]]^2] + 0.23 239.149[[2]] `