# Ticket #214: b2

File b2, 1.0 KB (added by hannes, 13 years ago)

example 2

Line
1timer = 1;
2option(mem);
3ring r = 0,(x(1),x(2),t(1),dx(1),dx(2),dt(1),u(1)),Dp;
4matrix C[7][7] = 0,1,1,1,1,1,1,0,0,1,1,1,1,1,0,0,0,1,1,1,1,0,0,0,0,1,1,1,0,0,0,0,0,1,1,0,0,0,0,0,0,1,0,0,0,0,0,0,0;
5matrix D[7][7] = 0,0,0,1,0,0,0,0,0,0,0,1,0,0,0,0,0,0,0,1,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0;
6def S = nc_algebra(C,D);
7setring S;
8basering;
9ideal I = x(1)^5*u(1)-5*x(1)^4*x(2)*u(1)+10*x(1)^3*x(2)^2*u(1)-10*x(1)^2*x(2)^3*u(1)+5*x(1)*x(2)^4*u(1)-x(2)^5*u(1)-x(1)^2*u(1)-2*x(1)*x(2)*u(1)-x(2)^2*u(1)+t(1)*u(1),
10-5*x(1)^4*dt(1)*u(1)+20*x(1)^3*x(2)*dt(1)*u(1)-30*x(1)^2*x(2)^2*dt(1)*u(1)+20*x(1)*x(2)^3*dt(1)*u(1)-5*x(2)^4*dt(1)*u(1)+2*x(1)*dt(1)*u(1)+2*x(2)*dt(1)*u(1)+dx(1)*u(1),
115*x(1)^4*dt(1)*u(1)-20*x(1)^3*x(2)*dt(1)*u(1)+30*x(1)^2*x(2)^2*dt(1)*u(1)-20*x(1)*x(2)^3*dt(1)*u(1)+5*x(2)^4*dt(1)*u(1)+2*x(1)*dt(1)*u(1)+2*x(2)*dt(1)*u(1)+dx(2)*u(1),
12-x(1)*x(2)*u(1)+x(1)*x(2);
13ideal J = eliminate(I,u(1));
14memory(0);memory(1);memory(2);
15listvar(all);
16kill J;
17memory(0);memory(1);memory(2);
18kill S,r;
19listvar(all);
20memory(0);memory(1);memory(2);
21quit;