#2;
pdsolve1(y(p_,t_)+b_*y(q_,x_),y_):= if(b<0, C_2*exp((-b)^(1/p)*t+x),0)
	+if(hasnot(b,y), if(q>0, C_1*(if(hasnot(b,t),if(hasnot(b,x),t^p/p!-x^q/q!/b, t^p/p!-int(1/b,x,q)), int(b,t,p)-x^q/q!))
	+if(p>1 and q>1, C_3*t*x,0)+gsolution(0,y,t,p)+gsolution(0,y,x,q),0) );
pdsolve1(y(1,t_)+b_*y(1,x_),y_):= if(has(b,y), if(has(b,x),if(has(b,t),solve(b+d(b,t)/d(b,x),y)), solve(b=(C_1+x)/(C_2+t),y) ),
	if(hasnot(b,t),C_1+C_2*F(t-int(1/b,x)), if(hasnot(b,x),C_1+C_2*F(x-int(b,t)) )));
#pdsolve1(y(1,t_)+a_*y(2,x_),y_):= if(hasnot(a,t) and hasnot(a,x)),
	int(1/int(1/a,y),y)=C_1-t+x and int(1/int(1/a,y),y)=C_1-t-x );
pdsolve1(y(p_,t_)+a_*b_*y(q_,x_),y_):=  if(hasnot(a,b,y), if(hasnot(a,x) and hasnot(b,t), C_1*(int(a,t,p)-int(1/b,x,q)), if(hasnot(b,x) and hasnot(a,t), C_1*(int(b,t,p)-int(1/a,x,q)) )),
	if(p<=1 and q<=1,
	if(has(b,y), if(hasnot(a,x),solve(b=(C_1+x)/(C_2+int(a,t)),y), if(hasnot(a,t),solve(b=(C_1+int(1/a,x))/(C_2+t),y) )),
	if(has(a,y), if(hasnot(b,x),solve(a=(C_1+x)/(C_2+int(b,t)),y), if(hasnot(b,t),solve(a=(C_1+int(1/b,x))/(C_2+t),y) )) )), 
	a*b=(C_1+x)/(C_2+t) ));
#pdsolve1(y(1,t_)+a_*b_*y(2,x_),y_):= if(hasnot(b,t) and hasnot(a,b,x)),
	int(1/int(1/b,y),y)=C_1-int(a,t)+x and int(1/int(1/b,y),y)=C_1-int(a,t)-x );


#3; 
pdsolve1(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)
	+psolution(1/a,q,0,d/a,y,t,x,p) );
pdsolve1(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(d,y), if(a<0,C_2*((-a)^(1/p)*t+x),0)+psolution(a,q,0,d,y,t,x,p)+if(q>0,C_1,0)+C_3*(if(hasnot(b,x),int(b,t,p),t^p/p!)-if(hasnot(b,t),int(1/b,x,q),x^q/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(isconstant(a),psolution(b/a,q,0,d/a,y,t,x,p),psolution(a,p,b,q,0,d,y,x,t)) ))); 
pdsolve1(a_*y(p_,t_)+b_*y(p_,x_)+d_,y_):= if(hasnot(a,b,d,y), pdsolve(a*y(p,t)+b*y(p,x),y)
	+if(isconstant(a),psolution(b/a,p,0,d/a,y,t,x,p),psolution(a,p,b,p,0,d,y,t,x,p)) ); 
#pdsolve1(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,x), if(isconstant(d), if(p<=1, int((1+b)^(1/p),y)=C_1-d*t-d*x,if(p==2, int(int(-d,y)^(-1/2),y)=sqrt(2/(C_2^2+b*C_1^2))*(C_3+C_2*t+C_1*x) )),
	replace(dsolve(-d/b,y,t,p),C_1,C_1+C_2*(x-b*t)) ), 
	replace(dsolve(-d/b,y,x,q),C_1,C_1+C_2*(x-b*t)) ), 
	if(p<1 and hasnot(d,x) and hasnot(d,x),replace(pdsolve(y(1,t)+b*y(q,x)+d,y),t,t^p/p!),
	if(p>1,dsolve(-d/b,y,x,q),
	if(q>1,dsolve(-d,y,t,p),
	dsolve(1-d/b,y,x,q)-t^p/p! )))));
