Hello, I am new to Julia and trying to convert one of the tasks that I used to do in Excel into Julia.

I need to find a matrix X. Matrix X represents discrete Markov model transition probabilities with the final column as a default state.

I know the default probabilities so year 1 to 30 so I am trying to solve for the Matrix X.

I am using JuMP and defining my variable matrix as follows:

```
@variable(model, 0 <= x[1:8, 1:8] <= 1)
```

Now, let’s say I want to define the probability of being in a default state after 3 years. It would be given by x ^ 3. However, Julia / JuMP does not like this operation. and I get the following error:

Ultimately, I want to use this to calculate the NL objective to solve.

On a side note, x*x works (x^2 doesn’t)

If you are interested in the problem I am trying to solve - it’s here: