1 (* Title: Test for rational equations
2 Author: Richard Lang 2009
3 (c) copyright due to lincense terms.
6 "-----------------------------------------------------------------";
7 "table of contents -----------------------------------------------";
8 "-----------------------------------------------------------------";
9 "----------- pbl: rational, univariate, equation ----------------";
10 "----------- solve (1/x = 5, x) by me ---------------------------";
11 "----------- S.68, Bsp.: 40, ((x)/(x - 8) + (x - 8)/(x) = 26/5)--";
12 "----------- x / (x ^ 2 - 6 * x + 9) - 1 / (x ^ 2 - 3 * x) = 1 /x";
13 "----------- ((5*x)/(x - 2) - x/(x+2)=(4::real)), incl. refine ---------------------------------";
14 "-----------------------------------------------------------------";
15 "-----------------------------------------------------------------";
17 val thy = @{theory "RatEq"};
18 val ctxt = Proof_Context.init_global thy;
20 "------------ pbl: rational, univariate, equation ----------------";
21 "------------ pbl: rational, univariate, equation ----------------";
22 "------------ pbl: rational, univariate, equation ----------------";
23 val t = TermC.parseNEW' ctxt "(1/b+1/x=1) is_ratequation_in x";
24 val SOME (t_, _) = rewrite_set_ ctxt false RatEq_prls t;
25 val result = UnparseC.term t_;
26 if result <> "True" then error "rateq.sml: new behaviour 1:" else ();
28 val t = TermC.parseNEW' ctxt "(sqrt(x)=1) is_ratequation_in x";
29 val SOME (t_, _) = rewrite_set_ ctxt false RatEq_prls t;
30 val result = UnparseC.term t_;
31 if result <> "False" then error "rateq.sml: new behaviour 2:" else ();
33 val t = TermC.parseNEW' ctxt "(x=- 1) is_ratequation_in x";
34 val SOME (t_,_) = rewrite_set_ ctxt false RatEq_prls t;
35 val result = UnparseC.term t_;
36 if result <> "False" then error "rateq.sml: new behaviour 3:" else ();
38 val t = TermC.parseNEW' ctxt "(3 + x \<up> 2 + 1/(x \<up> 2+3)=1) is_ratequation_in x";
39 val SOME (t_,_) = rewrite_set_ ctxt false RatEq_prls t;
40 val result = UnparseC.term t_;
41 if result <> "True" then error "rateq.sml: new behaviour 4:" else ();
43 val result = M_Match.match_pbl ["equality (x=(1::real))", "solveFor x", "solutions L"]
44 (Problem.from_store ["rational", "univariate", "equation"]);
45 case result of M_Match.NoMatch' _ => () | _ => error "rateq.sml: new behaviour: 5";
47 val result = M_Match.match_pbl ["equality (3 + x \<up> 2 + 1/(x \<up> 2+3)=1)", "solveFor x", "solutions L"]
48 (Problem.from_store ["rational", "univariate", "equation"]);
49 case result of M_Match.Matches' _ => () | _ => error "rateq.sml: new behaviour: 6";
51 "------------ solve (1/x = 5, x) by me ---------------------------";
52 "------------ solve (1/x = 5, x) by me ---------------------------";
53 "------------ solve (1/x = 5, x) by me ---------------------------";
54 val fmz = ["equality (1/x=(5::real))", "solveFor x", "solutions L"];
55 val (dI',pI',mI') = ("RatEq",["univariate", "equation"],["no_met"]);
56 val (p,_,f,nxt,_,pt) = CalcTreeTEST [(fmz, (dI',pI',mI'))];
57 (* val (p,_,f,nxt,_,pt) = me nxt p [1] pt;------------- now Refine_Tacitly*)
58 (* nxt = ("Model_Problem",Model_Problem ["rational", "univariate", "equation"]) *)
59 val (p,_,f,nxt,_,pt) = me nxt p [1] pt;val (p,_,f,nxt,_,pt) = me nxt p [1] pt;
60 val (p,_,f,nxt,_,pt) = me nxt p [1] pt;val (p,_,f,nxt,_,pt) = me nxt p [1] pt;
61 val (p,_,f,nxt,_,pt) = me nxt p [1] pt;val (p,_,f,nxt,_,pt) = me nxt p [1] pt;
62 val (p,_,f,nxt,_,pt) = me nxt p [1] pt;val (p,_,f,nxt,_,pt) = me nxt p [1] pt;
64 case nxt of (Rewrite_Set "RatEq_eliminate") => () | _ => ((*not checked before*));
65 val (p,_,f,nxt,_,pt) = me nxt p [1] pt;
67 WN120317.TODO dropped rateq: here "x ~= 0 should TermC.sub_at to ctxt, but it does not:
68 --- repair NO asms from rls RatEq_eliminate --- shows why.
69 so it needs more effort to find out, how Check_elementwise worked in 2002, see below.
72 (* val nxt = (_,Subproblem ("RatEq",["univariate", "equation"] ======= *)
73 val (p,_,f,nxt,_,pt) = me nxt p [1] pt;
74 (* val (p,_,f,nxt,_,pt) = me nxt p [1] pt;------------- now Refine_Tacitly*)
75 (*val nxt = ("Model_Problem", Model_Problem ["normalise", "polynomial", "univariate", "equation"])*)
76 val (p,_,f,nxt,_,pt) = me nxt p [1] pt;val (p,_,f,nxt,_,pt) = me nxt p [1] pt;
77 val (p,_,f,nxt,_,pt) = me nxt p [1] pt;val (p,_,f,nxt,_,pt) = me nxt p [1] pt;
78 val (p,_,f,nxt,_,pt) = me nxt p [1] pt;val (p,_,f,nxt,_,pt) = me nxt p [1] pt;
79 val (p,_,f,nxt,_,pt) = me nxt p [1] pt;
80 (*val nxt = Apply_Method ["PolyEq", "normalise_poly"])*)
81 val (p,_,f,nxt,_,pt) = me nxt p [1] pt;
82 val (p,_,f,nxt,_,pt) = me nxt p [1] pt;val (p,_,f,nxt,_,pt) = me nxt p [1] pt;
83 (* val nxt = (_,Subproblem ("PolyEq",["polynomial", "univariate", "equation"]=======*)
84 val (p,_,f,nxt,_,pt) = me nxt p [1] pt;
85 val (p,_,f,nxt,_,pt) = me nxt p [1] pt;
86 (* ("Model_Problem", Model_Problem ["degree_1", "polynomial", "univariate", "equation"])*)
87 val (p,_,f,nxt,_,pt) = me nxt p [1] pt;val (p,_,f,nxt,_,pt) = me nxt p [1] pt;
88 val (p,_,f,nxt,_,pt) = me nxt p [1] pt;val (p,_,f,nxt,_,pt) = me nxt p [1] pt;
89 val (p,_,f,nxt,_,pt) = me nxt p [1] pt;val (p,_,f,nxt,_,pt) = me nxt p [1] pt;
90 val (p,_,f,nxt,_,pt) = me nxt p [1] pt;val (p,_,f,nxt,_,pt) = me nxt p [1] pt;
91 val (p,_,f,nxt,_,pt) = me nxt p [1] pt;val (p,_,f,nxt,_,pt) = me nxt p [1] pt;
92 val (p,_,f,nxt,_,pt) = me nxt p [1] pt;val (p,_,f,nxt,_,pt) = me nxt p [1] pt;
93 val (p''',_,f,nxt''',_,pt''') = me nxt p [1] pt;
94 f2str f = "[x = 1 / 5]";
95 case nxt of (Check_elementwise "Assumptions") => () | _ => ((*not checked before*));
96 "~~~~~ fun me, args:"; val (tac, (p:pos'), _, (pt:ctree)) = (nxt, p, c, pt);
97 val (pt, p) = case Step.by_tactic tac (pt,p) of
98 ("ok", (_, _, ptp)) => ptp | _ => error "--- solve (1/x = 5.. Step.by_tactic";
99 "~~~~~ fun Step.do_next, args:"; val (ip as (_,p_), (ptp as (pt,p), tacis)) = (p, ((pt, e_pos'), []))
100 val pIopt = get_pblID (pt,ip); (*= SOME ["rational", "univariate", "equation"]
101 1- 1 associated to metID ["RatEq", "solve_rat_equation"]*)
103 member op = [Pbl,Met] p_; (*= false*)
104 "~~~~~ fun do_next, args:"; val (ptp as (pt, pos as (p, p_))) = (pt, ip);
105 val thy' = get_obj g_domID pt (par_pblobj pt p);
106 val (is, sc) = resume_prog thy' (p,p_) pt; (*is: which ctxt?*)
107 "~~~~~ fun find_next_step, args:"; val () = ();
108 (*----- outcommented during cleanup of args in lucas-interpreter.sml ------------------------\\* )
109 "~~~~~ fun go_scan_up, args:"; val (thy, ptp, (Prog sc), E, l, ay, a, v) =
110 (thy, ptp, sc, E, l, true, a, v);
111 1 < length l; (*true*)
112 val up = drop_last l;
113 TermC.sub_at up sc; (* = Const (\<^const_name>\<open>Let\<close>, *)
114 "~~~~~ fun scan_up, args:"; val (thy, ptp, (scr as (Prog sc)), E, l, ay,
115 (t as Const (\<^const_name>\<open>Let\<close>,_) $ _), a, v) = (thy, ptp, (Prog sc), E, up, ay, (TermC.sub_at up sc), a, v);
116 ay = Napp_; (*false*)
117 val up = drop_last l;
118 val (Const (\<^const_name>\<open>Let\<close>,_) $ e $ (Abs (i,T,body))) = TermC.sub_at up sc; (*Const (\<^const_name>\<open>SubProblem\<close>,..*)
119 val i = mk_Free (i, T);
120 val E = Env.update E (i, v);
121 "~~~~~ fun scan_dn, args:"; val ((thy as (th,sr)), (pt, p), E, l, t, a, v) =
122 (thy, ptp, E, (up@[R,D]), body, a, v);
123 "~~~~~ fun check_leaf, args:"; val (call, thy, srls, (E, (a, v)), t) = ("next ", th, sr, (E, (a, v)), t);
124 "~~~~~ fun eval_leaf, args:"; val (E, a, v,
125 (t as (Const (\<^const_name>\<open>Check_elementwise\<close>,_) $ _ $ _ ))) = (E, a, v, t);
126 val Program.Tac tm = Program.Tac (subst_atomic E t);
127 UnparseC.term tm = "Check_elementwise [x = 1 / 5] {v_v. Assumptions}";
128 (* ------ \<up> ----- ? "x" ?*)
129 "~~~~~ to check_leaf return val:"; val ((Program.Tac stac, a')) = ((Program.Tac (subst_atomic E t), NONE));
130 val stac' = eval_prog_expr (ThyC.get_theory thy) srls (subst_atomic (upd_env_opt E (a,v)) stac);
131 UnparseC.term stac' = "Check_elementwise [x = 1 / 5] {v_v. Assumptions}";
132 "~~~~~ to scan_dn return val:"; val ((a', Program.Tac stac)) = ((a', Program.Tac stac'));
133 val m = LItool.tac_from_prog pt (ThyC.get_theory th) stac;
134 case m of Check_elementwise "Assumptions" => () | _ => (); (*m' = Empty_Tac_ ???!??? *);
135 val (p''''', pt''''', m''''') = (p, pt, m);
136 "~~~~~ fun check , args:"; val ((p,p_), pt, (m as Check_elementwise pred)) = (p, pt, m);
137 member op = [Pbl,Met] p_; (* = false*)
138 val pp = par_pblobj pt p;
139 val thy' = (get_obj g_domID pt pp):theory';
140 val thy = ThyC.get_theory thy'
141 val metID = (get_obj g_metID pt pp)
142 val {crls,...} = MethodC.from_store metID
143 val (f,asm) = case p_ of Frm => (get_obj g_form pt p , [])
144 | Res => get_obj g_result pt p;
145 UnparseC.term f = "[x = 1 / 5]"; (*the current formula*)
146 val vp = (Proof_Context.init_global thy, pred) |-> parseNEW |> the |> mk_set thy pt p f;
147 val (bdv, asms) = vp;
149 UnparseC.term bdv = "x";
150 UnparseC.terms asms = (* asms from rewriting are missing : vvv *)
151 ("[\"~ (matches (?a = 0) (1 = 5 * x)) | ~ lhs (1 = 5 * x) is_poly_in x\",\"x = 1 / 5\", " ^
152 "\"lhs (1 + -5 * x = 0) is_poly_in x\",\"lhs (1 + -5 * x = 0) has_degree_in x = 1\", " ^
153 "\"1 / x = 5 is_ratequation_in x\"]");
155 WN120317.TODO dropped rateq: ctxt should contain "x ~= 0 here, but it does not, see above.
158 val Applicable.Yes (Check_elementwise' (curr_form, pred, (res, asms))) = check p''''' pt''''' m''''';
159 UnparseC.term curr_form = "[x = 1 / 5]";
160 pred = "Assumptions";
161 res = TermC.str2term "[]::bool list";
164 val (p,_,f,nxt,_,pt) = me nxt''' p''' [] pt'''; (*<<<----- this caused the error*)
166 val (p,_,f,nxt,_,pt) = me nxt p [1] pt;val (p,_,f,nxt,_,pt) = me nxt p [1] pt;
168 if p = ([], Res) andalso f2str f = "[x = 1 / 5]"
169 then case nxt of ("End_Proof'", End_Proof') => ()
170 | _ => error "rateq.sml: new behaviour: [x = 1 / 5] 1"
171 else error "rateq.sml: new behaviour: [x = 1 / 5] 2";
172 ( *----- outcommented during cleanup of args in lucas-interpreter.sml ------------------------//*)
174 "------------ S.68, Bsp.: 40, ((x)/(x - 8) + (x - 8)/(x) = 26/5)--";
175 "------------ S.68, Bsp.: 40, ((x)/(x - 8) + (x - 8)/(x) = 26/5)--";
176 "------------ S.68, Bsp.: 40, ((x)/(x - 8) + (x - 8)/(x) = 26/5)--";
177 (*EP Schalk_II_p68_n40*)
178 val fmz = ["equality ((x)/(x - 8) + (x - 8)/(x) = 26/(5::real))", "solveFor x", "solutions L"];
179 val (dI',pI',mI') = ("RatEq",["univariate", "equation"],["no_met"]);
180 val (p,_,f,nxt,_,pt) = CalcTreeTEST [(fmz, (dI',pI',mI'))];
181 val (p,_,f,nxt,_,pt) = me nxt p [1] pt;
182 (* nxt = ("Model_Problem",Model_Problem ["rational", "univariate", "equation"])*)
183 val (p,_,f,nxt,_,pt) = me nxt p [1] pt;val (p,_,f,nxt,_,pt) = me nxt p [1] pt;
184 val (p,_,f,nxt,_,pt) = me nxt p [1] pt;val (p,_,f,nxt,_,pt) = me nxt p [1] pt;
185 val (p,_,f,nxt,_,pt) = me nxt p [1] pt;val (p,_,f,nxt,_,pt) = me nxt p [1] pt;
186 val (p,_,f,nxt,_,pt) = me nxt p [1] pt;val (p,_,f,nxt,_,pt) = me nxt p [1] pt;
187 val (p,_,f,nxt,_,pt) = me nxt p [1] pt;val (p,_,f,nxt,_,pt) = me nxt p [1] pt;
188 val (p,_,f,nxt,_,pt) = me nxt p [1] pt;val (p,_,f,nxt,_,pt) = me nxt p [1] pt;
189 (* nxt = ("Subproblem",Subproblem ("RatEq",["univariate", "equation"]))*)
190 val (p,_,f,nxt,_,pt) = me nxt p [1] pt;val (p,_,f,nxt,_,pt) = me nxt p [1] pt;
191 (* nxt = ("Model_Problem", Model_Problem ["normalise", "polynomial", "univariate", "equation"])*)
192 val (p,_,f,nxt,_,pt) = me nxt p [1] pt;val (p,_,f,nxt,_,pt) = me nxt p [1] pt;
193 val (p,_,f,nxt,_,pt) = me nxt p [1] pt;val (p,_,f,nxt,_,pt) = me nxt p [1] pt;
194 val (p,_,f,nxt,_,pt) = me nxt p [1] pt;val (p,_,f,nxt,_,pt) = me nxt p [1] pt;
195 val (p,_,f,nxt,_,pt) = me nxt p [1] pt;val (p,_,f,nxt,_,pt) = me nxt p [1] pt;
196 val (p,_,f,nxt,_,pt) = me nxt p [1] pt;val (p,_,f,nxt,_,pt) = me nxt p [1] pt;
198 if p = ([4, 3], Pbl) then ()
201 (Add_Given "equality (320 + 128 * x + - 16 * x \<up> 2 = 0)") =>
203 Test_Out.PpcKF (Test_Out.Problem [], {
204 Find = [Incompl "solutions []"], Given = [Incompl "equality", Incompl "solveFor"],
205 Relate = [], Where = [False "matches (?a + ?v_ \<up> 2 = 0) e_e \<or>\nmatches (?a + ?b * ?v_ \<up> 2 = 0) e_e"],
207 | _ => error ("S.68, Bsp.: 40 PblObj changed"))
208 | _ => error ("S.68, Bsp.: 40 changed nxt =" ^ Tactic.input_to_string (nxt)));
210 val (p,_,f,nxt,_,pt) = me nxt p [1] pt;
211 (* ("Subproblem", Subproblem ("PolyEq",["polynomial", "univariate", "equation"])) *)
212 val (p,_,f,nxt,_,pt) = me nxt p [1] pt;val (p,_,f,nxt,_,pt) = me nxt p [1] pt;
213 (* nxt = ("Model_Problem", Model_Problem
214 ["abcFormula", "degree_2", "polynomial", "univariate", "equation"])*)
215 val (p,_,f,nxt,_,pt) = me nxt p [1] pt;val (p,_,f,nxt,_,pt) = me nxt p [1] pt;
216 val (p,_,f,nxt,_,pt) = me nxt p [1] pt;val (p,_,f,nxt,_,pt) = me nxt p [1] pt;
217 val (p,_,f,nxt,_,pt) = me nxt p [1] pt;val (p,_,f,nxt,_,pt) = me nxt p [1] pt;
218 val (p,_,f,nxt,_,pt) = me nxt p [1] pt;val (p,_,f,nxt,_,pt) = me nxt p [1] pt;
219 val (p,_,f,nxt,_,pt) = me nxt p [1] pt;val (p,_,f,nxt,_,pt) = me nxt p [1] pt;
220 val (p,_,f,nxt,_,pt) = me nxt p [1] pt;val (p,_,f,nxt,_,pt) = me nxt p [1] pt;
221 val (p,_,f,nxt,_,pt) = me nxt p [1] pt;val (p,_,f,nxt,_,pt) = me nxt p [1] pt;
222 if p = ([], Res) andalso f2str f = "[x = - 2, x = 10]" then ()
223 else error "rateq.sml: new behaviour: [x = - 2, x = 10]";
225 "----------- remove x = 0 from [x = 0, x = 6 / 5] ----------------------------------------------";
226 "----------- remove x = 0 from [x = 0, x = 6 / 5] ----------------------------------------------";
227 "----------- remove x = 0 from [x = 0, x = 6 / 5] ----------------------------------------------";
228 (*ER-7*) (*Schalk I s.87 Bsp 55b*)
229 val fmz = ["equality (x/(x \<up> 2 - 6*x+9) - 1/(x \<up> 2 - 3*x) =1/x)",
230 "solveFor x", "solutions L"];
231 val spec = ("RatEq",["univariate", "equation"],["no_met"]);
232 val (p,_,f,nxt,_,pt) = CalcTreeTEST [(fmz, spec)]; (* 0. specify-phase *)
233 val (p,_,f,nxt,_,pt) = me nxt p [] pt; val (p,_,f,nxt,_,pt) = me nxt p [] pt;
234 val (p,_,f,nxt,_,pt) = me nxt p [] pt; val (p,_,f,nxt,_,pt) = me nxt p [] pt;
235 val (p,_,f,nxt,_,pt) = me nxt p [] pt; val (p,_,f,nxt,_,pt) = me nxt p [] pt;
237 (*[], Met*)val (p,_,f,nxt,_,pt) = me nxt p [] pt;
238 (*+*)case nxt of Apply_Method ["RatEq", "solve_rat_equation"] => ()
239 (*+*)| _ => error "55b root specification broken";
241 val (p,_,f,nxt,_,pt) = me nxt p [] pt; (* 0. solve-phase*)
242 val (p,_,f,nxt,_,pt) = me nxt p [] pt; val (p,_,f,nxt,_,pt) = me nxt p [] pt;
243 val (p,_,f,nxt,_,pt) = me nxt p [] pt; f2str f = "(3 + - 1 * x + x \<up> 2) * x = 1 * (9 * x + -6 * x \<up> 2 + x \<up> 3)";
245 (*+*)if eq_set op = (Ctree.get_assumptions pt p |> map UnparseC.term,
246 (*+*) ["x \<noteq> 0", "9 * x + - 6 * x \<up> 2 + x \<up> 3 \<noteq> 0", "x / (x \<up> 2 - 6 * x + 9) - 1 / (x \<up> 2 - 3 * x) =\n1 / x is_ratequation_in x"])
247 (*+*)then () else error "assumptions before 1. Subproblem CHANGED";
248 (*+*)if p = ([3], Res) andalso f2str f = "(3 + - 1 * x + x \<up> 2) * x = 1 * (9 * x + - 6 * x \<up> 2 + x \<up> 3)"
250 (*+*) ((case nxt of Subproblem ("PolyEq", ["normalise", "polynomial", "univariate", "equation"]) => ()
251 (*+*) | _ => error ("S.68, Bsp.: 40 nxt =" ^ Tactic.input_to_string nxt)))
252 (*+*)else error "1. Subproblem -- call changed";
254 val (p,_,f,nxt,_,pt) = me nxt p [] pt; (* 1. specify-phase *)
255 val (p,_,f,nxt,_,pt) = me nxt p [] pt; val (p,_,f,nxt,_,pt) = me nxt p [] pt;
256 val (p,_,f,nxt,_,pt) = me nxt p [] pt; val (p,_,f,nxt,_,pt) = me nxt p [] pt;
257 val (p,_,f,nxt,_,pt) = me nxt p [] pt; val (p,_,f,nxt,_,pt) = me nxt p [] pt;
259 (*[4], Met*)val (p,_,f,nxt,_,pt) = me nxt p [] pt;
260 case nxt of Apply_Method ["PolyEq", "normalise_poly"] => ()
261 | _ => error "55b normalise_poly specification broken 1";
263 val (p,_,f,nxt,_,pt) = me nxt p [] pt; (* 1. solve-phase *)
264 val (p,_,f,nxt,_,pt) = me nxt p [] pt; val (p,_,f,nxt,_,pt) = me nxt p [] pt;
265 val (p,_,f,nxt,_,pt) = me nxt p [] pt; f2str f = "-6 * x + 5 * x \<up> 2 = 0";
267 if p = ([4, 3], Res) andalso f2str f = "- 6 * x + 5 * x \<up> 2 = 0"
269 ((case nxt of Subproblem ("PolyEq", ["bdv_only", "degree_2", "polynomial", "univariate", "equation"]) => ()
270 | _ => error ("S.68, Bsp.: 40 nxt =" ^ Tactic.input_to_string nxt)))
272 val (p,_,f,nxt,_,pt) = me nxt p [] pt; (* 2. specify-phase *)
273 val (p,_,f,nxt,_,pt) = me nxt p [] pt; val (p,_,f,nxt,_,pt) = me nxt p [] pt;
274 val (p,_,f,nxt,_,pt) = me nxt p [] pt; val (p,_,f,nxt,_,pt) = me nxt p [] pt;
275 val (p,_,f,nxt,_,pt) = me nxt p [] pt; val (p,_,f,nxt,_,pt) = me nxt p [] pt;
277 (*[4, 4], Met*)val (p,_,f,nxt,_,pt) = me nxt p [] pt;(*\<rightarrow>*)
278 case nxt of Apply_Method ["PolyEq", "solve_d2_polyeq_bdvonly_equation"] => ()
279 | _ => error "55b normalise_poly specification broken 2";
281 val (p,_,f,nxt,_,pt) = me nxt p [] pt; (*f = "-6 * x + 5 * x \<up> 2 = 0"*) (* 2. solve-phase *)
282 val (p,_,f,nxt,_,pt) = me nxt p [] pt; val (p,_,f,nxt,_,pt) = me nxt p [] pt;
284 (*[4, 4, 3], Res*)val (p,_,f,nxt,_,pt) = me nxt p [] pt;(*\<rightarrow>Or_to_List*)
285 (*[4, 4, 4], Res*)val (p,_,f,nxt,_,pt) = me nxt p [] pt; f2str f = "[x = 0, x = 6 / 5]";
286 (*[4, 4, 5], Res*)val (p,_,f,nxt,_,pt) = me nxt p [] pt;(*\<rightarrow>2. Check_Postcond ["bdv_only", "degree_2", "polynomial", "univariate", "equation"]*)
288 (* *)if eq_set op = ((Ctree.get_assumptions pt p |> map UnparseC.term), [
289 (*0.pre*) "x / (x \<up> 2 - 6 * x + 9) - 1 / (x \<up> 2 - 3 * x) =\n1 / x is_ratequation_in x",
290 (*1.pre*) "\<not> matches (?a = 0)\n ((3 + - 1 * x + x \<up> 2) * x =\n 1 * (9 * x + - 6 * x \<up> 2 + x \<up> 3)) \<or>\n"
291 (*1.pre*) ^ "\<not> lhs ((3 + - 1 * x + x \<up> 2) * x =\n 1 * (9 * x + - 6 * x \<up> 2 + x \<up> 3)) is_poly_in x",
292 (*2.pre*) "lhs (- 6 * x + 5 * x \<up> 2 = 0) is_poly_in x",
293 (*2.pre*) "lhs (- 6 * x + 5 * x \<up> 2 = 0) has_degree_in x = 2",
294 (*0.asm*) "x \<noteq> 0",
295 (*0.asm*) "9 * x + - 6 * x \<up> 2 + x \<up> 3 \<noteq> 0"
297 (* *)then () else error "assumptions at end 2. Subproblem CHANGED";
299 (*[4, 4], Res*)val (p''''',_,f,nxt''''',_,pt''''') = me nxt p [] pt;(*\<rightarrow>1. Check_Postcond ["normalise", "polynomial", "univariate", "equation"]*)
301 (*/--------- step into 2. Check_Postcond SEE .. ----------------------------------------------\*)
302 "----------- rat-equ: remove x = 0 from [x = 0, x = 6 / 5] due to contexts ---------------------";
303 (*\--------- step into 2. Check_Postcond -----------------------------------------------------/*)
305 (*[4], Res*)val (p,_,f,nxt,_,pt) = me nxt''''' p''''' [] pt''''';(*\<rightarrow>IDLE LEGACY: Check_elementwise "Assumptions"*)
306 (*[], Res*)val (p,_,f,nxt,_,pt) = me nxt''''' p''''' [] pt''''';(*\<rightarrow>End_Proof'*)
308 (*/-------- final test -----------------------------------------------------------------------\*)
309 if f2str f = "[x = 6 / 5]" andalso eq_set op = (map UnparseC.term (Ctree.get_assumptions pt p),
311 "lhs (- 6 * x + 5 * x \<up> 2 = 0) is_poly_in x",
312 "lhs (- 6 * x + 5 * x \<up> 2 = 0) has_degree_in x = 2",
313 "\<not> matches (?a = 0)\n ((3 + - 1 * x + x \<up> 2) * x =\n 1 * (9 * x + - 6 * x \<up> 2 + x \<up> 3)) \<or>\n\<not> lhs ((3 + - 1 * x + x \<up> 2) * x =\n 1 * (9 * x + - 6 * x \<up> 2 + x \<up> 3)) is_poly_in x",
315 "9 * x + - 6 * x \<up> 2 + x \<up> 3 \<noteq> 0",
316 "x / (x \<up> 2 - 6 * x + 9) - 1 / (x \<up> 2 - 3 * x) =\n1 / x is_ratequation_in x"]
317 ) then () else error "test CHANGED";
320 "----------- ((5*x)/(x - 2) - x/(x+2)=(4::real)), incl. refine ---------------------------------";
321 "----------- ((5*x)/(x - 2) - x/(x+2)=(4::real)), incl. refine ---------------------------------";
322 "----------- ((5*x)/(x - 2) - x/(x+2)=(4::real)), incl. refine ---------------------------------";
323 (*was in test/../usecases.sml*)
324 val fmz = ["equality ((5*x)/(x - 2) - x/(x+2)=(4::real))", "solveFor x", "solutions L"];
325 val (dI',pI',mI') = ("RatEq", ["univariate", "equation"], ["no_met"]);
326 val (p,_,f,nxt,_,pt) = CalcTreeTEST [(fmz, (dI',pI',mI'))];
327 (*[], Pbl*)val (p,_,f,nxt,_,pt) = me nxt p [1] pt; val (p,_,f,nxt,_,pt) = me nxt p [1] pt;
328 (**)val (p,_,f,nxt,_,pt) = me nxt p [1] pt; val (p,_,f,nxt,_,pt) = me nxt p [] pt;(**)
329 (**)val (p,_,f,nxt,_,pt) = me nxt p [1] pt; val (p,_,f,nxt,_,pt) = me nxt p [] pt;(**)
330 (*[], Met*)val (p,_,f,nxt,_,pt) = me nxt p [1] pt;(*Apply_Method ["RatEq", "solve_rat_equation"]*)
331 (*[1], Frm*)val (p,_,f,nxt,_,pt) = me nxt p [] pt;(*Rewrite_Set "RatEq_simplify":*)
333 (*+*)if (get_istate_LI pt p |> Istate.string_of) (* still specify-phase: found_accept = false ---------------------------------> vvvvv*)
334 (*+*) = "Pstate ([\"\n(e_e, 5 * x / (x - 2) - x / (x + 2) = 4)\", \"\n(v_v, x)\"], [], empty, NONE, \n??.empty, ORundef, false, true)"
335 (*+*)then () else error "rat-eq + subpbl: istate in specify-phase";
337 (*[1], Res*)val (p,_,f,nxt,_,pt) = me nxt p [1] pt;(*Rewrite_Set "norm_Rational"*)
339 (*+*)if (get_istate_LI pt p |> Istate.string_of) (* solve-phase: found_accept = true -----------------------------------------------------------------------------------------------> vvvvv*)
340 (*+*) = "Pstate ([\"\n(e_e, 5 * x / (x - 2) - x / (x + 2) = 4)\", \"\n(v_v, x)\"], [R, L, R, L, L, R, R, R], empty, SOME e_e, \n5 * x / (x + - 1 * 2) + - 1 * x / (x + 2) = 4, ORundef, true, true)"
341 (*+*)then () else error "rat-eq + subpbl: istate after found_accept";
343 (*[2], Res*)val (p,_,f,nxt,_,pt) = me nxt p [] pt;(*Rewrite_Set "RatEq_eliminate"*)
346 (*[3], Res*)val (p,_,f,nxt,_,pt) = me nxt p [1] pt;(*Subproblem ("PolyEq", ["normalise", "polynomial", "univariate", "equation"])*)
347 (*[4], Pbl*)val (p,_,f,nxt,_,pt) = me nxt p [] pt;(*Model_Problem*)
349 (**)val (p,_,f,nxt,_,pt) = me nxt p [1] pt; val (p,_,f,nxt,_,pt) = me nxt p [] pt;(**)
350 (**)val (p,_,f,nxt,_,pt) = me nxt p [1] pt; val (p,_,f,nxt,_,pt) = me nxt p [] pt;(**)
351 (**)val (p,_,f,nxt,_,pt) = me nxt p [1] pt; val (p,_,f,nxt,_,pt) = me nxt p [] pt;(**)
352 (**)val (p,_,f,nxt,_,pt) = me nxt p [1] pt; val (p,_,f,nxt,_,pt) = me nxt p [] pt;(**)
354 (**)val (p,_,f,nxt,_,pt) = me nxt p [1] pt; val (p,_,f,nxt,_,pt) = me nxt p [] pt;(**)
355 (**)val (p,_,f,nxt,_,pt) = me nxt p [1] pt; val (p,_,f,nxt,_,pt) = me nxt p [] pt;(**)
356 (**)val (p,_,f,nxt,_,pt) = me nxt p [1] pt; val (p,_,f,nxt,_,pt) = me nxt p [] pt;(**)
357 (**)val (p,_,f,nxt,_,pt) = me nxt p [1] pt; val (p,_,f,nxt,_,pt) = me nxt p [] pt;(**)
359 (*[4,4], Pbl*)val (p,_,f,nxt,_,pt) = me nxt p [1] pt;(*Specify_Problem ["degree_1", "polynomial", "univariate", "equation"]*)
360 (*[4,4], Pbl*)val (p,_,f,nxt,_,pt) = me nxt p [1] pt;(*Specify_Method ["PolyEq", "solve_d1_polyeq_equation"]*)
362 (*[4,4], Met*)val (p,_,f,nxt,_,pt) = me nxt p [1] pt;(*Apply_Method ["PolyEq", "solve_d1_polyeq_equation"]*)
363 (*[4,4,1], Frm*)val (p,_,f,nxt,_,pt) = me nxt p [1] pt;(*Rewrite_Set_Inst (["(''bdv'', x)"], "d1_polyeq_simplify")*)
364 (*[4,4,1], Res*)val (p,_,f,nxt,_,pt) = me nxt p [1] pt;(*Rewrite_Set "polyeq_simplify"*)
365 (*[4,4,2], Res*)val (p,_,f,nxt,_,pt) = me nxt p [] pt;(*Or_to_List*)
366 (* f = Test_Out.FormKF "[x = -4 / 3]" *)
367 (*[4,4,3], Res*)val (p,_,f,nxt,_,pt) = me nxt p [] pt;(*Check_elementwise "Assumptions"*)
368 (*[4,4,4], Res*)val (p,_,f,nxt,_,pt) = me nxt p [] pt;(*Check_Postcond ["degree_1", "polynomial", "univariate", "equation"]*)
369 (*[4,4], Res*)val (p,_,f,nxt,_,pt) = me nxt p [] pt;(*Check_Postcond ["normalise", "polynomial", "univariate", "equation"]*)
370 (*[4], Res*)val (p,_,f,nxt,_,pt) = me nxt p [] pt;(*Check_elementwise "Assumptions"*)
371 (*[5], Res*)val (p,_,f,nxt,_,pt) = me nxt p [] pt;(*Check_Postcond ["rational", "univariate", "equation"]*)
372 (*[], Res*)val (p,_,f,nxt,_,pt) = me nxt p [] pt;(*End_Proof'*)
374 if p = ([], Res) andalso f2str f = "[x = - 4 / 3]"
375 then case nxt of End_Proof' => () | _ => error "rat-eq + subpbl: end CHANGED 1"
376 else error "rat-eq + subpbl: end CHANGED 2";