pdsolve1(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<1 and hasnot(d,t),replace(pdsolve(y(1,t)+b*y(q,x)+d,y),t,t^p/p!),
	if(p>1,dsolve(-d/b,y,x,q),
	if(q>1,dsolve(-d,y,t,p),
	dsolve(1-d/b,y,x,q)-t^p/p! ))));
pdsolve1(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<1 and hasnot(d,t),replace(pdsolve(y(1,t)+y(q,x)+d,y),t,t^p/p!),
	if(hasnot(d,x),dsolve(-d,y,t,p),
	if(hasnot(d,t),dsolve(-d,y,x,q),
	dsolve(1-d,y,x,q)-t^p/p! ))));
#pdsolve1(y(p_,t_)+b_*y(p_,x_)+d_,y_):= if(hasnot(b,d,y), pdsolve(y(p,t)+b*y(p,x),y)+psolution(b,p,0,d,y,t,x,p), 
	if(hasnot(d,t) and b<> -1, block(f1:=dsolve(-d/(1+b),y,x,p),if(has(f1,C_1+x),replace(f1,C_1+x,C_1+t+x),replace(f1,x,C_2+t+x))),  
	if(hasnot(d,x), block(f1:=dsolve(-d,y,t,p),if(hasnot(f1,C_1+t),replace(f1,C_1,C_1+C_2*(-b*t+x)),replace(f1,C_1,C_1+C_2*(-b*t^p+x^p))))   
	))); 
pdsolve1(y(p_,t_)+b_*y(p_,x_)+d_,y_):= if(hasnot(d,y), pdsolve(y(p,t)-y(p,x),y)+psolution(b,p,0,d,y,t,x,p), 
	if(hasnot(d,t), block(f1:=dsolve(-d/b,y,x,p),if(has(f1,C_1+x),replace(f1,C_1,C_1+C_2*(x-b*t)),replace(f1,C_1,C_1+C_2*(x^p-b*t^p)))),  
	if(hasnot(d,x), block(f1:=dsolve(-d/b,y,t,p),if(has(f1,C_1+t),replace(f1,C_1,C_1+C_2*(b*t-x)),replace(f1,C_1,C_1+C_2*(b*t^p-x^p))))   
	))); 
pdsolve1(y(p_,t_)+y(p_,x_)+d_,y_):= if(hasnot(d,y), pdsolve(y(p,t)+y(p,x),y)+psolution(1,p,0,d,y,t,x,p), 
	if(hasnot(d,t), block(f1:=dsolve(-d,y,x,p),if(has(f1,C_1+x),replace(f1,C_1,C_1+C_2*(x-t)),replace(f1,C_1,C_1+C_2*(x^p-t^p)))),  
	if(hasnot(d,x), block(f1:=dsolve(-d,y,t,p),if(has(f1,C_1+t),replace(f1,C_1,C_1+C_2*(t-x)),replace(f1,C_1,C_1+C_2*(t^p-x^p)))) 
	))); 


#pdsolve1(a_*b_+y(p_,t_)+c_*y(p_,x_),y_):= if(p<=1 and has(b,y) and hasnot(a,y), int(1/b,y,p)=pdsolve(y(p,t)+c*y(p,x),y)+psolution(c,p,0,a,y,t,x,p)*p!,
	if(p<=1 and has(a,y) and hasnot(b,y), int(1/a,y,p)=pdsolve(y(p,t)+c*y(p,x),y)+psolution(c,p,0,b,y,t,x,p)*p!,
	if(hasnot(b,y), pdsolve(y(p,t)+c*y(p,x),y)+psolution(c,p,0,a*b,y,t,x,p) )));
#pdsolve1(a_*b_+y(p_,t_)+y(p_,x_),y_):= if(p<=1 and has(b,y) and hasnot(a,y), int(1/b,y,p)=C_2+C_1*(t^p-x^p)+psolution(1,p,0,a,y,t,x,p)*p!,
	if(p<=1 and has(a,y) and hasnot(b,y), int(1/a,y,p)=C_2+C_1*(t^p-x^p)+psolution(1,p,0,b,y,t,x,p)*p!,
	if(hasnot(b,y), C_2+C_1*(t^p-x^p)+psolution(1,p,0,a*b,y,t,x,p) )));

