FIX: call to intbymat4 removed - no use other than computing i's for each variable.
Included support for continuous case and changed Pxy,Ps