data = (
{
  class = "mpdf_harness";
  mpdf = {
    class = "mprod";
    mpdfs = ( 
    {
      class = "mgamma";
      beta = 2.0;
      k = 1.0;
      rv : 
      {
	class = "RV";
	names = ( "b" );
      };
      rvc : 
      {
	class = "RV";
	names = ( "a" );
      };
    }, 
    {
	class = "egamma";
	alpha = 2.0;
	beta = 2.0;
	rv : 
	{
	  class = "RV";
	  names = ( "a" );
	};
    });
  };
  mean = [ 1.0, 1.0 ];
  R = ( "matrix", 2, 2, [ 2.0, 0.5, 0.5, 0.5 ] );
  variance = [ 1.2, 5.0 ];
  tolerance = 0.3;
},
{
  class = "mpdf_harness";
  mpdf = {
    class = "mprod";
    mpdfs = ( 
    {
	class = "enorm<ldmat>";
	mu = [ 1.5, 1.7 ];
	R = ( "matrix", 2, 2, [ 1.2, 0.3, 0.3, 5.0 ] );
	rv : 
	{
	  class = "RV";
	  names = ( "mprodx" );
	  sizes = [ 2 ];
	};
    }, 
    {
	class = "egamma";
	alpha = [ 100000, 10000 ];
	beta = [ 10000, 1000 ];
	rv = {
	  class = "RV";
	  names = ( "mprody" );
	  sizes = [ 2 ];
	};
    });
  };
  mean = [ 1.5, 1.7, 10.0, 10.0 ];
  variance = [ 5.0, 5.0, 5.0, 5.0 ];
});