version 1.7, 2000/12/14 01:38:37 |
version 1.8, 2000/12/14 09:13:37 |
|
|
* DEVELOPER SHALL HAVE NO LIABILITY IN CONNECTION WITH THE USE, |
* DEVELOPER SHALL HAVE NO LIABILITY IN CONNECTION WITH THE USE, |
* PERFORMANCE OR NON-PERFORMANCE OF THE SOFTWARE. |
* PERFORMANCE OR NON-PERFORMANCE OF THE SOFTWARE. |
* |
* |
* $OpenXM: OpenXM_contrib2/asir2000/lib/bfct,v 1.6 2000/12/13 05:37:31 noro Exp $ |
* $OpenXM: OpenXM_contrib2/asir2000/lib/bfct,v 1.7 2000/12/14 01:38:37 noro Exp $ |
*/ |
*/ |
/* requires 'primdec' */ |
/* requires 'primdec' */ |
|
|
|
|
def ann(F) |
def ann(F) |
{ |
{ |
V = vars(F); |
V = vars(F); |
W = append([y1,y2,t],V); |
|
N = length(V); |
N = length(V); |
B = [1-y1*y2,t-y1*F]; |
D = newvect(N); |
|
|
|
for ( I = 0; I < N; I++ ) |
|
D[I] = [deg(F,V[I]),V[I]]; |
|
qsort(D,compare_first); |
|
for ( V = [], I = N-1; I >= 0; I-- ) |
|
V = cons(D[I][1],V); |
|
|
for ( I = N-1, DV = []; I >= 0; I-- ) |
for ( I = N-1, DV = []; I >= 0; I-- ) |
DV = cons(strtov("d"+rtostr(V[I])),DV); |
DV = cons(strtov("d"+rtostr(V[I])),DV); |
|
|
|
W = append([y1,y2,t],V); |
DW = append([dy1,dy2,dt],DV); |
DW = append([dy1,dy2,dt],DV); |
|
|
|
B = [1-y1*y2,t-y1*F]; |
for ( I = 0; I < N; I++ ) { |
for ( I = 0; I < N; I++ ) { |
B = cons(DV[I]+y1*diff(F,V[I])*dt,B); |
B = cons(DV[I]+y1*diff(F,V[I])*dt,B); |
} |
} |
|
|
|
|
def indicial1(F,V) |
def indicial1(F,V) |
{ |
{ |
V = vars(F); |
|
W = append([y1,t],V); |
W = append([y1,t],V); |
N = length(V); |
N = length(V); |
B = [t-y1*F]; |
B = [t-y1*F]; |