Question: Error, unable to compute coeff

hi

I want to equating coefficients of like powers of  in expr732vc and expr732wc, but

this error is appeared.please help me.

Error, unable to compute coeff
3.mw
 

restart;intSet:={int1,int2}

{int1, int2}

(1)

expr732 := 120000*(diff(v[0](x, T[0], T[1])+epsilon*v[1](x, T[0], T[1]), x, x))+1200*(diff(w[0](x, T[0], T[1])+epsilon*w[1](x, T[0], T[1]), x, x))*(diff(w[0](x, T[0], T[1])+epsilon*w[1](x, T[0], T[1]), x))+96.0000*(diff(diff(v[0](x, T[0], T[1])+epsilon*v[1](x, T[0], T[1]), T[0])+epsilon*(diff(v[0](x, T[0], T[1])+epsilon*v[1](x, T[0], T[1]), T[1])), x, x))+.9600*(diff(diff(w[0](x, T[0], T[1])+epsilon*w[1](x, T[0], T[1]), T[0])+epsilon*(diff(w[0](x, T[0], T[1])+epsilon*w[1](x, T[0], T[1]), T[1])), x))*(diff(w[0](x, T[0], T[1])+epsilon*w[1](x, T[0], T[1]), x, x))+.9600*(diff(diff(w[0](x, T[0], T[1])+epsilon*w[1](x, T[0], T[1]), T[0])+epsilon*(diff(w[0](x, T[0], T[1])+epsilon*w[1](x, T[0], T[1]), T[1])), x, x))*(diff(w[0](x, T[0], T[1])+epsilon*w[1](x, T[0], T[1]), x))-(diff(v[0](x, T[0], T[1])+epsilon*v[1](x, T[0], T[1]), T[0], T[0]))-epsilon*(diff(v[0](x, T[0], T[1])+epsilon*v[1](x, T[0], T[1]), T[1], T[0]))-epsilon*(diff(v[0](x, T[0], T[1])+epsilon*v[1](x, T[0], T[1]), T[1], T[0])+epsilon*(diff(v[0](x, T[0], T[1])+epsilon*v[1](x, T[0], T[1]), T[1], T[1]))), .9600*(diff(diff(v[0](x, T[0], T[1])+epsilon*v[1](x, T[0], T[1]), T[0])+epsilon*(diff(v[0](x, T[0], T[1])+epsilon*v[1](x, T[0], T[1]), T[1])), x, x))*(diff(w[0](x, T[0], T[1])+epsilon*w[1](x, T[0], T[1]), x))+.9600*(diff(diff(v[0](x, T[0], T[1])+epsilon*v[1](x, T[0], T[1]), T[0])+epsilon*(diff(v[0](x, T[0], T[1])+epsilon*v[1](x, T[0], T[1]), T[1])), x))*(diff(w[0](x, T[0], T[1])+epsilon*w[1](x, T[0], T[1]), x, x))+0.192e-1*(diff(diff(w[0](x, T[0], T[1])+epsilon*w[1](x, T[0], T[1]), T[0])+epsilon*(diff(w[0](x, T[0], T[1])+epsilon*w[1](x, T[0], T[1]), T[1])), x))*(diff(w[0](x, T[0], T[1])+epsilon*w[1](x, T[0], T[1]), x, x))*(diff(w[0](x, T[0], T[1])+epsilon*w[1](x, T[0], T[1]), x))+0.96e-2*(diff(w[0](x, T[0], T[1])+epsilon*w[1](x, T[0], T[1]), x))^2*(diff(diff(w[0](x, T[0], T[1])+epsilon*w[1](x, T[0], T[1]), T[0])+epsilon*(diff(w[0](x, T[0], T[1])+epsilon*w[1](x, T[0], T[1]), T[1])), x, x))+1200*(diff(v[0](x, T[0], T[1])+epsilon*v[1](x, T[0], T[1]), x, x))*(diff(w[0](x, T[0], T[1])+epsilon*w[1](x, T[0], T[1]), x))+1200*(diff(w[0](x, T[0], T[1])+epsilon*w[1](x, T[0], T[1]), x, x))*(diff(v[0](x, T[0], T[1])+epsilon*v[1](x, T[0], T[1]), x))+18*(diff(w[0](x, T[0], T[1])+epsilon*w[1](x, T[0], T[1]), x))^2*(diff(w[0](x, T[0], T[1])+epsilon*w[1](x, T[0], T[1]), x, x))-1.3943*(diff(w[0](x, T[0], T[1])+epsilon*w[1](x, T[0], T[1]), x, x, x, x))-0.111544e-2*(diff(diff(w[0](x, T[0], T[1])+epsilon*w[1](x, T[0], T[1]), T[0])+epsilon*(diff(w[0](x, T[0], T[1])+epsilon*w[1](x, T[0], T[1]), T[1])), x, x, x, x))+10-(diff(w[0](x, T[0], T[1])+epsilon*w[1](x, T[0], T[1]), T[0], T[0]))-epsilon*(diff(w[0](x, T[0], T[1])+epsilon*w[1](x, T[0], T[1]), T[1], T[0]))-epsilon*(diff(w[0](x, T[0], T[1])+epsilon*w[1](x, T[0], T[1]), T[1], T[0])+epsilon*(diff(w[0](x, T[0], T[1])+epsilon*w[1](x, T[0], T[1]), T[1], T[1]))):

