|
neuper@42289
|
1 |
(* Title: partial fraction decomposition
|
|
neuper@42289
|
2 |
Author: Jan Rocnik
|
|
neuper@42289
|
3 |
(c) due to copyright terms
|
|
neuper@42289
|
4 |
*)
|
|
neuper@42289
|
5 |
|
|
neuper@42289
|
6 |
"--------------------------------------------------------";
|
|
neuper@42289
|
7 |
"table of contents --------------------------------------";
|
|
neuper@42289
|
8 |
"--------------------------------------------------------";
|
|
neuper@42305
|
9 |
"----------- why helpless here ? ------------------------";
|
|
neuper@42289
|
10 |
"--------------------------------------------------------";
|
|
neuper@42289
|
11 |
"--------------------------------------------------------";
|
|
neuper@42289
|
12 |
"--------------------------------------------------------";
|
|
neuper@42289
|
13 |
|
|
neuper@42289
|
14 |
|
|
neuper@42305
|
15 |
"----------- why helpless here ? ------------------------";
|
|
neuper@42305
|
16 |
"----------- why helpless here ? ------------------------";
|
|
neuper@42305
|
17 |
"----------- why helpless here ? ------------------------";
|
|
neuper@42305
|
18 |
val fmz = ["filterExpression (X z = 3 / (z - 1/4 + -1/8 * (1/(z::real))))",
|
|
neuper@42305
|
19 |
"stepResponse (x[n::real]::bool)"];
|
|
neuper@42305
|
20 |
val (dI,pI,mI) = ("Isac", ["inverse", "Z_Transform", "SignalProcessing"],
|
|
neuper@42305
|
21 |
["SignalProcessing","Z_Transform","inverse"]);
|
|
neuper@42305
|
22 |
val (p,_,f,nxt,_,pt) = CalcTreeTEST [(fmz, (dI,pI,mI))];
|
|
neuper@42305
|
23 |
val (p,_,f,nxt,_,pt) = me nxt p [] pt; "Add_Given";
|
|
neuper@42305
|
24 |
val (p,_,f,nxt,_,pt) = me nxt p [] pt; "Add_Find";
|
|
neuper@42305
|
25 |
val (p,_,f,nxt,_,pt) = me nxt p [] pt; "Specify_Theory";
|
|
neuper@42305
|
26 |
val (p,_,f,nxt,_,pt) = me nxt p [] pt; "Specify_Problem";
|
|
neuper@42305
|
27 |
val (p,_,f,nxt,_,pt) = me nxt p [] pt; "Specify_Method";
|
|
neuper@42305
|
28 |
val (p,_,f,nxt,_,pt) = me nxt p [] pt; "nxt = Apply_Method";
|
|
neuper@42305
|
29 |
val (p,_,f,nxt,_,pt) = me nxt p [] pt;
|
|
neuper@42305
|
30 |
"nxt = Rewrite (ruleZY, Inverse_Z_Transform.ruleZY) --> X z = 3 / (z - 1 / 4 + -1 / 8 * (1 / z))";
|
|
neuper@42305
|
31 |
val (p,_,f,nxt,_,pt) = me nxt p [] pt;
|
|
neuper@42305
|
32 |
"nxt = Rewrite_Set norm_Rational --> X' z = 3 / (z * (z - 1 / 4 + -1 / 8 * (1 / z)))";
|
|
neuper@42305
|
33 |
"~~~~~ fun me, args:"; val ((_,tac), (p:pos'), _, (pt:ptree)) = (nxt, p, [], pt);
|
|
neuper@42305
|
34 |
val ("ok", (_, _, ptp)) = locatetac tac (pt,p)
|
|
neuper@42305
|
35 |
val (pt, p) = ptp;
|
|
neuper@42305
|
36 |
"~~~~~ fun step, args:"; val (((ip as (_,p_)):pos'), ((ptp as (pt,p), tacis):calcstate)) =
|
|
neuper@42305
|
37 |
(p, ((pt, e_pos'),[]));
|
|
neuper@42305
|
38 |
val pIopt = get_pblID (pt,ip);
|
|
neuper@42305
|
39 |
ip = ([],Res); "false";
|
|
neuper@42305
|
40 |
tacis; " = []";
|
|
neuper@42305
|
41 |
pIopt; (* = SOME ["inverse", "Z_Transform", "SignalProcessing"]*)
|
|
neuper@42305
|
42 |
member op = [Pbl,Met] p_ andalso is_none (get_obj g_env pt (fst p)); "false";
|
|
neuper@42305
|
43 |
(*nxt_solve_ (pt,ip); "WAS isalist2list applied to NON-list 'no_meth'"
|
|
neuper@42305
|
44 |
THIS MEANS: replace no_meth by [no_meth] in Script.*)
|
|
neuper@42305
|
45 |
(*WAS val ("helpless",_) = step p ((pt, e_pos'),[]) *)
|
|
neuper@42305
|
46 |
(*WAS val (p,_,f,nxt,_,pt) = me nxt p [] pt; "Empty_Tac instead SubProblem";*)
|
|
neuper@42289
|
47 |
|