Hi,

I understand that the INTERFACE fzero allows to automatically use Broydn when the input is not a scalar.

I’ve created a copy of the FOC function in the stochastic growth model returning a vector of residuals and receiving a vector of guesses

```
function foc2(x_both)
use toolbox
implicit none
real*8, intent(in), dimension(2) :: x_both
real*8 :: foc2(2)
real*8 :: kplus, cplus,vplus(NS),Ev,xin1,xin2
integer :: is_p
[....]
! get first order condition
foc2(1) = xin1 - (beta*foc2(1))**(-gamma)
foc2(2) = xin2 - (1-beta)*x_in**(1d0-1d0/gamma)/(1d0-1d0/gamma)+beta*Ev**(1/theta)
end function
```

I then call it as done for fzero in the original code

```
call fzero(x_both, foc2, check) !
```

Already at that stage I get the warning from vscode (then error in compilation with ifort) that

The characteristics of the associated actual function result differ from the characteristics of the dummy function result

Yet I don’t see any discrepancy in declarations (real*8, dimension(2) for both input and function return).

Any suggestion/example on how to deal with this problem?

Thanks

Gianni