// Reference simulation with “simul_backward” and default “solve_algo” value

@#include "simul_backward_common.inc"

s = simul_backward_model(initialconditions, @{simlength}, innovations);

save simul_backward_ref.mat s