expr732va := 1.20000*10^5*(diff(v[0](x, T[0], T[1]), x, x))+1200.*(diff(w[0](x, T[0], T[1]), x, x))*(diff(w[0](x, T[0], T[1]), x))+96.0000*(diff(v[0](x, T[0], T[1]), T[0], x, x))+.9600*(diff(w[0](x, T[0], T[1]), T[0], x))*(diff(w[0](x, T[0], T[1]), x, x))+.9600*(diff(w[0](x, T[0], T[1]), T[0], x, x))*(diff(w[0](x, T[0], T[1]), x))-1.*(diff(v[0](x, T[0], T[1]), T[0], T[0]))+1.20000*10^5*(diff(v[1](x, T[0], T[1]), x, x))+1200.*(diff(w[0](x, T[0], T[1]), x, x))*(diff(w[1](x, T[0], T[1]), x))+1200.*(diff(w[1](x, T[0], T[1]), x, x))*(diff(w[0](x, T[0], T[1]), x))+96.0000*(diff(v[1](x, T[0], T[1]), T[0], x, x))+96.0000*(diff(v[0](x, T[0], T[1]), T[1], x, x))+.9600*(diff(w[0](x, T[0], T[1]), T[0], x))*(diff(w[1](x, T[0], T[1]), x, x))+(.9600*(diff(w[1](x, T[0], T[1]), T[0], x)+diff(w[0](x, T[0], T[1]), T[1], x)))*(diff(w[0](x, T[0], T[1]), x, x))+.9600*(diff(w[0](x, T[0], T[1]), T[0], x, x))*(diff(w[1](x, T[0], T[1]), x))+(.9600*(diff(w[1](x, T[0], T[1]), T[0], x, x)+diff(w[0](x, T[0], T[1]), T[1], x, x)))*(diff(w[0](x, T[0], T[1]), x))-1.*(diff(v[1](x, T[0], T[1]), T[0], T[0]))-2.*(diff(v[0](x, T[0], T[1]), T[1], T[0])):

expr732vb1:=select(has,expr732[1],intSet)

0

(2)

expr732vb:=convert(series(map(s->remove(has,s,intSet)*subs(epsilon=0,select(has,s,intSet)),expr732vb1),epsilon,2),polynom);

0 = 0

(3)

