data = (
{
  class = "epdf_harness";
  epdf = {
    class = "egiw";
    fV = ( "matrix", 2, 2, [ 13.1, 11.0, 11.0, 10.0 ] );
// m=1.1; var_e = 0.1;
    nu = 10.0;
    dimx = 1;
    rv = {
      class = "RV";
      names = ( "a", "b" );
    };
  };
  mean = [ 1.1, 0.2 ];
  variance = [ 0.02, 0.01333 ];
  R = ( "matrix", 2, 2, [ 0.02, 0.0, 0.0, 0.022 ] );
  support = {
  	class = "rectangular_support";
  	ranges = ( [ -2.0, 4.0 ], [0.01, 2.0 ] );
  	gridsizes = [ 100, 100 ];
	};
  nbins = [ 100, 200 ];
  tolerance = 0.01;
},
{
  class = "epdf_harness";
  epdf = {
    class = "egiw";
    fV = ( "matrix", 2, 2, [ 200.0, 80.0, 80.0, 40.0 ] );
	// m=2; var_e = 1;
    nu = 40.0;
    dimx = 1;
    rv : 
    {
      class = "RV";
      names = ( "a", "b" );
    };
  };
  mean = [2.0, 1.14286];
  variance = [0.0285714, 0.0395795];
  R = ( "matrix", 2, 2, [ 0.028, 0.0, 0.0, 0.078 ] );
  support = {
  	class = "rectangular_support";
  	ranges = ( [ 0.0, 4.0 ], [0.01, 3.0 ] );
  	gridsizes = [ 100, 100 ];
	};
  tolerance = 0.01;
},
{
  class = "epdf_harness";
  epdf = {
    class = "egiw";
    fV = ( "matrix", 2, 2, [ 10.0, 5.0, 5.0, 15.0 ] );
    nu = 16.0;
    dimx = 2;
    rv : 
    {
      class = "RV";
      names = ( "a", "b" );
    };
  };
  mean = [1.0, 0.5, 0.5, 1.5];// V/(delta-2-1)=V/(nu-2-1-2-1)
  variance = [0.2500, 0.2045, 0.2045, 0.5625]; // m=13; var(i,j)=(12Vij^2+10Vii*Vjj)/(11*100*8)
  R = ("matrix",4,4,[0.2500, 0.0, 0.0, 0.0, 0.0, 0.2045, 0.0, 0.0, 0.0, 0.0, 0.2045, 0.0, 0.0, 0.0, 0.0, 0.5625]); 
  support = {
  	class = "discrete_support";
  	points = "s@egiw_support.cfg";
	volumes = "v@egiw_support.cfg";
	};
  tolerance = 0.03;
}
);