Hi all,

I have this system

> system1D := H = alpha*gamma[2, 2]*d[2, 1]-beta*d[1, 2]*gamma[1, 2]^2-gamma*d[1, 2]*gamma[2, 1]^2+alpha*gamma[2, 2]^2*d[2, 2]-beta*d[2, 2]*gamma[2, 2]^2-gamma*d[2, 2]*gamma[2, 2]^2, E = alpha*gamma[2, 1]*d[1, 1]-beta*d[1, 2]*gamma[1, 1]-gamma*d[1, 1]*gamma[2, 1]+alpha*gamma[2, 1]^2*d[1, 2]-beta*d[2, 2]*gamma[2, 1]-gamma*d[2, 1]*gamma[2, 2], B = alpha*gamma[1, 1]*d[2, 1]-beta*d[1, 1]*gamma[1, 1]^2-gamma*d[1, 1]*gamma[1, 1]^2+alpha*gamma[1, 1]^2*d[2, 2]-beta*d[2, 1]*gamma[2, 1]^2-gamma*d[2, 1]*gamma[1, 2]^2, D = alpha*gamma[1, 2]*d[2, 1]-beta*d[1, 1]*gamma[1, 2]^2-gamma*d[1, 2]*gamma[1, 1]^2+alpha*gamma[1, 2]^2*d[2, 2]-beta*d[2, 1]*gamma[2, 2]^2-gamma*d[2, 2]*gamma[1, 2]^2, A = alpha*gamma[1, 1]*d[1, 1]-beta*d[1, 1]*gamma[1, 1]-gamma*d[1, 1]*gamma[1, 1]+alpha*gamma[1, 1]^2*d[1, 2]-beta*d[2, 1]*gamma[2, 1]-gamma*d[2, 1]*gamma[1, 2], C = alpha*gamma[1, 2]*d[1, 1]-beta*d[1, 1]*gamma[1, 2]-gamma*d[1, 2]*gamma[1, 1]+alpha*gamma[1, 2]^2*d[1, 2]-beta*d[2, 1]*gamma[2, 2]-gamma*d[2, 2]*gamma[1, 2], F = alpha*gamma[2, 1]*d[2, 1]-beta*d[1, 2]*gamma[1, 1]^2-gamma*d[1, 1]*gamma[2, 1]^2+alpha*gamma[2, 1]^2*d[2, 2]-beta*d[2, 2]*gamma[2, 1]^2-gamma*d[2, 1]*gamma[2, 2]^2, G = alpha*gamma[2, 2]*d[1, 1]-beta*d[1, 2]*gamma[1, 2]-gamma*d[1, 2]*gamma[2, 1]+alpha*gamma[2, 2]^2*d[1, 2]-beta*d[2, 2]*gamma[2, 2]-gamma*d[2, 2]*gamma[2, 2], H = alpha*delta[2, 2]*d[2, 1]-beta*d[1, 2]*delta[1, 2]^2-gamma*d[1, 2]*delta[2, 1]^2+alpha*delta[2, 2]^2*d[2, 2]-beta*d[2, 2]*delta[2, 2]^2-gamma*d[2, 2]*delta[2, 2]^2, E = alpha*delta[2, 1]*d[1, 1]-beta*d[1, 2]*delta[1, 1]-gamma*d[1, 1]*delta[2, 1]+alpha*delta[2, 1]^2*d[1, 2]-beta*d[2, 2]*delta[2, 1]-gamma*d[2, 1]*delta[2, 2], B = alpha*delta[1, 1]*d[2, 1]-beta*d[1, 1]*delta[1, 1]^2-gamma*d[1, 1]*delta[1, 1]^2+alpha*delta[1, 1]^2*d[2, 2]-beta*d[2, 1]*delta[2, 1]^2-gamma*d[2, 1]*delta[1, 2]^2, D = alpha*delta[1, 2]*d[2, 1]-beta*d[1, 1]*delta[1, 2]^2-gamma*d[1, 2]*delta[1, 1]^2+alpha*delta[1, 2]^2*d[2, 2]-beta*d[2, 1]*delta[2, 2]^2-gamma*d[2, 2]*delta[1, 2]^2, A = alpha*delta[1, 1]*d[1, 1]-beta*d[1, 1]*delta[1, 1]-gamma*d[1, 1]*delta[1, 1]+alpha*delta[1, 1]^2*d[1, 2]-beta*d[2, 1]*delta[2, 1]-gamma*d[2, 1]*delta[1, 2], C = alpha*delta[1, 2]*d[1, 1]-beta*d[1, 1]*delta[1, 2]-gamma*d[1, 2]*delta[1, 1]+alpha*delta[1, 2]^2*d[1, 2]-beta*d[2, 1]*delta[2, 2]-gamma*d[2, 2]*delta[1, 2], F = alpha*delta[2, 1]*d[2, 1]-beta*d[1, 2]*delta[1, 1]^2-gamma*d[1, 1]*delta[2, 1]^2+alpha*delta[2, 1]^2*d[2, 2]-beta*d[2, 2]*delta[2, 1]^2-gamma*d[2, 1]*delta[2, 2]^2, G = alpha*delta[2, 2]*d[1, 1]-beta*d[1, 2]*delta[1, 2]-gamma*d[1, 2]*delta[2, 1]+alpha*delta[2, 2]^2*d[1, 2]-beta*d[2, 2]*delta[2, 2]-gamma*d[2, 2]*delta[2, 2];

> subs({A = 0, B = 0, C = 0, D = 0, E = 0, F = 0, G = 0, H = 0, delta[1, 1] = 1, delta[1, 2] = 0, delta[2, 1] = 0, delta[2, 2] = 0, gamma[1, 1] = 1, gamma[1, 2] = 0, gamma[2, 1] = 0, gamma[2, 2] = 0, delta[1, 1]^2 = 0, delta[1, 2]^2 = 0, delta[2, 1]^2 = 1, delta[2, 2]^2 = 0, gamma[1, 1]^2 = 0, gamma[1, 2]^2 = 1, gamma[2, 1]^2 = 0, gamma[2, 2]^2 = 0}, {system1D});

The problem is: there is any simple way to use command "subs" when some expression such that delta[1,1]=1, gamma[1,1]=1, gamma[1,2]^2=1 have value and others are zero.

Can someone please advice and help me on this?

thanks

witribm