Sigma is now workig (test.m compute the numerical approximation and everything is ok).
First really working version.