LM multivariate test of heteroskedasticity
[lmf,pvalue]=hetero_sq0_multi(res,np)
res = results tlist from a first stage estimation
np = unused argument (but put here for compatibility with other testing functions)
lmf = value of the statistic
pvalue = its p-value