(* ::Package:: *)

(**********************************************************************)
(* data file lp_Q1-with-error-3.m *)
(*   *)
(* Last Updated: on June 26, 2017 by Hereman *)
(* Introduced an error for testing purposes: qqq in s factor should have been q, also KKK in laxPairMatrixL should have been k *)
(* or worse, later replace a x[0,0,0] by x[1,0,0] in s *)
(***  Q1 equation (ABS classification 2003) - with errors in both s and L  ***)

ddEQ = {p*(x[0,0,0] - x[0,1,0])*(x[1,0,0] - x[1,1,0])
        - q*(x[0,0,0] - x[1,0,0])*(x[0,1,0] - x[1,1,0])
        + d^2*p*q*(p-q)};
 
nameINPUT = "Q1 equation (ABS classification 2003) - with errors in given s and L";
ddeEquationListINPUT = ddEQ;

(* *)
laxPairMatrixL = {{KKK*(x[0,0,0]-x[1,0,0]) + p*x[1,0,0], -p*(d^2*k*(p-k) + x[0,0,0]*x[1,0,0])},
                  {p, k*(x[0,0,0]-x[1,0,0])-p*x[0,0,0]}};

laxPairMatrixM = {{k*(x[0,0,0]-x[0,1,0]) + q*x[0,1,0], -q*(d^2*k*(q-k) + x[0,0,0]*x[0,1,0])},
                  {q, k*(x[0,0,0]-x[0,1,0])-q*x[0,0,0]}};

explicitScalars = True;

tFunc = 1/(d*p + (x[0,0,0] - x[1,0,0]));
sFunc = 1/(d*q + (x[0,0,0] - x[0,1,0]));
(* correct sFunc = 1/(d*q + (x[0,0,0] - x[0,1,0])); *)
(* sFunc = 1/(d*qqq + (x[0,0,0] - x[0,1,0])); *)
(* x[1,0,0] in sFunc should have been x[0,0,0] *)
(* very bad sFunc = 1/(d*q + (x[0,0,0] - x[0,1,0])); *)

(*  Also work: 
tFunc = 1/(d*p - (x[0,0,0] - x[1,0,0]));
sFunc = 1/(d*q - (x[0,0,0] - x[0,1,0])); 
*)
(* *)

(* weightedParametersINPUT = {}; *)