#pdsolve1(a_*b_+y(p_,t_)+c_*y(p_,x_),y_):= if(p<=1 and has(b,y) and hasnot(a,y), solve(int(1/b,y,p)=pdsolve(y(p,t)+c*y(p,x),y)+psolution(c,p,0,a,y,t,x,p)*p!,y),
	if(p<=1 and has(a,y) and hasnot(b,y), solve(int(1/a,y,p)=pdsolve(y(p,t)+c*y(p,x),y)+psolution(c,p,0,b,y,t,x,p)*p!,y),
	if(hasnot(b,y), pdsolve(y(p,t)+c*y(p,x),y)+psolution(c,p,0,a*b,y,t,x,p) )));
#pdsolve1(a_*b_+y(p_,t_)+y(p_,x_),y_):= if(p<=1 and has(b,y) and hasnot(a,y), solve(int(1/b,y,p)=C_2+C_1*(t^p-x^p)+psolution(1,p,0,a,y,t,x,p)*p!,y),
	if(p<=1 and has(a,y) and hasnot(b,y), solve(int(1/a,y,p)=C_2+C_1*(t^p-x^p)+psolution(1,p,0,b,y,t,x,p)*p!,y),
	if(hasnot(b,y), C_2+C_1*(t^p-x^p)+psolution(1,p,0,a*b,y,t,x,p) )));

pdsolve(a_+y(1,t_)+b_*y_*y(1,x_),y_):= if(isconstant(a,b), y^2/2+y+C_1+a*t+a*x/b=0 );
pdsolve(a_+y(1,t_)+b_*log(y_)*y(1,x_),y_):= if(isconstant(a,b), log(y)*y=C_1-a*t-a*x/b );
pdsolve(a_+y(1,t_)+log(y_)*y(1,x_),y_):= if(hasnot(a,y), log(y)*y=C_1-a*t-a*x );
pdsolve(a_+y(1,t_)+y_*y(1,x_),y_):= if(hasnot(a,y), y^2/2+y+C_1+a*t+a*x=0 );

pdsolve1((a_+y_)^n_+z1_+z2_,y_):= if(d(a,x)==0 and d(a,t)==0, pdsolve(y^n+z1+z2,y)-a );
pdsolve1((a_+b_*y_)^n_+z1_+z2_,y_):= if(isconstant(a,b), pdsolve(y^n+z1+z2,y)/b^n-a/b );
pdsolve1(c_*(a_+y_)^n_+z1_+z2_,y_):= if(hasnot(c,y) and d(a,x)==0 and d(a,t)==0, pdsolve(c*y^n+z1+z2,y)-a );
pdsolve1(c_*(a_+y_)+z1_+z2_,y_):= if(hasnot(c,y) and d(a,x)==0 and d(a,t)==0, pdsolve(c*y+z1+z2,y)-a );
pdsolve1(y(p_,t_)+z1_+z2_,y_):= dsolve(z1+z2,y)+if(hasnot(z1,z2,y),gsolution(0,y,t,p),0);


pdsolve1(y(p1_,t_)+y(p2_,u_)+y(p3_,x_),y_):= C_1+C_2*(2*t^p1/p1!-u^p2/p2!-x^p3/p3!);
pdsolve1(y(p1_,t_)+y(p2_,u_)+c_*y(p3_,x_),y_):= if(hasnot(c,t) and hasnot(c,u), C_1+C_2*(2*t^p1/p1!-u^p2/p2!-int(1/c,x,p3)) );
pdsolve1(y(p1_,t_)+b_*y(p2_,u_)+y(p3_,x_),y_):= if(hasnot(b,t) and hasnot(b,x), C_1+C_2*(2*t^p1/p1!-int(1/b,u,p2)-x^p3/p3!) );
pdsolve1(y(p1_,t_)+b_*y(p2_,u_)+c_*y(p3_,x_),y_):= if(hasnot(b,t) and hasnot(b,x) and hasnot(c,t) and hasnot(c,u), C_1+C_2*(2*t^p1/p1!-int(1/b,u,p2)-int(1/c,x,p3)) );
#pdsolve1(y(p1_,t_)+b_*y(p2_,u_)+c_*y(p3_,x_),y_):= if(isconstant(b,c), C_1+C_2*(2*t^p1/p1!-u^p2/p2!/b-x^p3/p3!/c) );
#pdsolve1(y(p_,t_)+b_*y(2,x_)+c_*y(2,u_),y_):= if(isconstant(b,c), C_1+C_3*(2t^p/p!-x*x/2/b-u*u/2/c)+C_2*E(p,(b+c)*t^p)*sin(x)*sin(u) );
pdsolve1(y(1,t_)+u_*y(1,x_)-x_*y(1,u_),y_):= C_1+C_2*(x^2 + u^2)-t + atan(x/u);

