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