expr732vc := 1.20000*10^5*(diff(v[0](x, T[0], T[1]), x, x))+1200.*(diff(w[0](x, T[0], T[1]), x, x))*(diff(w[0](x, T[0], T[1]), x))+96.0000*(diff(v[0](x, T[0], T[1]), T[0], x, x))+.9600*(diff(w[0](x, T[0], T[1]), T[0], x))*(diff(w[0](x, T[0], T[1]), x, x))+.9600*(diff(w[0](x, T[0], T[1]), T[0], x, x))*(diff(w[0](x, T[0], T[1]), x))-1.*(diff(v[0](x, T[0], T[1]), T[0], T[0]))+(1.20000*10^5*(diff(v[1](x, T[0], T[1]), x, x))+1200.*(diff(w[0](x, T[0], T[1]), x, x))*(diff(w[1](x, T[0], T[1]), x))+1200.*(diff(w[1](x, T[0], T[1]), x, x))*(diff(w[0](x, T[0], T[1]), x))+96.0000*(diff(v[1](x, T[0], T[1]), T[0], x, x))+96.0000*(diff(v[0](x, T[0], T[1]), T[1], x, x))+.9600*(diff(w[0](x, T[0], T[1]), T[0], x))*(diff(w[1](x, T[0], T[1]), x, x))+(.9600*(diff(w[1](x, T[0], T[1]), T[0], x)+diff(w[0](x, T[0], T[1]), T[1], x)))*(diff(w[0](x, T[0], T[1]), x, x))+.9600*(diff(w[0](x, T[0], T[1]), T[0], x, x))*(diff(w[1](x, T[0], T[1]), x))+(.9600*(diff(w[1](x, T[0], T[1]), T[0], x, x)+diff(w[0](x, T[0], T[1]), T[1], x, x)))*(diff(w[0](x, T[0], T[1]), x))-1.*(diff(v[1](x, T[0], T[1]), T[0], T[0]))-2.*(diff(v[0](x, T[0], T[1]), T[1], T[0])))*epsilon = 1.20000*10^5*(diff(v[0](x, T[0], T[1]), x, x))+1200.*(diff(w[0](x, T[0], T[1]), x, x))*(diff(w[0](x, T[0], T[1]), x))+96.0000*(diff(v[0](x, T[0], T[1]), T[0], x, x))+.9600*(diff(w[0](x, T[0], T[1]), T[0], x))*(diff(w[0](x, T[0], T[1]), x, x))+.9600*(diff(w[0](x, T[0], T[1]), T[0], x, x))*(diff(w[0](x, T[0], T[1]), x))-1.*(diff(v[0](x, T[0], T[1]), T[0], T[0]))+(1.20000*10^5*(diff(v[1](x, T[0], T[1]), x, x))+1200.*(diff(w[0](x, T[0], T[1]), x, x))*(diff(w[1](x, T[0], T[1]), x))+1200.*(diff(w[1](x, T[0], T[1]), x, x))*(diff(w[0](x, T[0], T[1]), x))+96.0000*(diff(v[1](x, T[0], T[1]), T[0], x, x))+96.0000*(diff(v[0](x, T[0], T[1]), T[1], x, x))+.9600*(diff(w[0](x, T[0], T[1]), T[0], x))*(diff(w[1](x, T[0], T[1]), x, x))+(.9600*(diff(w[1](x, T[0], T[1]), T[0], x)+diff(w[0](x, T[0], T[1]), T[1], x)))*(diff(w[0](x, T[0], T[1]), x, x))+.9600*(diff(w[0](x, T[0], T[1]), T[0], x, x))*(diff(w[1](x, T[0], T[1]), x))+(.9600*(diff(w[1](x, T[0], T[1]), T[0], x, x)+diff(w[0](x, T[0], T[1]), T[1], x, x)))*(diff(w[0](x, T[0], T[1]), x))-1.*(diff(v[1](x, T[0], T[1]), T[0], T[0]))-2.*(diff(v[0](x, T[0], T[1]), T[1], T[0])))*epsilon:

expr732wa:=convert(series(remove(has,expr732[2],intSet),epsilon,2),polynom):

expr732wb1:=select(has,expr732[2],intSet)

0

(4)

expr732wb:=convert(series(map(s->remove(has,s,intSet)*subs(epsilon=0,select(has,s,intSet)),expr732wb1),epsilon,2),polynom)

0 = 0

(5)