pdsolve1(y(p_,t_)+y(m_,x_)+b_*y(q_,x_),y_):= if(hasnot(b,y), gsolution(-1/b,m,0,y,x,q)+gsolution(0,y,t,p)+C_3*(t^p/p!-x^m/m!) );
pdsolve1(y(p_,t_)+a_*y(m_,x_)+y(q_,x_),y_):= if(hasnot(a,y), gsolution(-a,m,0,y,x,q)+gsolution(0,y,t,p)+C_3*(t^p/p!-if(isconstant(a),x^m/m!/a,psolution(-a,m,0,1,y,x,q))), 
	dsolve(-a*y(m,x),y,x,q) );
pdsolve1(y(p_,t_)+a_*y(m_,x_)+b_*y(q_,x_),y_):= if(m==q, pdsolve(y(p,t)+(a+b)*y(m,x),y),
	if(hasnot(a,b,y), gsolution(-a/b,m,0,y,x,q)+gsolution(0,y,t,p)+C_3*(t^p/p!-if(isconstant(a),x^m/m!/a,psolution(-a/b,m,0,1/b,y,x,q))),
	replace(dsolve(-a*y(m,x),y,x,q),C_1,C_1+(-b+1)^(1/p)*t) ));


#4;
pdsolve1(x_*y(1,t_)+t_*y(1,x_)+b_*t_^2+b_*x_^2,y_):= C_1-b*t*x+C_2*(t^2-x^2);

pdsolve1(a_+y_+z1_+z2_,y_):= if(isconstant(a), pdsolve(y+z1+z2,y)-a );
pdsolve1(a_+b_*y_+z1_+z2_,y_):= if(isconstant(a,b), pdsolve(b*y+z1+z2,y)-a/b );

pdsolve1(y(p_,t_)+b_*y(p_,x_)+z1_+z2_,y_):= if(hasnot(b,y), replace(dsolve(z1+z2),x,x-b*t) );
pdsolve1(y(p_,t_)+y(p_,x_)+z1_+z2_,y_):= replace(dsolve(z1+z2),x,x-t);
pdsolve1(y(1,t_)+y_*y(1,x_)+y(p_,x_)+z2_,y_):=if(p>1,(C_1+x)/(C_2+t) );
pdsolve1(y(1,t_)+b_*y_*y(1,x_)+y(p_,x_)+z2_,y_):=if(p>1,(C_1+x)/(C_2+b*t) );

pdsolve1(y(p_,t_)+y(q_,x_)^n_+d_+f_,y_):= if(hasnot(d,x) and hasnot(f,t), 
	t^p/p!-int(d,t,p)+(-1)^(1+n)*dsolve(y(q,x)^n+f+1,y)+gsolution(0,y,t,p) );
pdsolve1(y(p_,t_)+b_*y(q_,x_)^n_+d_+f_,y_):= if(hasnot(d,x) and hasnot(f,t), 
	t^p/p!-int(d,t,p)+(-1)^(1+n)*dsolve(y(q,x)^n+f+1,y)*b^(-1/n)+gsolution(0,y,t,p) );

pdsolve1(y(p_,t_)+y(q_,x_)^n_+b_*y(m_,x_)+d_,y_):= if(hasnot(d,b,t),if(isconstant(d), dsolve(y(q,x)^n+b*y(m,x),y)-d*t^p/p!+gsolution(0,y,t,p),
	dsolve(y(q,x)^n+b*y(m,x)+d,y)+gsolution(-d,y,t,p) ));
