restart;intSet:={int1,int2}

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)

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

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)

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

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

Error, (in pdsolve) invalid input: `pdsolve/numeric` expects its 2nd argument, IBCs, to be of type {list, set}, but received (diff(sigma(r, z), r)+r*(diff(diff(sigma(r, z), r), r)))/r+(diff(theta(r, z), r)+r*(diff(diff(theta(r, z), r), r)))/r = 0