Function files:
[a, bd, gas, nopt] = function inputdata(args)
. . . code . . .
end
[outputs] = function algorithm(a, bd, gas, nopt)
. . . code . . .
end
then in your main script:
[a, bd, gas, nopt] = inputdata(args);
[outputs] = algorithm(a, bd, gas, nopt);