pdsolve1(y(p_,t_)+a_*y(q_,x_)^n_+b_*y(m_,x_)+d_,y_):= if(hasnot(a,d,b,t),
	if(isconstant(d), dsolve(y(q,x)^n+b/a*y(m,x),y)-d*t^p/p!+gsolution(0,y,t,p),dsolve(y(q,x)^n+b/a*y(m,x)+d/a,y)+gsolution(-d,y,t,p) ));

pdsolve1(y(m_,t_)+b_*y(p_,x_)+c_*y(q_,x_)+d_,y_):= if(m<=1 and hasnot(b,c,d,t), replace(dsolve(-d/c,y,x,q),C_1+x,C_1+C_3*(x^p/p!-b*t^m/m!)) );
pdsolve1(y(m_,t_)+b_*y(p_,x_)+y(q_,x_)+d_,y_):= if(m<=1 and hasnot(b,d,t), replace(dsolve(-d,y,x,q),C_1+x,C_1+C_3*(x^p/p!-b*t^m/m!)) );
pdsolve1(y(m_,t_)+y(p_,x_)+c_*y(q_,x_)+d_,y_):= if(m<=1 and hasnot(c,d,t), replace(dsolve(-d/c,y,x,q),C_1+x,C_1+C_3*(x^p/p!-t^m/m!)) );
pdsolve1(y(m_,t_)+y(p_,x_)+y(q_,x_)+d_,y_):= if(m<=1 and hasnot(d,t), replace(dsolve(-d,y,x,q),C_1+x,C_1+C_3*(x^p/p!-t^m/m!)) );

pdsolve1(y(q_,t_)+b_*y(p_,x_)+c_*y(q_,x_)+d_,y_):= if(hasnot(b,c,d,t), replace(dsolve(-d/b,y,x,p),x,C_2*(x+(-c)^(1/q)*t)) );
pdsolve1(y(q_,t_)+b_*y(p_,x_)+y(q_,x_)+d_,y_):= if(hasnot(b,d,t), replace(dsolve(-d/b,y,x,p),x,C_2*(x+(-1)^(1/q)*t)) );
pdsolve1(y(q_,t_)+y(p_,x_)+c_*y(q_,x_)+d_,y_):= if(hasnot(c,d,t), replace(dsolve(-d,y,x,p),x,C_2*(x+(-c)^(1/q)*t)) );
pdsolve1(y(q_,t_)+y(p_,x_)+y(q_,x_)+d_,y_):= if(hasnot(d,t), replace(dsolve(-d,y,x,p),x,C_2*(x+(-1)^(1/q)*t)) );

pdsolve1(y(p_,t_)+b_*y(p_,x_)+c_*y(q_,x_)+d_,y_):= if(hasnot(b,c,d,t), replace(dsolve(-d/c,y,x,q),C_1+x,C_1+C_3*(x^p-b*t^p)) );
pdsolve1(y(p_,t_)+b_*y(p_,x_)+y(q_,x_)+d_,y_):= if(hasnot(b,d,t), replace(dsolve(-d,y,x,q),C_1+x,C_1+C_3*(x^p-b*t^p)) );
pdsolve1(y(p_,t_)+y(p_,x_)+c_*y(q_,x_)+d_,y_):= if(hasnot(c,d,t), replace(dsolve(-d/c,y,x,q),C_1+x,C_1+C_3*(x^p-t^p)) );
pdsolve1(y(p_,t_)+y(p_,x_)+y(q_,x_)+d_,y_):= if(hasnot(d,t), replace(dsolve(-d,y,x,q),C_1+x,C_1+C_3*(x^p-t^p)) );

#0;
pdsolve1(a_):=pdsolve1(toy(a),y);
pdsolve1(a_=b_):=pdsolve1(toy(a-b),y);
pdsolve1(ds(y_,t_)=a_+b_*y_):=if(hasnot(a,b,t), pdsolve(a+b*y)+C_3*exp(b*t),pdsolve(y(1,t)-a-b*y));
pdsolve1(ds(y_,t_,p_)=a_+b_*y_):=if(hasnot(a,b,t), pdsolve(a+b*y)+C_3*E(p,b*t^p),pdsolve(y(p,t)-a-b*y));