# data are from Bonner and Schwarz (2006)
list(X=structure(.Data = c(1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,0,1,1,1,1,1,1,1,0,1,0,0,0,1,1,1,1,1,1,1,0,1,1,1,1,1,1,1,1,1,1,0,0,1,1,1,1,1,0,0,0,1,1,0,0,1,0,0,0,1,1,1,0,1,0,0,0,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,0,0,0,1,1,1,1,1,1,1,1,1,1,1,0,1,1,1,0,0,1,1,0,1,0,0,0,1,0,0,0,1,0,0,0,1,0,0,0,1,1,1,1,1,0,0,0,1,1,1,1,0,0,1,1,1,1,1,1,1,1,0,0,1,0,0,0,1,0,0,0,1,0,1,0,1,0,0,0,1,1,1,1,1,1,1,0,1,0,1,0,1,1,1,1,1,1,1,1,1,0,0,0,1,1,1,1,1,0,0,0,1,0,0,0,1,1,1,1,1,1,1,1,1,0,0,0,1,1,0,0,1,1,1,0,1,1,1,0,1,1,1,1,1,0,0,0,1,1,1,1,0,1,1,1,1,1,1,0,1,1,1,1,1,0,0,0,0,1,1,1,0,1,1,1,1,0,0,0,1,0,0,0,1,0,0,0,1,0,0,0,0,1,1,1,1,0,0,0,0,0,1,1,0,1,1,0,0,0,1,1,1,1,1,1,1,0,0,0,1,0,0,0,0,1,1,1,1,1,1,1,1,0,0,0,0,1,1,0,1,0,0,0,0,1,1,0,0,1,1,1,1,0,0,0,1,0,0,0,1,0,0,0,1,1,0,0,0,1,1,1,0,1,0,0,0,1,1,0,1,1,0,0,0,1,1,1,0,1,1,1,1,0,0,0,1,0,0,0,0,1,1,1,1,0,0,0,1,0,0,0,1,1,0,1,1,0,0,1,0,0,1,1,0,0,1,0,0,0,1,0,0,1,1,1,0,1,1,1,0,1,1,1,0,1,0,0,0,0,1,1,0,0,1,1,0,1,1,1,0,1,1,1,0,1,0,0,0,1,1,1,0,0,1,1,0,1,1,1,0,0,1,1,0,0,1,1,0,0,1,0,0,0,1,1,0,1,1,1,0,0,1,1,0,0,1,1,0,0,1,1,0,0,1,0,0,0,1,1,0,1,0,0,0,1,1,1,0,1,1,1,0,0,1,1,0,1,1,1,0,0,1,1,0,0,1,1,0,1,0,1,0,0,1,1,0,0,1,1,0,0,1,1,0,0,1,1,0,0,1,1,0,1,0,0,0,0,1,1,0,1,0,0,0,0,1,1,0,1,1,1,0,0,1,1,0,1,1,0,0,0,1,1,0,0,1,1,0,0,1,1,0,0,1,1,0,0,1,1,0,0,1,0,0,0,1,1,0,0,1,0,0,0,1,1,0,0,1,1,0,0,1,1,0,0,1,1,1,1,1,1,0,0,1,0,0,0,1,1,0,0,1,0,0,0,1,0,0,0,1,1,0,0,1,0,0,0,1,1,0,0,1,1,0,0,1,0,0,0,1,1,0,0,1,1,0,0,1,0,0,0,1,0,0,0,1,1,0,0,1,1,0,0,1,0,1,0,1,0,0,0,1,1,0,0,1,0,0,0,1,1,0,0,1,0,1,1,1,1,0,0,1,1,1,1,1,1,1,1,1,0,1,0,0,0,1,1,1,1,1,1,1,1,1,0,0,0,1,1,1,1,1,1,1,1,1,1,0,0,0,1,1,0,1,0,1,1,0,1,0,1,0,1,1,1,0,1,1,1),.Dim=c(199,4)),
z=structure(.Data = c(2.219700,0.78011,0.57446,0.67729,2.2197,2.2197,2.2197,2.2197,2.2197,2.2197,2.2197,2.2197,2.2197,2.2197,2.2197,2.2197,0.88294,2.2197,2.2197,NA,0.57446,0.47163,0.26598,0.16315,0.26598,0.16315,-0.24816,NA,1.0886,NA,NA,NA,1.1914,0.98577,1.0886,0.78011,0.57446,0.57446,-0.042502,NA,0.16315,0.98577,-0.14533,0.36881,0.57446,0.16315,0.47163,0.16315,0.16315,0.060325,NA,NA,1.1914,0.88294,0.67729,0.78011,0.78011,NA,NA,NA,0.26598,0.26598,NA,NA,0.47163,NA,NA,NA,0.78011,1.0886,0.16315,NA,0.67729,NA,NA,NA,1.1914,1.0886,1.3971,0.57446,0.57446,0.060325,0.16315,0.060325,-0.042502,-0.45381,-0.14533,-0.76229,1.1914,0.78011,0.78011,0.57446,2.2197,1.3971,0.57446,0.67729,0.78011,NA,NA,NA,1.0886,0.78011,0.78011,0.78011,-0.14533,0.57446,-0.042502,-0.14533,0.47163,0.26598,0.16315,NA,1.6027,2.2197,0.88294,NA,NA,0.26598,0.36881,NA,1.2942,NA,NA,NA,0.67729,NA,NA,NA,0.47163,NA,NA,NA,0.98577,NA,NA,NA,1.0886,0.47163,0.78011,0.57446,-0.55664,NA,NA,NA,0.47163,0.67729,1.0886,0.98577,NA,NA,1.0886,2.2197,0.67729,0.98577,0.16315,0.36881,1.0886,1.2942,NA,NA,1.1914,NA,NA,NA,-0.042502,NA,NA,NA,0.88294,NA,0.36881,NA,2.2197,NA,NA,NA,-0.14533,-0.35098,-0.14533,0.16315,2.2197,2.2197,2.2197,NA,1.0886,NA,0.78011,NA,0.16315,-0.24816,-0.45381,-0.55664,0.67729,0.78011,0.67729,0.47163,2.2197,NA,NA,NA,1.1914,2.2197,2.2197,2.2197,-0.86512,NA,NA,NA,0.67729,NA,NA,NA,0.36881,0.26598,-0.24816,-0.24816,-0.042502,0.26598,-0.042502,-0.35098,-0.35098,NA,NA,NA,0.98577,0.57446,NA,NA,0.78011,0.57446,-0.24816,NA,0.57446,0.16315,-0.042502,NA,-0.86512,-0.35098,-0.35098,-0.24816,0.98577,NA,NA,NA,-1.0708,-0.86512,-0.96794,-1.0708,NA,-0.65946,-0.55664,-0.35098,-0.042502,0.16315,-0.042502,NA,-0.14533,-0.042502,-0.35098,-0.24816,-0.76229,NA,NA,NA,NA,-0.55664,-0.55664,-0.24816,NA,-0.96794,-0.96794,-0.55664,2.2197,NA,NA,NA,-0.86512,NA,NA,NA,-1.5849,NA,NA,NA,0.16315,NA,NA,NA,NA,-1.3793,-1.2764,-1.0708,-0.45381,NA,NA,NA,NA,NA,-1.1736,-1.3793,NA,-1.2764,-0.86512,NA,NA,NA,-0.45381,-0.65946,-0.76229,-0.14533,-0.35098,-0.14533,-0.14533,NA,NA,NA,-1.1736,NA,NA,NA,NA,-0.96794,-0.76229,-0.86512,1.1914,2.2197,0.67729,1.1914,0.57446,NA,NA,NA,NA,-0.45381,-0.35098,NA,1.3971,NA,NA,NA,NA,-1.2764,-1.1736,NA,NA,-1.2764,-1.3793,-0.96794,-1.5849,NA,NA,NA,-1.1736,NA,NA,NA,0.47163,NA,NA,NA,1.0886,1.2942,NA,NA,NA,-0.96794,-0.65946,-0.86512,NA,-1.4821,NA,NA,NA,-1.2764,-0.96794,NA,-1.2764,-0.35098,NA,NA,NA,-1.3793,-1.1736,-0.65946,NA,-1.1736,-0.76229,-0.65946,-0.14533,NA,NA,NA,1.0886,NA,NA,NA,NA,-0.86512,-0.65946,1.9112,-0.76229,NA,NA,NA,0.78011,NA,NA,NA,-0.65946,-0.35098,NA,-0.042502,0.16315,NA,NA,0.57446,NA,NA,-0.96794,0.16315,NA,NA,-0.86512,NA,NA,NA,-1.1736,NA,NA,-0.65946,-0.45381,0.26598,NA,-0.65946,0.16315,0.16315,NA,-0.65946,-0.45381,-0.042502,NA,-0.042502,NA,NA,NA,NA,-0.86512,-1.0708,NA,NA,-1.1736,0.060325,NA,-1.1736,-0.86512,-0.65946,NA,-1.2764,-1.0708,-1.0708,NA,-0.86512,NA,NA,NA,-1.5849,-0.96794,-0.86512,NA,NA,-0.96794,-0.96794,NA,-0.96794,-0.76229,-0.76229,NA,NA,-1.0708,-0.55664,NA,NA,-0.76229,-0.35098,NA,NA,-0.86512,NA,NA,NA,-1.2764,-1.0708,NA,-0.14533,0.26598,0.67729,NA,NA,-0.55664,-0.24816,NA,NA,-0.76229,-0.86512,NA,NA,-1.1736,-0.45381,NA,NA,-1.2764,NA,NA,NA,-1.1736,-0.76229,NA,-1.1736,NA,NA,NA,-0.86512,-0.14533,-0.35098,NA,-1.0708,-0.86512,-1.0708,NA,NA,-0.86512,-0.35098,NA,-1.4821,-0.55664,-0.86512,NA,NA,-0.55664,-0.45381,NA,NA,0.26598,-0.76229,NA,-0.65946,NA,1.8084,NA,NA,-0.96794,-0.96794,NA,NA,-0.96794,-0.96794,NA,NA,-1.0708,-0.96794,NA,NA,-1.0708,-0.86512,NA,NA,-0.76229,-0.65946,NA,0.47163,NA,NA,NA,NA,-1.3793,-1.2764,NA,-1.2764,NA,NA,NA,NA,-1.2764,-0.86512,NA,-1.5849,-0.76229,-0.96794,NA,NA,-1.0708,-0.86512,NA,-1.5849,-0.24816,NA,NA,NA,-1.1736,0.16315,NA,NA,-0.96794,-0.35098,NA,NA,-0.96794,-0.65946,NA,NA,-0.76229,-0.65946,NA,NA,-1.2764,-1.2764,NA,NA,-1.2764,NA,NA,NA,-1.4821,-1.4821,NA,NA,-0.45381,NA,NA,NA,-0.86512,0.67729,NA,NA,-1.3793,-0.45381,NA,NA,-0.65946,-0.45381,NA,NA,-1.1736,-0.35098,0.88294,0.78011,0.26598,0.98577,NA,NA,0.16315,NA,NA,NA,-0.96794,0.36881,NA,NA,-0.96794,NA,NA,NA,-0.96794,NA,NA,NA,-0.86512,-0.45381,NA,NA,-0.45381,NA,NA,NA,-1.1736,1.1914,NA,NA,0.67729,0.67729,NA,NA,-0.76229,NA,NA,NA,-0.76229,-0.14533,NA,NA,-1.2764,-1.2764,NA,NA,0.36881,NA,NA,NA,-0.35098,NA,NA,NA,-1.1736,-1.1736,NA,NA,-0.96794,-1.0708,NA,NA,-0.65946,NA,0.47163,NA,1.0886,NA,NA,NA,-0.55664,-0.35098,NA,NA,-0.55664,NA,NA,NA,-1.3793,-0.65946,NA,NA,-0.14533,NA,0.16315,0.47163,0.16315,-0.042502,NA,NA,-0.86512,-0.76229,1.1914,0.98577,1.2942,0.67729,2.2197,0.57446,0.67729,NA,0.88294,NA,NA,NA,2.2197,2.2197,2.2197,0.98577,2.2197,0.98577,0.78011,0.78011,0.78011,NA,NA,NA,1.2942,0.78011,0.98577,0.78011,0.26598,0.26598,-0.24816,-0.86512,0.78011,0.16315,NA,NA,NA,-1.3793,-0.76229,NA,-0.24816,NA,-0.24816,-0.76229,NA,-1.2764,NA,-1.2764,NA,-1.3793,-1.0708,-1.1736,NA,-0.86512,-0.65946,-0.35098),.Dim=c(199,4)),
a=structure(.Data = c(1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,NA,1,1,1,1,1,1,1,NA,1,NA,NA,NA,1,1,1,1,1,1,1,NA,1,1,1,1,1,1,1,1,1,1,NA,NA,1,1,1,1,1,NA,NA,NA,1,1,NA,NA,1,NA,NA,NA,1,1,1,NA,1,NA,NA,NA,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,NA,NA,NA,1,1,1,1,1,1,1,1,1,1,1,NA,1,1,1,NA,NA,1,1,NA,1,NA,NA,NA,1,NA,NA,NA,1,NA,NA,NA,1,NA,NA,NA,1,1,1,1,1,NA,NA,NA,1,1,1,1,NA,NA,1,1,1,1,1,1,1,1,NA,NA,1,NA,NA,NA,1,NA,NA,NA,1,1,1,NA,1,NA,NA,NA,1,1,1,1,1,1,1,NA,1,1,1,NA,1,1,1,1,1,1,1,1,1,NA,NA,NA,1,1,1,1,1,NA,NA,NA,1,NA,NA,NA,1,1,1,1,1,1,1,1,1,NA,NA,NA,1,1,NA,NA,1,1,1,NA,1,1,1,NA,1,1,1,1,1,NA,NA,NA,1,1,1,1,NA,1,1,1,1,1,1,NA,1,1,1,1,1,NA,NA,NA,NA,1,1,1,NA,1,1,1,1,NA,NA,NA,1,NA,NA,NA,1,NA,NA,NA,1,NA,NA,NA,NA,1,1,1,1,NA,NA,NA,NA,NA,1,1,NA,1,1,NA,NA,NA,1,1,1,1,1,1,1,NA,NA,NA,1,NA,NA,NA,NA,1,1,1,1,1,1,1,1,NA,NA,NA,NA,1,1,NA,1,NA,NA,NA,NA,1,1,NA,NA,1,1,1,1,NA,NA,NA,1,NA,NA,NA,1,NA,NA,NA,1,1,NA,NA,NA,1,1,1,NA,1,NA,NA,NA,1,1,NA,1,1,NA,NA,NA,1,1,1,NA,1,1,1,1,NA,NA,NA,1,NA,NA,NA,NA,1,1,1,1,NA,NA,NA,1,NA,NA,NA,1,1,1,1,1,1,1,1,NA,NA,1,1,NA,NA,1,NA,NA,NA,1,NA,NA,1,1,1,NA,1,1,1,NA,1,1,1,NA,1,NA,NA,NA,NA,1,1,NA,NA,1,1,NA,1,1,1,NA,1,1,1,NA,1,NA,NA,NA,1,1,1,NA,NA,1,1,NA,1,1,1,NA,NA,1,1,NA,NA,1,1,NA,NA,1,NA,NA,NA,1,1,NA,1,1,1,NA,NA,1,1,NA,NA,1,1,NA,NA,1,1,NA,NA,1,NA,NA,NA,1,1,NA,1,NA,NA,NA,1,1,1,NA,1,1,1,NA,NA,1,1,NA,1,1,1,NA,NA,1,1,NA,NA,1,1,NA,1,1,1,NA,NA,1,1,NA,NA,1,1,NA,NA,1,1,NA,NA,1,1,NA,NA,1,1,NA,1,NA,NA,NA,NA,1,1,NA,1,NA,NA,NA,NA,1,1,NA,1,1,1,NA,NA,1,1,NA,1,1,NA,NA,NA,1,1,NA,NA,1,1,NA,NA,1,1,NA,NA,1,1,NA,NA,1,1,NA,NA,1,NA,NA,NA,1,1,NA,NA,1,NA,NA,NA,1,1,NA,NA,1,1,NA,NA,1,1,NA,NA,1,1,1,1,1,1,NA,NA,1,NA,NA,NA,1,1,NA,NA,1,NA,NA,NA,1,NA,NA,NA,1,1,NA,NA,1,NA,NA,NA,1,1,NA,NA,1,1,NA,NA,1,NA,NA,NA,1,1,NA,NA,1,1,NA,NA,1,NA,NA,NA,1,NA,NA,NA,1,1,NA,NA,1,1,NA,NA,1,NA,1,1,1,NA,NA,NA,1,1,NA,NA,1,NA,NA,NA,1,1,NA,NA,1,NA,1,1,1,1,NA,NA,1,1,1,1,1,1,1,1,1,NA,1,NA,NA,NA,1,1,1,1,1,1,1,1,1,NA,NA,NA,1,1,1,1,1,1,1,1,1,1,NA,NA,NA,1,1,NA,1,1,1,1,NA,1,1,1,NA,1,1,1,NA,1,1,1),.Dim=c(199,4)),
first=c(1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,2,1,1,1,1,1,1,1,3,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,2,1,1,1,2,2,1,1,1,1,2,1,3,2,3,1,1,1,2,1,1,2,1,2,2,1,1,1,1,2,2,2,1,2,2,1,1,2,1,1,1,1,3,3,3,2,2,2,2,3,3,2,2,2,2,3,2,3,3,3,3,2,3,3,3,3,3,2,2,2,3,2,3,3,2,3,3,3,3,3,2,3,2,3,2,3,2,3,3,3,3,3,3,3,3,3,3,3,3,1,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,1,3,3,3,3,1,3,1,1,1,1,1,1,1,1,1,2,1,2,2,2),
udot=199,k=4)