Dear Tanya, If my understanding is correct your function returns more than one output argument. This is not allowed with the external_function syntax. The simplest approach is to call this function and update a1, a2, a3 and a4 in the steadystate file. I posted an example in this thread last year.
Best, Stéphane.