#3;
pdsolve2(y(p_,t_)^m_+y(q_,x_)^n_+d_,y_):= if(isconstant(d), C_1+C_2*t^p/p!+(-d-C_2^m)^(1/n)*x^q/q!);
pdsolve2(y(p_,t_)^m_+y(p_,x_)^m_+c_,y_):= if(hasnot(c,t), t+dsolve(y(p,x)^m+c+1,y) );
#pdsolve2(y(p_,t_)^m_+b_*y(p_,x_)^m_+c_,y_):= if(hasnot(b,c,t), t^p/p!+int((-1/b-c/b)^(1/m),x)+gsolution(0,y,x,p) );
pdsolve2(y(p_,t_)^m_+b_*y(p_,x_)^m_+c_,y_):= if(hasnot(c,t), t^p/p!+dsolve(y(p,x)^m+c/b+1/b,y) );

pdsolve2(y(p_,t_)+y(q_,x_)^n_+d_,y_):= if(hasnot(d,t) and hasnot(d,y),if(p>1,dsolve(y(q,x)^n+d,y)+C_1*t, dsolve(y(q,x)^n+d-1,y)-t^p/p! ));

pdsolve2(y(p_,t_)+y(q_,x_)+d_,y_):= if(hasnot(d,y), pdsolve(y(p,t)+y(q,x),y)+psolution(1,q,0,d,y,t,x,p), 
	if(p==q,if(hasnot(d,t),replace(replace(dsolve(-d/2,y,x,q),C_1,C_1+t),c_1,c_1+t),  
	if(has(d,x),replace(dsolve(-d/2,y,x,q),C_1,C_1+t), replace(dsolve(-d/2,y,x,q),x,C_2*(t+x)))),
	if(q>1,dsolve(-d,y,t,p),
	if(p>1,dsolve(-d,y,x,q),
	dsolve(1-d,y,x,q)-t^p/p! ))));
pdsolve2(y(p_,t_)+b_*y(q_,x_)+d_,y_):= if(hasnot(b,d,y), pdsolve(y(p,t)+b*y(q,x),y)+psolution(b,q,0,d,y,t,x,p),
	if(p==q,if(hasnot(b,d,t),if(b+1==0, replace(dsolve(-d,y,x,p),x,2^(1/p)*t+x), replace(dsolve(-d,y,x,q),x,(t+x)/(1+b)^(1/p)) ),
	if(has(d,x),replace(dsolve(-d/(1+b),y,x,q),C_1,C_1+t),replace(dsolve(expand(-d/(1+b)),y,x,q),x,C_2*(t+x)))), 
	if(q>1,dsolve(-d,y,t,p),
	if(p>1,dsolve(-d/b,y,x,q),
	dsolve(1-d/b,y,x,q)-t^p/p! ))));
pdsolve2(a_*y(p_,t_)+y(q_,x_)+d_,y_):= if(hasnot(a,d,y), if(hasnot(a,x),C_1*(int(1/a,t,p)-x^q/q!),0)
	+if(hasnot(a,t),C_1*(t^p/p!-int(a,x,q)),0)+gsolution(0,y,t,p)+gsolution(0,y,x,q)
	+if(has(d/a,t),if(has(d/a,x), psolution(1/a,q,0,d/a,y,t,x,p), -int(d/a,t,p) ), -int(d,x,q) ), 
	if(hasnot(d,t), dsolve(-d,y,x,q)) );
pdsolve2(a_*y(p_,t_)+b_*y(q_,x_)+d_,y_):= if(hasnot(a,b,d,y), if(q<0,dsolve(d(a,x,-q)*y(p,t)+b*y+d(d,x,-q)),
	if(hasnot(a,x) and hasnot(b,t),C_1*(int(1/a,t,p)-int(1/b,x,q)),0)
	+if(hasnot(a,t) and hasnot(b,x),C_1*(int(b,t,p)-int(a,x,q)),0)
	+gsolution(0,y,t,p)+gsolution(0,y,x,q)
	+if(has(d/b,t),if(has(d/b,x), psolution(b/a,q,0,d/a,y,t,x,p), -int(d/a,t,p) ), -int(d/b,x,q) )), 
	if(hasnot(d/b,t), dsolve(-d/b,y,x,q)) );


#4;
pdsolve2(_f(y(1,x_))+b_*y_+y(p_,t_)+z_,y_):=dsolve(b*y+_f(y(1,x))+z,y)+gsolution(-b,y,t,p);


#5;
pdsolve2(a_+y_+z1_+z2_+z3_,y_):= if(isconstant(a) and hasnot(z1,y), pdsolve(y+z1+z2+z3,y)-a );
pdsolve2(a_+b_*y_+z1_+z2_+z3_,y_):= if(isconstant(a,b) and hasnot(z1,y), pdsolve(b*y+z1+z2+z3,y)-a/b );
pdsolve2(b_*y_+z1_+z2_+z3_+z4_,y_):= if(hasnot(z1,z2,x) and hasnot(z3,z4,t) and hasnot(z1,y), dsolve(b*y+z1+z2,y)+dsolve(b*y+z3+z4,y) );


#6;
pdsolve2(a_+y_+z1_+z2_+z3_+z4_,y_):= if(isconstant(a) and hasnot(z1,y), pdsolve(y+z1+z2+z3+z4,y)-a );
pdsolve2(a_+b_*y_+z1_+z2_+z3_+z4_,y_):= if(isconstant(a,b) and hasnot(z1,y), pdsolve(b*y+z1+z2+z3+z4,y)-a/b );


#0;
pdsolve2(a_):=pdsolve2(replace(toy(a),y(x),y),y);
pdsolve2(a_=b_):=pdsolve2(toy(a-b),y);