test/Tools/isac/Knowledge/eqsystem.sml
branchisac-update-Isa09-2
changeset 37991 028442673981
parent 37984 972a73d7c50b
child 37995 fac82f29f143
equal deleted inserted replaced
37990:24609758d219 37991:028442673981
   863 "----------- me [EqSystem,normalize,2x2] -------------------------";
   863 "----------- me [EqSystem,normalize,2x2] -------------------------";
   864 val fmz = ["equalities [0 = -1 * q_0 * 0 ^^^ 2 / 2 + 0 * c + c_2,\
   864 val fmz = ["equalities [0 = -1 * q_0 * 0 ^^^ 2 / 2 + 0 * c + c_2,\
   865 	               \0 = -1 * q_0 * L ^^^ 2 / 2 + L * c + c_2]", 
   865 	               \0 = -1 * q_0 * L ^^^ 2 / 2 + L * c + c_2]", 
   866 	   "solveForVars [c, c_2]", "solution L"];
   866 	   "solveForVars [c, c_2]", "solution L"];
   867 val (dI',pI',mI') =
   867 val (dI',pI',mI') =
   868   ("Biegelinie.thy",["normalize", "2x2", "linear", "system"],
   868   ("Biegelinie",["normalize", "2x2", "linear", "system"],
   869    ["EqSystem","normalize","2x2"]);
   869    ["EqSystem","normalize","2x2"]);
   870 val p = e_pos'; val c = []; 
   870 val p = e_pos'; val c = []; 
   871 val (p,_,f,nxt,_,pt) = CalcTreeTEST [(fmz, (dI',pI',mI'))];
   871 val (p,_,f,nxt,_,pt) = CalcTreeTEST [(fmz, (dI',pI',mI'))];
   872 val (p,_,f,nxt,_,pt) = me nxt p c pt;val (p,_,f,nxt,_,pt) = me nxt p c pt;
   872 val (p,_,f,nxt,_,pt) = me nxt p c pt;val (p,_,f,nxt,_,pt) = me nxt p c pt;
   873 val (p,_,f,nxt,_,pt) = me nxt p c pt;val (p,_,f,nxt,_,pt) = me nxt p c pt;
   873 val (p,_,f,nxt,_,pt) = me nxt p c pt;val (p,_,f,nxt,_,pt) = me nxt p c pt;
   881 val (p,_,f,nxt,_,pt) = me nxt p c pt;f2str f(*["(bdv_1, c)", "(bdv_2, hd (tl [c, c_2] ... corrected srls; ran only AFTER use"RCODE-root.sml", store_met was NOT SUFFICIENT*);
   881 val (p,_,f,nxt,_,pt) = me nxt p c pt;f2str f(*["(bdv_1, c)", "(bdv_2, hd (tl [c, c_2] ... corrected srls; ran only AFTER use"RCODE-root.sml", store_met was NOT SUFFICIENT*);
   882 val (p,_,f,nxt,_,pt) = me nxt p c pt;f2str f;
   882 val (p,_,f,nxt,_,pt) = me nxt p c pt;f2str f;
   883 val (p,_,f,nxt,_,pt) = me nxt p c pt;f2str f;
   883 val (p,_,f,nxt,_,pt) = me nxt p c pt;f2str f;
   884 val (p,_,f,nxt,_,pt) = me nxt p c pt;f2str f;
   884 val (p,_,f,nxt,_,pt) = me nxt p c pt;f2str f;
   885 case nxt of
   885 case nxt of
   886     (_, Subproblem ("Biegelinie.thy", ["triangular", "2x2", "linear",_])) => ()
   886     (_, Subproblem ("Biegelinie", ["triangular", "2x2", "linear",_])) => ()
   887   | _ => raise error "eqsystem.sml me [EqSystem,normalize,2x2] SubProblem";
   887   | _ => raise error "eqsystem.sml me [EqSystem,normalize,2x2] SubProblem";
   888 
   888 
   889 val (p,_,f,nxt,_,pt) = me nxt p c pt;val (p,_,f,nxt,_,pt) = me nxt p c pt;
   889 val (p,_,f,nxt,_,pt) = me nxt p c pt;val (p,_,f,nxt,_,pt) = me nxt p c pt;
   890 val (p,_,f,nxt,_,pt) = me nxt p c pt;val (p,_,f,nxt,_,pt) = me nxt p c pt;
   890 val (p,_,f,nxt,_,pt) = me nxt p c pt;val (p,_,f,nxt,_,pt) = me nxt p c pt;
   891 val (p,_,f,nxt,_,pt) = me nxt p c pt;val (p,_,f,nxt,_,pt) = me nxt p c pt;
   891 val (p,_,f,nxt,_,pt) = me nxt p c pt;val (p,_,f,nxt,_,pt) = me nxt p c pt;
   931      \                                            -1 * q_0 / 24 * 0 ^^^ 4),\
   931      \                                            -1 * q_0 / 24 * 0 ^^^ 4),\
   932      \ 0 = c_2 + c * L + 1 / EI * (L * q_0 / 12 * L ^^^ 3 +                \
   932      \ 0 = c_2 + c * L + 1 / EI * (L * q_0 / 12 * L ^^^ 3 +                \
   933      \                                            -1 * q_0 / 24 * L ^^^ 4)]",
   933      \                                            -1 * q_0 / 24 * L ^^^ 4)]",
   934      "solveForVars [c, c_2]", "solution L"];
   934      "solveForVars [c, c_2]", "solution L"];
   935 val (dI',pI',mI') =
   935 val (dI',pI',mI') =
   936   ("Biegelinie.thy",["linear", "system"], ["no_met"]);
   936   ("Biegelinie",["linear", "system"], ["no_met"]);
   937 val p = e_pos'; val c = []; 
   937 val p = e_pos'; val c = []; 
   938 val (p,_,f,nxt,_,pt) = CalcTreeTEST [(fmz, (dI',pI',mI'))];
   938 val (p,_,f,nxt,_,pt) = CalcTreeTEST [(fmz, (dI',pI',mI'))];
   939 val (p,_,f,nxt,_,pt) = me nxt p c pt;val (p,_,f,nxt,_,pt) = me nxt p c pt;
   939 val (p,_,f,nxt,_,pt) = me nxt p c pt;val (p,_,f,nxt,_,pt) = me nxt p c pt;
   940 val (p,_,f,nxt,_,pt) = me nxt p c pt;val (p,_,f,nxt,_,pt) = me nxt p c pt;
   940 val (p,_,f,nxt,_,pt) = me nxt p c pt;val (p,_,f,nxt,_,pt) = me nxt p c pt;
   941 val (p,_,f,nxt,_,pt) = me nxt p c pt;val (p,_,f,nxt,_,pt) = me nxt p c pt;
   941 val (p,_,f,nxt,_,pt) = me nxt p c pt;val (p,_,f,nxt,_,pt) = me nxt p c pt;
   950 val (p,_,f,nxt,_,pt) = me nxt p c pt;f2str f;
   950 val (p,_,f,nxt,_,pt) = me nxt p c pt;f2str f;
   951 if f2str f = 
   951 if f2str f = 
   952 "[c_2 = 0, L * c + c_2 = -1 * q_0 * L ^^^ 4 / (24 * EI)]"
   952 "[c_2 = 0, L * c + c_2 = -1 * q_0 * L ^^^ 4 / (24 * EI)]"
   953 then () else raise error "eqsystem.sml me simpl. before SubProblem b";
   953 then () else raise error "eqsystem.sml me simpl. before SubProblem b";
   954 case nxt of
   954 case nxt of
   955     (_, Subproblem ("Biegelinie.thy", ["triangular", "2x2", "linear",_])) => ()
   955     (_, Subproblem ("Biegelinie", ["triangular", "2x2", "linear",_])) => ()
   956   | _ => raise error "eqsystem.sml me [linear,system] SubProblem b";
   956   | _ => raise error "eqsystem.sml me [linear,system] SubProblem b";
   957 
   957 
   958 val (p,_,f,nxt,_,pt) = me nxt p c pt;val (p,_,f,nxt,_,pt) = me nxt p c pt;
   958 val (p,_,f,nxt,_,pt) = me nxt p c pt;val (p,_,f,nxt,_,pt) = me nxt p c pt;
   959 val (p,_,f,nxt,_,pt) = me nxt p c pt;val (p,_,f,nxt,_,pt) = me nxt p c pt;
   959 val (p,_,f,nxt,_,pt) = me nxt p c pt;val (p,_,f,nxt,_,pt) = me nxt p c pt;
   960 val (p,_,f,nxt,_,pt) = me nxt p c pt;val (p,_,f,nxt,_,pt) = me nxt p c pt;
   960 val (p,_,f,nxt,_,pt) = me nxt p c pt;val (p,_,f,nxt,_,pt) = me nxt p c pt;
  1002 "------- Bsp 7.27";
  1002 "------- Bsp 7.27";
  1003 states:=[];
  1003 states:=[];
  1004 CalcTree [(["Traegerlaenge L","Streckenlast q_0","Biegelinie y",
  1004 CalcTree [(["Traegerlaenge L","Streckenlast q_0","Biegelinie y",
  1005 	     "Randbedingungen [y 0 = 0, y L = 0, M_b 0 = 0, M_b L = 0]",
  1005 	     "Randbedingungen [y 0 = 0, y L = 0, M_b 0 = 0, M_b L = 0]",
  1006 	     "FunktionsVariable x"],
  1006 	     "FunktionsVariable x"],
  1007 	    ("Biegelinie.thy", ["Biegelinien"],
  1007 	    ("Biegelinie", ["Biegelinien"],
  1008 		     ["IntegrierenUndKonstanteBestimmen2"]))];
  1008 		     ["IntegrierenUndKonstanteBestimmen2"]))];
  1009 moveActiveRoot 1;
  1009 moveActiveRoot 1;
  1010 (*
  1010 (*
  1011 trace_script := true; autoCalculate 1 CompleteCalc; trace_script := false;
  1011 trace_script := true; autoCalculate 1 CompleteCalc; trace_script := false;
  1012 ##7.27##          ordered           substs
  1012 ##7.27##          ordered           substs
  1043 states:=[];
  1043 states:=[];
  1044 CalcTree [(["Traegerlaenge L","Momentenlinie (-q_0 / L * x^3 / 6)",
  1044 CalcTree [(["Traegerlaenge L","Momentenlinie (-q_0 / L * x^3 / 6)",
  1045 	    "Biegelinie y",
  1045 	    "Biegelinie y",
  1046 	    "Randbedingungen [y L = 0, y' L = 0]",
  1046 	    "Randbedingungen [y L = 0, y' L = 0]",
  1047 	    "FunktionsVariable x"],
  1047 	    "FunktionsVariable x"],
  1048 	   ("Biegelinie.thy", ["vonMomentenlinieZu","Biegelinien"],
  1048 	   ("Biegelinie", ["vonMomentenlinieZu","Biegelinien"],
  1049 	    ["Biegelinien", "AusMomentenlinie"]))];
  1049 	    ["Biegelinien", "AusMomentenlinie"]))];
  1050 moveActiveRoot 1;
  1050 moveActiveRoot 1;
  1051 (*
  1051 (*
  1052 trace_script := true; autoCalculate 1 CompleteCalc; trace_script := false;
  1052 trace_script := true; autoCalculate 1 CompleteCalc; trace_script := false;
  1053 *)
  1053 *)
  1055 "------- Bsp 7.69";
  1055 "------- Bsp 7.69";
  1056 states:=[];
  1056 states:=[];
  1057 CalcTree [(["Traegerlaenge L","Streckenlast q_0","Biegelinie y",
  1057 CalcTree [(["Traegerlaenge L","Streckenlast q_0","Biegelinie y",
  1058 	     "Randbedingungen [y 0 = 0, y L = 0, y' 0 = 0, y' L = 0]",
  1058 	     "Randbedingungen [y 0 = 0, y L = 0, y' 0 = 0, y' L = 0]",
  1059 	     "FunktionsVariable x"],
  1059 	     "FunktionsVariable x"],
  1060 	    ("Biegelinie.thy", ["Biegelinien"],
  1060 	    ("Biegelinie", ["Biegelinien"],
  1061 	     ["IntegrierenUndKonstanteBestimmen2"] ))];
  1061 	     ["IntegrierenUndKonstanteBestimmen2"] ))];
  1062 moveActiveRoot 1;
  1062 moveActiveRoot 1;
  1063 (*
  1063 (*
  1064 trace_script := true; autoCalculate 1 CompleteCalc; trace_script := false;
  1064 trace_script := true; autoCalculate 1 CompleteCalc; trace_script := false;
  1065 ##7.69##          ordered           subst                   2x2
  1065 ##7.69##          ordered           subst                   2x2
  1075 "------- Bsp 7.70";
  1075 "------- Bsp 7.70";
  1076 states:=[];
  1076 states:=[];
  1077 CalcTree [(["Traegerlaenge L","Streckenlast q_0","Biegelinie y",
  1077 CalcTree [(["Traegerlaenge L","Streckenlast q_0","Biegelinie y",
  1078 	     "Randbedingungen [Q 0 = q_0 * L, M_b L = 0, y 0 = 0, y' 0 = 0]",
  1078 	     "Randbedingungen [Q 0 = q_0 * L, M_b L = 0, y 0 = 0, y' 0 = 0]",
  1079 	     "FunktionsVariable x"],
  1079 	     "FunktionsVariable x"],
  1080 	    ("Biegelinie.thy", ["Biegelinien"],
  1080 	    ("Biegelinie", ["Biegelinien"],
  1081 	     ["IntegrierenUndKonstanteBestimmen2"] ))];
  1081 	     ["IntegrierenUndKonstanteBestimmen2"] ))];
  1082 moveActiveRoot 1;
  1082 moveActiveRoot 1;
  1083 (*
  1083 (*
  1084 trace_script := true; autoCalculate 1 CompleteCalc; trace_script := false;
  1084 trace_script := true; autoCalculate 1 CompleteCalc; trace_script := false;
  1085 ##7.70##        |subst
  1085 ##7.70##        |subst
  1123 		\ 0 = (2 * c_2 + 2 * L * c + -1 * L ^^^ 2 * q_0) / 2,\
  1123 		\ 0 = (2 * c_2 + 2 * L * c + -1 * L ^^^ 2 * q_0) / 2,\
  1124 		\ 0 = c_4,                           \
  1124 		\ 0 = c_4,                           \
  1125 		\ 0 = c_3]", 
  1125 		\ 0 = c_3]", 
  1126 	    "solveForVars [c, c_2, c_3, c_4]", "solution L"];
  1126 	    "solveForVars [c, c_2, c_3, c_4]", "solution L"];
  1127 val (dI',pI',mI') =
  1127 val (dI',pI',mI') =
  1128   ("Biegelinie.thy",["linear", "system"],["no_met"]);
  1128   ("Biegelinie",["linear", "system"],["no_met"]);
  1129 val p = e_pos'; val c = []; 
  1129 val p = e_pos'; val c = []; 
  1130 val (p,_,f,nxt,_,pt) = CalcTreeTEST [(fmz, (dI',pI',mI'))];
  1130 val (p,_,f,nxt,_,pt) = CalcTreeTEST [(fmz, (dI',pI',mI'))];
  1131 val (p,_,f,nxt,_,pt) = me nxt p c pt;val (p,_,f,nxt,_,pt) = me nxt p c pt;
  1131 val (p,_,f,nxt,_,pt) = me nxt p c pt;val (p,_,f,nxt,_,pt) = me nxt p c pt;
  1132 val (p,_,f,nxt,_,pt) = me nxt p c pt;val (p,_,f,nxt,_,pt) = me nxt p c pt;
  1132 val (p,_,f,nxt,_,pt) = me nxt p c pt;val (p,_,f,nxt,_,pt) = me nxt p c pt;
  1133 val (p,_,f,nxt,_,pt) = me nxt p c pt;val (p,_,f,nxt,_,pt) = me nxt p c pt;
  1133 val (p,_,f,nxt,_,pt) = me nxt p c pt;val (p,_,f,nxt,_,pt) = me nxt p c pt;
  1301 "----- 7.70 with met top_down_: me";
  1301 "----- 7.70 with met top_down_: me";
  1302 val fmz = ["equalities                                         \
  1302 val fmz = ["equalities                                         \
  1303 	    \[c = L * q_0, L * c + c_2 = q_0 * L ^^^ 2 / 2, c_3 = 0, c_4 = 0]",
  1303 	    \[c = L * q_0, L * c + c_2 = q_0 * L ^^^ 2 / 2, c_3 = 0, c_4 = 0]",
  1304 	    "solveForVars [c, c_2, c_3, c_4]", "solution L"];
  1304 	    "solveForVars [c, c_2, c_3, c_4]", "solution L"];
  1305 val (dI',pI',mI') =
  1305 val (dI',pI',mI') =
  1306   ("Biegelinie.thy",["linear", "system"],["no_met"]);
  1306   ("Biegelinie",["linear", "system"],["no_met"]);
  1307 val p = e_pos'; val c = []; 
  1307 val p = e_pos'; val c = []; 
  1308 val (p,_,f,nxt,_,pt) = CalcTreeTEST [(fmz, (dI',pI',mI'))];
  1308 val (p,_,f,nxt,_,pt) = CalcTreeTEST [(fmz, (dI',pI',mI'))];
  1309 val (p,_,f,nxt,_,pt) = me nxt p c pt;val (p,_,f,nxt,_,pt) = me nxt p c pt;
  1309 val (p,_,f,nxt,_,pt) = me nxt p c pt;val (p,_,f,nxt,_,pt) = me nxt p c pt;
  1310 val (p,_,f,nxt,_,pt) = me nxt p c pt;val (p,_,f,nxt,_,pt) = me nxt p c pt;
  1310 val (p,_,f,nxt,_,pt) = me nxt p c pt;val (p,_,f,nxt,_,pt) = me nxt p c pt;
  1311 val (p,_,f,nxt,_,pt) = me nxt p c pt;val (p,_,f,nxt,_,pt) = me nxt p c pt;
  1311 val (p,_,f,nxt,_,pt) = me nxt p c pt;val (p,_,f,nxt,_,pt) = me nxt p c pt;
  1326 "------- Bsp 7.71";
  1326 "------- Bsp 7.71";
  1327 states:=[];
  1327 states:=[];
  1328 CalcTree [(["Traegerlaenge L","Streckenlast q_0","Biegelinie y",
  1328 CalcTree [(["Traegerlaenge L","Streckenlast q_0","Biegelinie y",
  1329 	     "Randbedingungen [M_b L = 0, y 0 = 0, y L = 0, y' 0 = 0]",
  1329 	     "Randbedingungen [M_b L = 0, y 0 = 0, y L = 0, y' 0 = 0]",
  1330 	     "FunktionsVariable x"],
  1330 	     "FunktionsVariable x"],
  1331 	    ("Biegelinie.thy", ["Biegelinien"],
  1331 	    ("Biegelinie", ["Biegelinien"],
  1332 	     ["IntegrierenUndKonstanteBestimmen2"] ))];
  1332 	     ["IntegrierenUndKonstanteBestimmen2"] ))];
  1333 moveActiveRoot 1;
  1333 moveActiveRoot 1;
  1334 (*
  1334 (*
  1335 trace_script := true; autoCalculate 1 CompleteCalc; trace_script := false;
  1335 trace_script := true; autoCalculate 1 CompleteCalc; trace_script := false;
  1336 ##7.71##       |ordered       |subst.singles (recurs) |2x2       |diagonal
  1336 ##7.71##       |ordered       |subst.singles (recurs) |2x2       |diagonal
  1348 CalcTree [(["Traegerlaenge L",
  1348 CalcTree [(["Traegerlaenge L",
  1349 	    "Momentenlinie ((q_0 * L)/ 6 * x - q_0 /(6 * L) * x^3)",
  1349 	    "Momentenlinie ((q_0 * L)/ 6 * x - q_0 /(6 * L) * x^3)",
  1350 	    "Biegelinie y",
  1350 	    "Biegelinie y",
  1351 	    "Randbedingungen [y 0 = 0, y L = 0]",
  1351 	    "Randbedingungen [y 0 = 0, y L = 0]",
  1352 	    "FunktionsVariable x"],
  1352 	    "FunktionsVariable x"],
  1353 	   ("Biegelinie.thy", ["vonMomentenlinieZu","Biegelinien"],
  1353 	   ("Biegelinie", ["vonMomentenlinieZu","Biegelinien"],
  1354 	    ["Biegelinien", "AusMomentenlinie"]))];
  1354 	    ["Biegelinien", "AusMomentenlinie"]))];
  1355 moveActiveRoot 1;
  1355 moveActiveRoot 1;
  1356 (*
  1356 (*
  1357 trace_script := true; autoCalculate 1 CompleteCalc; trace_script := false;
  1357 trace_script := true; autoCalculate 1 CompleteCalc; trace_script := false;
  1358 *)
  1358 *)
  1360 "------- Bsp 7.72b";
  1360 "------- Bsp 7.72b";
  1361 states:=[];
  1361 states:=[];
  1362 CalcTree [(["Traegerlaenge L","Streckenlast (q_0 / L * x)","Biegelinie y",
  1362 CalcTree [(["Traegerlaenge L","Streckenlast (q_0 / L * x)","Biegelinie y",
  1363 	    "Randbedingungen [M_b 0 = 0, M_b L = 0, y 0 = 0, y L = 0]",
  1363 	    "Randbedingungen [M_b 0 = 0, M_b L = 0, y 0 = 0, y L = 0]",
  1364 	    "FunktionsVariable x"],
  1364 	    "FunktionsVariable x"],
  1365 	   ("Biegelinie.thy", ["Biegelinien"],
  1365 	   ("Biegelinie", ["Biegelinien"],
  1366 	    ["IntegrierenUndKonstanteBestimmen2"] ))];
  1366 	    ["IntegrierenUndKonstanteBestimmen2"] ))];
  1367 moveActiveRoot 1;
  1367 moveActiveRoot 1;
  1368 (*
  1368 (*
  1369 trace_script := true; autoCalculate 1 CompleteCalc; trace_script := false;
  1369 trace_script := true; autoCalculate 1 CompleteCalc; trace_script := false;
  1370 ##7.72b##      |ord. |subst.singles         |ord.triang.
  1370 ##7.72b##      |ord. |subst.singles         |ord.triang.
  1381 states:=[];
  1381 states:=[];
  1382 CalcTree [(["Traegerlaenge L","Momentenlinie ???",(*description unclear*)
  1382 CalcTree [(["Traegerlaenge L","Momentenlinie ???",(*description unclear*)
  1383 	    "Biegelinie y",
  1383 	    "Biegelinie y",
  1384 	    "Randbedingungen [y L = 0, y' L = 0]",
  1384 	    "Randbedingungen [y L = 0, y' L = 0]",
  1385 	    "FunktionsVariable x"],
  1385 	    "FunktionsVariable x"],
  1386 	   ("Biegelinie.thy", ["vonMomentenlinieZu","Biegelinien"],
  1386 	   ("Biegelinie", ["vonMomentenlinieZu","Biegelinien"],
  1387 	    ["Biegelinien", "AusMomentenlinie"]))];
  1387 	    ["Biegelinien", "AusMomentenlinie"]))];
  1388 moveActiveRoot 1;
  1388 moveActiveRoot 1;
  1389 (*
  1389 (*
  1390 trace_script := true; autoCalculate 1 CompleteCalc; trace_script := false;
  1390 trace_script := true; autoCalculate 1 CompleteCalc; trace_script := false;
  1391 *)
  1391 *)
  1401 	   \ 0 = c_4 + L * c_3 +(12 * L ^^^ 2 * c_2 + 4 * L ^^^ 3 * c + -1 * L ^^^ 4 * q_0) / (-24 * EI),                                       \
  1401 	   \ 0 = c_4 + L * c_3 +(12 * L ^^^ 2 * c_2 + 4 * L ^^^ 3 * c + -1 * L ^^^ 4 * q_0) / (-24 * EI),                                       \
  1402 	   \ 0 = c_2,                                           \
  1402 	   \ 0 = c_2,                                           \
  1403 	   \ 0 = (2 * c_2 + 2 * L * c + -1 * L ^^^ 2 * q_0) / 2]", 
  1403 	   \ 0 = (2 * c_2 + 2 * L * c + -1 * L ^^^ 2 * q_0) / 2]", 
  1404 	   "solveForVars [c, c_2, c_3, c_4]", "solution L"];
  1404 	   "solveForVars [c, c_2, c_3, c_4]", "solution L"];
  1405 val (dI',pI',mI') =
  1405 val (dI',pI',mI') =
  1406   ("Biegelinie.thy",["normalize", "4x4", "linear", "system"],
  1406   ("Biegelinie",["normalize", "4x4", "linear", "system"],
  1407    ["EqSystem","normalize","4x4"]);
  1407    ["EqSystem","normalize","4x4"]);
  1408 val p = e_pos'; val c = []; 
  1408 val p = e_pos'; val c = []; 
  1409 val (p,_,f,nxt,_,pt) = CalcTreeTEST [(fmz, (dI',pI',mI'))];
  1409 val (p,_,f,nxt,_,pt) = CalcTreeTEST [(fmz, (dI',pI',mI'))];
  1410 val (p,_,f,nxt,_,pt) = me nxt p c pt; val (p,_,f,nxt,_,pt) = me nxt p c pt;
  1410 val (p,_,f,nxt,_,pt) = me nxt p c pt; val (p,_,f,nxt,_,pt) = me nxt p c pt;
  1411 val (p,_,f,nxt,_,pt) = me nxt p c pt; val (p,_,f,nxt,_,pt) = me nxt p c pt;
  1411 val (p,_,f,nxt,_,pt) = me nxt p c pt; val (p,_,f,nxt,_,pt) = me nxt p c pt;