Dear UQ community,

I am trying to solve a Bayesian inversion problem in which a transfer function is treated as the model. The transfer function is given by \frac{8}{s^3+2s^2+10s+8}. .

Following this, the model coefficients for the function can be represented by `X(:,1) = 1, X(:,2) = 2, X(:,3) = 10 and X(:,4) = 8`

. Using this, the .mFile is written as:

```
function V = uq_dynamic(X)
s = tf('s');
sys_closed = (X(:,4))./((X(:,1).*(s^3))+(X(:,2).*(s^2))+(X(:,3).*(s))+(X(:,4)));
t = linspace(0,10,400);
y = step(sys_closed,t);
V = y;
```

In doing so, I get an error in line 3: `Undefined function 'rdivide' for input arguments of type 'tf'`

. Is there any mistake in my vector representation of the model?

Best,

Vinit