expr732wc := .9600*(diff(v[0](x, T[0], T[1]), T[0], x, x))*(diff(w[0](x, T[0], T[1]), x))+.9600*(diff(v[0](x, T[0], T[1]), T[0], x))*(diff(w[0](x, T[0], T[1]), x, x))+0.192e-1*(diff(w[0](x, T[0], T[1]), T[0], x))*(diff(w[0](x, T[0], T[1]), x, x))*(diff(w[0](x, T[0], T[1]), x))+0.96e-2*(diff(w[0](x, T[0], T[1]), x))^2*(diff(w[0](x, T[0], T[1]), T[0], x, x))+1200.*(diff(v[0](x, T[0], T[1]), x, x))*(diff(w[0](x, T[0], T[1]), x))+1200.*(diff(w[0](x, T[0], T[1]), x, x))*(diff(v[0](x, T[0], T[1]), x))+18.*(diff(w[0](x, T[0], T[1]), x))^2*(diff(w[0](x, T[0], T[1]), x, x))-1.3943*(diff(w[0](x, T[0], T[1]), x, x, x, x))-0.111544e-2*(diff(w[0](x, T[0], T[1]), T[0], x, x, x, x))+10.-1.*(diff(w[0](x, T[0], T[1]), T[0], T[0]))+(.9600*(diff(v[0](x, T[0], T[1]), T[0], x, x))*(diff(w[1](x, T[0], T[1]), x))+(.9600*(diff(v[1](x, T[0], T[1]), T[0], x, x)+diff(v[0](x, T[0], T[1]), T[1], x, x)))*(diff(w[0](x, T[0], T[1]), x))+.9600*(diff(v[0](x, T[0], T[1]), T[0], x))*(diff(w[1](x, T[0], T[1]), x, x))+(.9600*(diff(v[1](x, T[0], T[1]), T[0], x)+diff(v[0](x, T[0], T[1]), T[1], x)))*(diff(w[0](x, T[0], T[1]), x, x))+0.192e-1*(diff(w[0](x, T[0], T[1]), T[0], x))*(diff(w[0](x, T[0], T[1]), x, x))*(diff(w[1](x, T[0], T[1]), x))+(0.192e-1*((diff(w[0](x, T[0], T[1]), T[0], x))*(diff(w[1](x, T[0], T[1]), x, x))+(diff(w[1](x, T[0], T[1]), T[0], x)+diff(w[0](x, T[0], T[1]), T[1], x))*(diff(w[0](x, T[0], T[1]), x, x))))*(diff(w[0](x, T[0], T[1]), x))+0.96e-2*(diff(w[0](x, T[0], T[1]), x))^2*(diff(w[1](x, T[0], T[1]), T[0], x, x)+diff(w[0](x, T[0], T[1]), T[1], x, x))+0.192e-1*(diff(w[0](x, T[0], T[1]), x))*(diff(w[1](x, T[0], T[1]), x))*(diff(w[0](x, T[0], T[1]), T[0], x, x))+1200.*(diff(v[0](x, T[0], T[1]), x, x))*(diff(w[1](x, T[0], T[1]), x))+1200.*(diff(v[1](x, T[0], T[1]), x, x))*(diff(w[0](x, T[0], T[1]), x))+1200.*(diff(w[0](x, T[0], T[1]), x, x))*(diff(v[1](x, T[0], T[1]), x))+1200.*(diff(w[1](x, T[0], T[1]), x, x))*(diff(v[0](x, T[0], T[1]), x))+18.*(diff(w[0](x, T[0], T[1]), x))^2*(diff(w[1](x, T[0], T[1]), x, x))+36.*(diff(w[0](x, T[0], T[1]), x))*(diff(w[1](x, T[0], T[1]), x))*(diff(w[0](x, T[0], T[1]), x, x))-1.3943*(diff(w[1](x, T[0], T[1]), x, x, x, x))-0.111544e-2*(diff(w[1](x, T[0], T[1]), T[0], x, x, x, x))-0.111544e-2*(diff(w[0](x, T[0], T[1]), T[1], x, x, x, x))-1.*(diff(w[1](x, T[0], T[1]), T[0], T[0]))-2.*(diff(w[0](x, T[0], T[1]), T[1], T[0])))*epsilon = .9600*(diff(v[0](x, T[0], T[1]), T[0], x, x))*(diff(w[0](x, T[0], T[1]), x))+.9600*(diff(v[0](x, T[0], T[1]), T[0], x))*(diff(w[0](x, T[0], T[1]), x, x))+0.192e-1*(diff(w[0](x, T[0], T[1]), T[0], x))*(diff(w[0](x, T[0], T[1]), x, x))*(diff(w[0](x, T[0], T[1]), x))+0.96e-2*(diff(w[0](x, T[0], T[1]), x))^2*(diff(w[0](x, T[0], T[1]), T[0], x, x))+1200.*(diff(v[0](x, T[0], T[1]), x, x))*(diff(w[0](x, T[0], T[1]), x))+1200.*(diff(w[0](x, T[0], T[1]), x, x))*(diff(v[0](x, T[0], T[1]), x))+18.*(diff(w[0](x, T[0], T[1]), x))^2*(diff(w[0](x, T[0], T[1]), x, x))-1.3943*(diff(w[0](x, T[0], T[1]), x, x, x, x))-0.111544e-2*(diff(w[0](x, T[0], T[1]), T[0], x, x, x, x))+10.-1.*(diff(w[0](x, T[0], T[1]), T[0], T[0]))+(.9600*(diff(v[0](x, T[0], T[1]), T[0], x, x))*(diff(w[1](x, T[0], T[1]), x))+(.9600*(diff(v[1](x, T[0], T[1]), T[0], x, x)+diff(v[0](x, T[0], T[1]), T[1], x, x)))*(diff(w[0](x, T[0], T[1]), x))+.9600*(diff(v[0](x, T[0], T[1]), T[0], x))*(diff(w[1](x, T[0], T[1]), x, x))+(.9600*(diff(v[1](x, T[0], T[1]), T[0], x)+diff(v[0](x, T[0], T[1]), T[1], x)))*(diff(w[0](x, T[0], T[1]), x, x))+0.192e-1*(diff(w[0](x, T[0], T[1]), T[0], x))*(diff(w[0](x, T[0], T[1]), x, x))*(diff(w[1](x, T[0], T[1]), x))+(0.192e-1*((diff(w[0](x, T[0], T[1]), T[0], x))*(diff(w[1](x, T[0], T[1]), x, x))+(diff(w[1](x, T[0], T[1]), T[0], x)+diff(w[0](x, T[0], T[1]), T[1], x))*(diff(w[0](x, T[0], T[1]), x, x))))*(diff(w[0](x, T[0], T[1]), x))+0.96e-2*(diff(w[0](x, T[0], T[1]), x))^2*(diff(w[1](x, T[0], T[1]), T[0], x, x)+diff(w[0](x, T[0], T[1]), T[1], x, x))+0.192e-1*(diff(w[0](x, T[0], T[1]), x))*(diff(w[1](x, T[0], T[1]), x))*(diff(w[0](x, T[0], T[1]), T[0], x, x))+1200.*(diff(v[0](x, T[0], T[1]), x, x))*(diff(w[1](x, T[0], T[1]), x))+1200.*(diff(v[1](x, T[0], T[1]), x, x))*(diff(w[0](x, T[0], T[1]), x))+1200.*(diff(w[0](x, T[0], T[1]), x, x))*(diff(v[1](x, T[0], T[1]), x))+1200.*(diff(w[1](x, T[0], T[1]), x, x))*(diff(v[0](x, T[0], T[1]), x))+18.*(diff(w[0](x, T[0], T[1]), x))^2*(diff(w[1](x, T[0], T[1]), x, x))+36.*(diff(w[0](x, T[0], T[1]), x))*(diff(w[1](x, T[0], T[1]), x))*(diff(w[0](x, T[0], T[1]), x, x))-1.3943*(diff(w[1](x, T[0], T[1]), x, x, x, x))-0.111544e-2*(diff(w[1](x, T[0], T[1]), T[0], x, x, x, x))-0.111544e-2*(diff(w[0](x, T[0], T[1]), T[1], x, x, x, x))-1.*(diff(w[1](x, T[0], T[1]), T[0], T[0]))-2.*(diff(w[0](x, T[0], T[1]), T[1], T[0])))*epsilon:

eqEps:=seq(convert(map(s->coeff(s,epsilon,i)=0,[expr732vc,expr732wc]),exp),i=0..1)

Error, unable to compute coeff

 

``

``


 

Download 3.mw

 

Please Wait...