pdsolve(a_*b_,y_):=pdsolve(a,y) and pdsolve(b,y);


#2;
pdsolve(y(p_,t_)*y(q_,x_)+d_*y_^2,y_):= if(isconstant(d), C_1*exp((-d)^(1/2/p)*t+(-d)^(1/2/q)*x) );
pdsolve(y(p_,t_)*y(q_,x_)+d_*(y_+c_)^2,y_):= if(isconstant(c,d), -c+C_1*exp((-d)^(1/2/p)*t+(-d)^(1/2/q)*x) );
pdsolve(y(p_,t_)*y(q_,x_)+d_*(b_*y_+c_)^2,y_):= if(isconstant(b,c,d), -c/b+C_1*exp(b^(1/p)*(-d)^(1/2/p)*t+b^(1/q)*(-d)^(1/2/q)*x) );

pdsolve(y(1,t_)^n_+y_^m_*y(1,x_)^n_,y_):=(-1)^(1/m+n)*(C_1+m*x)^(n/m)/(C_2+m*t)^(n/m);
pdsolve(y(1,t_)^n_+y_*y(1,x_)^n_,y_):=(-1)^(n+1)*(C_1+x)^n/(C_2+t)^n;

pdsolve(y(p_,t_)+b_*y(q_,x_)^n_,y_):= if(isconstant(b),-C_2^n*b*t^p/p!+C_2*x^q/q!+gsolution(0,y,t,p)+gsolution(0,y,x,q) );
pdsolve(y(p_,t_)+y(q_,x_)^n_,y_):= -C_2^n*t^p/p!+C_2*x^q/q!+gsolution(0,y,t,p)+gsolution(0,y,x,q);
pdsolve(y(p_,t_)*y_+y(q_,x_)^2,y_):=C_1*exp((-1)^(1/p)*t+(-1)^(1/q)*x);
pdsolve(y(p_,t_)^2+y_*y(q_,x_),y_):=C_1*exp((-1)^(1/p)*t+(-1)^(1/q)*x);
pdsolve(y(p_,t_)*y_-y(q_,x_)^2,y_):=C_1*exp(t+x);
pdsolve(y(p_,t_)^2-y_*y(q_,x_),y_):=C_1*exp(t+x);
pdsolve(y(2,t_)*y_-y(1,x_)^2,y_):=C_1*exp(C_2*t+C_2*x);
pdsolve(-y(2,x_)*y_+y(1,t_)^2,y_):=C_1*exp(C_2*t+C_2*x);

#pdsolve(y(p_,t_)+b_*y(q_,x_),y_):= if(b<0, C_2*exp((-b)^(1/p)*t+x)+if(q>0,C_1+C_3*(t^p/p!-int(1/b,x,q)),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_2*t*x,0)+gsolution(0,y,t,p)+gsolution(0,y,x,q),0) )) ;
pdsolve(a_*y(p_,t_)+y(q_,x_),y_):= 
	if(hasnot(a,t), C_1*(t^p/p!-int(a,x,q)),0)
	+if(hasnot(a,x), C_1*(int(1/a,t,p)-x^q/q!),0)
	+if(p>1 and q>1, C_2*t*x,0)+gsolution(0,y,t,p)+gsolution(0,y,x,q) ;
pdsolve(a_*y(p_,t_)+b_*y(q_,x_),y_):= 
	if(hasnot(a,t) and hasnot(b,x), C_1*(int(b,t,p)-int(a,x,q)),
	if(hasnot(a,x) and hasnot(b,t), C_1*(int(1/a,t,p)-int(1/b,x,q)) ))
	+if(p>1 and q>1, C_2*t*x,0)+gsolution(0,y,t,p)+gsolution(0,y,x,q) ;
pdsolve(y(p_,t_)+y(q_,x_),y_):= if(p>0 and q>0, C_2*(t^p/p!-x^q/q!)+gsolution(0,y,t,p)+gsolution(0,y,x,q),0)
	+if(p>1 and q>1, C_2*t*x,0)
	+if(p<0 or q<0,C_1*exp(t-x),0) ;

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

#pdsolve(y(p_,t_)+a_*b_*y(q_,x_),y_):= if(p<=1 and q<=1,
	if(hasnot(a,y) and 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) and hasnot(b,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) ))
	)));

pdsolve(y(p_,t_)+a_*b_*c_*y(q_,x_),y_):= if(p<=1,
	if(has(c,y), if(hasnot(a,x) and hasnot(b,t), c=(int(1/b,x,q)+C_2)/(C_1+int(a,t,p)), if(hasnot(a,t) and hasnot(b,x), c=((int(1/a,x,q)+C_2)/(C_1+int(b,t,p))) )),
	if(has(b,y), if(hasnot(a,x) and hasnot(c,t),solve(b=(C_1+int(1/c,x,q))/(C_2+int(a,t)),y), if(hasnot(a,t) and hasnot(c,x),solve(b=(C_1+int(1/a,x))/(C_2+int(c,t,p)),y) )),
	if(has(a,y), if(hasnot(b,x) and hasnot(c,t),solve(a=(C_1+int(1/c,x,q))/(C_2+int(b,t)),y), if(hasnot(b,t) and hasnot(c,x),solve(a=(C_1+int(1/b,x))/(C_2+int(c,t,p)),y) 
	))))));


pdsolve(y(p_,t_)+b_*y_*y(q_,x_),y_):=  if(hasnot(b,y) and p<=1,
	if(hasnot(b,x), (x^q/q!+gsolution(0,y,x,q))/(C_1+int(b,t,p)),
	if(hasnot(b,t), (int(1/b,x,q)+gsolution(0,y,x,q))/(C_1+t^p/p!) 
	)));
pdsolve(y(p_,t_)+a_*b_*y_*y(q_,x_),y_):=  if(hasnot(a,b,y) and p<=1,
	if(hasnot(a,x) and hasnot(b,t), (int(1/b,x,q)+gsolution(0,y,x,q))/(C_1+int(a,t,p)), 
	if(hasnot(a,t) and hasnot(b,x), (int(1/a,x,q)+gsolution(0,y,x,q))/(C_1+int(b,t,p))
	)));

pdsolve(y(p_,t_)+c_*exp(a_+b_)*y(q_,x_),y_):= if(hasnot(a,c,x) and hasnot(b,y),C_1*(int(c*exp(a),t,p)-int(exp(-b),x,q)) );
pdsolve(y(p_,t_)+b_*exp(c_+y_)*y(q_,x_),y_):=  if(isconstant(b),if(d(c,t,p)==0 and d(c,x,q)==0, -c-p*log(C_1+b*t)+log((-1)^(p-q)*(x)^q*p!/q!),
	if(hasnot(c,t), log(C_1+(-1)^(p-q)*p!/q!*int(exp(-c),x,q))-log((if(p>1,0,C_2)+b*t)), 
	if(hasnot(c,x), log(if(q>1,0,C_2)+x^q/q!)-log(if(p>1,0,C_1)+b*int(exp(c),t,p))  ))));
pdsolve(y(p_,t_)+exp(a_+b_)*f1_*y(q_,x_),y_):= if(hasnot(a,f1,x) and hasnot(b,y),C_1*(int(f1*exp(a),t,p)-int(exp(-b),x,q)) );
pdsolve(y(p_,t_)+exp(y_)*f1_*y(q_,x_),y_):=  if(hasnot(f1,x), -log(if(p>1,0,C_1)+int(f1,t,p))+log((-1)^(p-q)*p!/q!*(if(q>1,0,C_2)+x)^q) );
pdsolve(y(p_,t_)+b_*exp(y_)*y(q_,x_),y_):=  if(hasnot(b,x), -log(if(p>1,0,C_1)+p!*int(b,t,p))+log((-1)^(p-q)*p!/q!*(if(q>1,0,C_2)+x)^q) );
#pdsolve(y(p_,t_)+b_*exp(y_)*y(q_,x_),y_):=  if(isconstant(b), -log(b*(if(p>1,0,C_1)+t)^p)+log((-1)^(p-q)*p!/q!*(if(q>1,0,C_2)+x)^q),
	if(hasnot(b,x), -log(if(p>1,0,C_1)+int(b,t,p))+log((-1)^(p-q)*p!/q!*(if(q>1,0,C_2)+x)^q) 
	));
pdsolve(y(p_,t_)+exp(a_+b_)*y_*y(q_,x_),y_):= if(hasnot(a,x) and hasnot(b,t), (-1)^p/int(exp(a),t,p)*(gsolution(0,y,x,q)-int(exp(-b),x,q)) );
#pdsolve(y(p_,t_)+b_*y_^n_*y(q_,x_),y_):= if(hasnot(b,x), (1/(C_1+int(b,t,p))*(gsolution(0,y,x,q)+x^q/q!))^(1/n), 
	if(hasnot(b,t), (1/(C_1+factorial(p)*t^p)*(gsolution(0,y,x,q)+int(1/b,x,q)))^(1/n)
	));
pdsolve(y(p_,t_)+b_*y_^n_*y(q_,x_),y_):=  if(hasnot(b,y),
	if(hasnot(b,t), (-1)^(1/n)*(C_1+fallingfactorial(-p/n,p)^(-1/p)*t)^(-p/n)*(gsolution(0,y,x,q)+fallingfactorial(q/n,q)^(-1/q)*int(1/b,x,q))^(q/n), 
	if(hasnot(b,x), (1/(C_1+int(b,t,p))*(gsolution(0,y,x,q)+x^q/q!))^(1/n)
	)));

pdsolve(y(p_,t_)+c_*exp(a_+b_)*y_*y(q_,x_),y_):= if(hasnot(a,c,x) and hasnot(b,c,t), 
	(-1)^p*(gsolution(0,y,x,q)-int(exp(-b),x,q))/int(exp(a),t,p)/c );
pdsolve(y(p_,t_)+b_*exp(y_)*f1_*y(q_,x_),y_):= if(hasnot(b,f1,x), -log(if(p>1,0,C_1)+b*int(f1,t,p))+log((-1)^(p-q)*p!/q!*(if(q>1,0,C_2)+x)^q) );
pdsolve(y(p_,t_)+a_*b_*y_^n_*y(q_,x_),y_):= if(p<=1 and hasnot(a,b,y),
	if(hasnot(a,x) and hasnot(b,t), ((int(1/b,x,q)+gsolution(0,y,x,q))/(C_1+int(a,t,p)))^(1/n), 
	if(hasnot(a,t) and hasnot(b,x), ((int(1/a,x,q)+gsolution(0,y,x,q))/(C_1+int(b,t,p)))^(1/n)
	)));

pdsolve(y(p_,t_)+exp(a_+b_)*y(q_,x_),y_):= if(hasnot(a,x) and hasnot(b,y),C_1*(int(exp(a),t,p)-int(exp(-b),x,q)) );
#pdsolve(y(p_,t_)+exp(c_+y_)*y(q_,x_),y_):= if(hasnot(c,t), log((-1)^(p-q)*p!/q!*int(exp(-c),x,q))-p*log((C_1+t)), 
	if(hasnot(c,x), log(x^q/q!)-log(C_2+int(exp(c),t,p))  ));
pdsolve(y(p_,t_)+exp(c_+y_)*y(q_,x_),y_):=  if(d(c,t,p)==0 and d(c,x,q)==0, -c-p*log(C_1+t)+log((-1)^(p-q)*(x)^q*p!/q!),
	if(hasnot(c,t), log(C_1+(-1)^(p-q)*p!/q!*int(exp(-c),x,q))-log((if(p>1,0,C_2)+t)), 
	if(hasnot(c,x), log(if(q>1,0,C_2)+x^q/q!)-log(if(p>1,0,C_1)+int(exp(c),t,p))  )));
pdsolve(y(p_,t_)+exp(y_)*y(q_,x_),y_):= -p*log(C_1+t)+log((-1)^(p-q)*(if(q>1,0,C_2)+x)^q*p!/q!);

pdsolve(y(p_,t_)+y_^n_*y(q_,x_),y_):= -(C_1+t)^(-p/n)*(C_2+x)^(q/n)*fallingfactorial(-p/n,p)^sgn(n)/fallingfactorial(q/n,q)^sgn(n);
pdsolve(y(p_,t_)+y_*y(q_,x_),y_):= -(C_1+t)^(-p)*(x^q+gsolution(0,y,x,q))*fallingfactorial(-p,p)/q!;
pdsolve(y(p_,t_)+(b_+y_)*y(q_,x_),y_):= if(d(b,t,p)==0 and d(b,x,q)==0, -b-fallingfactorial(-p,p)*(C_1+t)^(-p)*(x^q+gsolution(0,y,x,q))/q!);
pdsolve(y(p_,t_)+(b_+c_*y_)*y(q_,x_),y_):= if(isconstant(c) and d(b,t,p)==0 and d(b,x,q)==0, -b/c-1/c*fallingfactorial(-p,p)*(C_1+t)^(-p)*(x^q+gsolution(0,y,x,q))/q! );

	
#pdsolve(y(p_,t_)+a_*y(2,x_),y_):= if(hasnot(a,y), if(hasnot(a,x) and p<=1, C_1+C_2*(t^p/p!+x-x^2/2/a)+C_3*E(p,a*t^p)*(sin(x)+cos(x)), C_1+C_2*x+int(-1/a,x,2)+t^p/p! ));
#pdsolve(y(2,t_)+b_*y(2,x_),y_):= if(isconstant(b), C_1+C_3*(t+x+t^2-x^2/b)+C_2*(sin(x)+cos(x))*if(b>0,(exp(-sqrt(b)*t)+exp(sqrt(b)*t)),(sin(sqrt(-b)*t)+cos(sqrt(-b)*t)) ) );
#pdsolve(y(1,t_)+a_*y(1,x_),y_):= if(hasnot(a,y), if(hasnot(a,x), C_1+C_2*(t-x/a), C_1+t-int(1/a,x) ));
#pdsolve(y(1,t_)+b_*y(1,x_),y_):= if(hasnot(b,y), if(hasnot(b,x), C_1+C_2*(t-x/b), C_1+t-int(1/b,x) ), if(hasnot(b,t) and hasnot(b,x),solve(b=(C_1+x)/(C_2+t),y) ));
#pdsolve(a_*y(1,t_)+y(1,x_),y_):= if(has(a,y), if(hasnot(a,x),solve(a=(C_1+t)/(C_2+x),y),
	if(hasnot(a,t),solve(a=(C_1+t)/x,y) )),C_1+C_2*(x-int(1/a,t)) );
#pdsolve(a_*y(1,t_)+b_*y(1,x_),y_):= if(hasnot(a,t) and hasnot(b,x), solve(int(a,x)+int(b,t)+d(b,t)/d(b,x),y),
	if(hasnot(a,x) and hasnot(b,t), C_1*(int(1/a,t)-int(1/b,x)) ));

pdsolve(y(2,t_)+a_*y(1,x_),y_):= if(has(a,y), int(1/int(a,y),y)=C_1+t-x and int(1/int(a,y),y)=C_1-t-x );
pdsolve(y(1,t_)+a_*y(2,x_),y_):= if(has(a,y), int(1/int(1/a,y),y)=C_1-t+x and int(1/int(1/a,y),y)=C_1-t-x );
pdsolve(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)) )));

pdsolve(y(2,t_)+a_*b_*y(1,x_),y_):= if(hasnot(a,y),
	if(has(b,y), if(hasnot(a,t), int(1/int(b,y),y)=C_1-int(1/a,x)+t and int(1/int(b,y),y)=C_1-int(1/a,x)-t),
	if(has(a,y), if(hasnot(b,t), int(1/int(a,y),y)=C_1-int(1/b,x)+t and int(1/int(a,y),y)=C_1-int(1/b,x)-t) 
	)));
pdsolve(y(1,t_)+a_*b_*y(2,x_),y_):= if(hasnot(a,y),
	if(has(b,y), if(hasnot(a,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),
	if(has(a,y), if(hasnot(b,x), int(1/int(1/a,y),y)=C_1-int(b,t)+x and int(1/int(1/a,y),y)=C_1-int(b,t)-x) 
	)));
pdsolve(y(2,t_)+y_*y(1,x_),y_):= 2/(C_1+t+x) and -2/(C_1+t-x);
pdsolve(y(p_,t_)+b_*y(m_,x_)*y(q_,x_),y_):=gsolution(0,y,t,p)+gsolution(0,y,x,q)+if(p>1 and q>1,C_3*t*x,0);

pdsolve(y(1,t_)+(x_+b_)*y(1,x_),y_):= C_1*F(exp(-t)*x-int(exp(-t)*b,t))+C_2;
pdsolve(y(1,t_)+(x_+t_)*y(1,x_),y_):= C_1*F(exp(-t)*(t+x+1))+C_2;
pdsolve(y(1,t_)+(a_*x_+b_)*y(1,x_),y_):= C_1*F(exp(-int(a,t))*x-int(exp(-int(a,t))*b,t))+C_2;
pdsolve(y(1,t_)+(a_*t_+b_*x_)*y(1,x_),y_):= C_1*F(replace(int(1/(a+b*v),v)-t,v,a*t+b*x))+C_2;
pdsolve(y(1,t_)+(a_*t_+b_*x_+c_)*y(1,x_),y_):= C_1*F(replace(int(1/(a+b*v),v)-t,v,a*t+b*x+c))+C_2;
pdsolve(x_*y(1,t_)+t_*y(1,x_),y_):= C_1+C_2*F(t^2-x^2);
pdsolve(t_*y(1,t_)+x_*y(1,x_),y_):= C_1+C_2*F(t/x);


#3;
pdsolve(y(p_,t_)+y(q_,x_)^n_+c_*y_,y_):= if(isconstant(c), dsolve((-c*y)^(1/n),y,x,q)+gsolution(-c,y,t,p) );
pdsolve(y(p_,t_)+y(q_,x_)^n_+y_,y_):= dsolve((-y)^(1/n),y,x,q)+gsolution(-1,y,t,p);

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

pdsolve(y(p_,t_)+a_*y(1,x_)^2+b_*y(2,x_),y_):= if(isconstant(a), if(hasnot(b,y),dsolve(b*y(2,x)+a*y(1,x)^2-1,y)-t^p/p!+gsolution(0,y,t,p),
	dsolve(b*y(2,x)+a*y(1,x)^2,y)+gsolution(0,y,t,p) ) and C_2*x-a*C_2^2*t^p/p!+gsolution(0,y,t,p), 
	if(d(b,y)==a, if(p==1, solve(int(b/y,y)=C_1-t+x,y) if(p==2, solve(b=-(C_1+x)^2/(C_2+t)^2,y) ))));
pdsolve(y(p_,t_)+exp(y_)*y(1,x_)^2+exp(y_)*y(2,x_),y_):= -p*log(C_1+t)+log((-1)^(p+1)*p!*x^2/2);
pdsolve(y(p_,t_)+b_*exp(y_)*y(1,x_)^2+b_*exp(y_)*y(2,x_),y_):= -p*log(C_1+t)+log((-1)^(p+1)*p!*x^2/2/b);

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

pdsolve(y(p_,t_)+a_*y(1,x_)^m_+y_*y(q_,x_),y_):= if(isconstant(a),C_2*x-C_2^m*a*t^p/p!+gsolution(0,y,t,p) );
pdsolve(y(p_,t_)+a_*y(1,x_)^m_+b_*y_*y(q_,x_),y_):= if(isconstant(a),C_2*x-C_2^m*a*t^p/p!+gsolution(0,y,t,p) );
pdsolve(y(p_,t_)+y(1,x_)^2+y_*y(1,x_),y_):= C_2-exp(C_1-x+t*C_2)/C_2;
pdsolve(y(p_,t_)+y(1,x_)^2+a_*y_*y(1,x_),y_):= if(isconstant(a),C_2/a-exp((C_1-a*x+a^(1/p)*t*C_2))/C_2/a );
pdsolve(y(p_,t_)+b_*y(1,x_)^2+y_*y(1,x_),y_):= if(isconstant(b),-C_2-exp(C_1-1/b*x)*E(p,-1/b*t^p*C_2)/C_2 );
pdsolve(y(p_,t_)+b_*y(1,x_)^2+a_*y_*y(1,x_),y_):= if(isconstant(a,b),-C_2/a-exp(C_1-a/b*x)*E(p,-a/b*t^p*C_2)/C_2/a );

#pdsolve(y_+y(p_,t_)+y(q_,x_),y_):=if(p<=1 and q<=1,C_2+t^p/p!-x^q/q!,1)*gsolution(-1,y,t,p)*if(p<=1 and q>1,gsolution(0,y,x,q),1)
	+if(q<=1 and p<=1,C_2+t^p/p!-x^q/q!,1)*gsolution(-1,y,x,q)*if(q<=1 and p>1,gsolution(0,y,t,p),1);
#pdsolve(y_+y(p_,t_)+c_*y(q_,x_),y_):=if(hasnot(c,y),if(p<0 or q<0 and isconstant(c), C_1*exp((-1/2)^(1/p)*t+(-1/2/c)^(1/q)*x),
	if(p<=1 and q<=1,if(hasnot(c,x),C_2+int(c,t,p)-x^q/q!,C_2+t^p/p!-int(1/c,x,q)),1)*gsolution(-1,y,t,p)
	+if(p<=1 and q<=1,C_2+t^p/p!-int(1/c,x,q),1)*if(hasnot(c,t),gsolution(-1/c,y,x,q),0) ));
#pdsolve(a_*y_+y(p_,t_)+c_*y(q_,x_),y_):=if(hasnot(a,x),gsolution(-a,y,t,p),0)*if(p<=1 and q<=1,C_2+if(hasnot(c,x),int(c,t,p),t^p/p!)-if(has(c,x),int(1/c,x,q),x^q/q!/q),1)
	+if(hasnot(a,t),gsolution(-a/c,y,x,q),0)*if(p<=1 and q<=1,C_2+if(hasnot(c,x),int(c,t,p),t^p/p!)-if(has(c,x),int(1/c,x,q),x^q/q!/q),1); 
#pdsolve(a_*y_+y(p_,t_)+c_*y(q_,x_),y_):=if(p<=1 and q<=1 and isconstant(a),if(hasnot(c,x),C_2+int(c,t,p)-x^q/q!,C_2+t^p/p!-int(1/c,x,q)),1)*if(hasnot(a,x),gsolution(-a,y,t,p),0)*if(p<=1 and q>1,gsolution(0,y,x,q),1)
	+if(p==1 and q==1 and isconstant(c),C_1*exp(-int(replace(replace(a,t,z),x,x-c*t+c*z),z,0,t))*F(t-x/c),0)
	+if(p<=1 and q<=1 and isconstant(c),if(hasnot(c,x),C_2+int(c,t,p)-x^q/q!,C_2+t^p/p!-int(1/c,x,q),1)*if(hasnot(a/c,t),gsolution(-a/c,y,x,q),0),0);
#pdsolve(a_*y_+y(p_,t_)+c_*y(q_,x_),y_):=if(p<0 or q<0 and isconstant(a,c), C_1*exp((-a/2)^(1/p)*t+(-a/2/c)^(1/q)*x),
	if(p<=1 and q<=1 and isconstant(a),if(hasnot(c,x),C_2+int(c,t,p)-x^q/q!,C_2+t^p/p!-int(1/c,x,q)),1)*if(hasnot(a,x),gsolution(-a,y,t,p),0)*if(p<=1 and q>1,gsolution(0,y,x,q),1)
	+if(p==1 and q==1 and isconstant(c),C_1*exp(-int(replace(replace(a,t,z),x,x-c*t+c*z),z,0,t))*(t-x/c),0)
	+if(p<=1 and q<=1 and isconstant(c),if(hasnot(c,x),C_2+int(c,t,p)-x^q/q!,C_2+t^p/p!-int(1/c,x,q))*if(hasnot(a/c,t),gsolution(-a/c,y,x,q),0),0) );
#pdsolve(a_*y_+y(p_,t_)+y(q_,x_),y_):=if(p<=1 and q<=1,C_2+t^p/p!-x^q/q!,1)*if(hasnot(a,x),gsolution(-a,y,t,p),0)
	+if(p<=1,C_2+t^p/p!-x^q/q!,1)*if(hasnot(a,x),gsolution(-a,y,t,p),0)
	+if(q<=1,C_2+t^p/p!-x^q/q!,1)*if(hasnot(a,t),gsolution(-a,y,x,q),0)
	+if(p<=1 and q<=1,(t^p/p!-x^q/q!),1)*if(has(a,t) and has(a,x),gsolution(-replace(a,x,0),y,t,p)*gsolution(-replace(a,t,0),y,x,q),0);
#pdsolve(a_*y_+y(p_,t_)+y(q_,x_),y_):= C_1*exp(psolution(1,q,0,a,y,t,x,p))*(x^q/q!-t^p/p!);

pdsolve(y_+b_*y(p_,t_)+c_*y(q_,x_),y_):= if(p<0 or q<0 and isconstant(b,c), C_1*exp((-1/2/b)^(1/p)*t+(-1/2/c)^(1/q)*x),
	if(hasnot(b,x),gsolution(-1/b,y,t,p),0)
	+if(hasnot(c,t),gsolution(-1/c,y,x,q),0) ); 
pdsolve(a_*y_+b_*y(p_,t_)+c_*y(q_,x_),y_):=	if(hasnot(a,y),if(p<0 or q<0 and isconstant(a,b,c), C_1*exp((-a/2/b)^(1/p)*t+(-a/2/c)^(1/q)*x),
	if(hasnot(a,b,x),gsolution(-a/b,y,t,p),0)*if(has(a,t),t/x,1)
	+if(hasnot(a,c,t),gsolution(-a/c,y,x,q),0)*if(has(a,x),t/x,1) 
	)); 
pdsolve(y_+y(p_,t_)+c_*y(q_,x_),y_):= if(hasnot(c,y),
	if(p<0 or q<0 and isconstant(c), C_1*exp((-1/2)^(1/p)*t+(-1/2/c)^(1/q)*x),
	if(hasnot(c,t),gsolution(-1,y,t,p)+gsolution(-1/c,y,x,q),
	if(p<=1 and q<=1,if(hasnot(c,x),C_2+int(c,t,p)-x^q/q!,C_2+t^p/p!-int(1/c,x,q)),1)*gsolution(-1,y,t,p)
	+if(p<=1 and q<=1 and isconstant(c),if(hasnot(c,x),C_2+int(c,t,p)-x^q/q!,C_2+t^p/p!-int(1/c,x,q))*if(hasnot(c,t),gsolution(-1/c,y,x,q),0),0) 
	)));
pdsolve(a_*y_+y(p_,t_)+c_*y(q_,x_),y_):= if(hasnot(a,y),if(hasnot(a,c,y),
	if(p<0 or q<0 and isconstant(a,c), C_1*exp((-a/2)^(1/p)*t+(-a/2/c)^(1/q)*x),
	if(isconstant(a) and hasnot(c,t),gsolution(-a,y,t,p)+gsolution(-a/c,y,x,q),
	if(p<=1 and q<=1 and isconstant(a),if(hasnot(c,x),C_2+int(c,t,p)-x^q/q!,C_2+t^p/p!-int(1/c,x,q)),1)*if(hasnot(a,x),gsolution(-a,y,t,p),0)
	+if(p<=1 and q<=1 and isconstant(c),if(hasnot(c,x),C_2+int(c,t,p)-x^q/q!,C_2+t^p/p!-int(1/c,x,q))*if(hasnot(a/c,t),gsolution(-a/c,y,x,q),0),0) 
	))));
pdsolve(a_*y_+y(p_,t_)+y(q_,x_),y_):= if(hasnot(a,y),if(isconstant(a), 
	if(p<0 or q<0, C_1*exp((-a/2)^(1/p)*t+(-a/2)^(1/q)*x), gsolution(-a,y,t,p)+gsolution(-a,y,x,q)),
	 (C_2+t^p/p!-x^q/q!)*if(hasnot(a,x),gsolution(-a,y,t,p),0)
	+(C_2+t^p/p!-x^q/q!)*if(hasnot(a,t),gsolution(-a,y,x,q),0) ));
pdsolve(y(p_,t_)+y(p_,x_)+c_*y_,y_):= if(hasnot(c,y),if(isconstant(c), 
	if(p>0,gsolution(-c,y,x,p)+gsolution(-c,y,t,p),
	C_1*exp((-c/2)^(1/p)*t+(-c/2)^(1/p)*x)+if(p>0,C_2*exp((-c)^(1/p)*t),0)+if(p>0,C_2*exp((-c)^(1/p)*x),0) ),
	if(hasnot(c,x),C_1*E(p,-int(c,t,p)*p!+C_2+t^p-x^p),
	if(hasnot(c,t),C_1*E(p,-int(c,x,p)*p!+C_2+t^p-x^p),
	C_1*E(p,C_2+psolution(1,p,0,c,y,t,x,p)*p!)
	))));
pdsolve(y(p_,t_)+y(q_,x_)+y_,y_):= if(p<0 or q<0, C_1*exp((-1/2)^(1/p)*t+(-1/2)^(1/q)*x), gsolution(-1,y,t,p)+gsolution(-1,y,x,q) );

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

pdsolve(c_*y_^3+y(1,t_)+y(2,x_),y_):= sqrt(-2/c)*(2C_1*x+C_2)/(C_1*x^2+C_2*x-6C_1*t+C_3);
pdsolve(c_*y_^3+y(1,t_)-y(2,x_),y_):= sqrt(2/c)*(2C_1*x+C_2)/(C_1*x^2+C_2*x+6C_1*t+C_3);
pdsolve(a_*log(y_)*y_+y(1,t_)+y(2,x_),y_):= exp(C_1*exp(-a*t)-(x+C_1)^2/4/a+1/2);
pdsolve(a_*log(y_)*y_+y(1,t_)+b_*y(2,x_),y_):= if(isconstant(b), exp(C_1*exp(-a*t)-(x+C_1)^2/4/a/b+1/2) );
pdsolve(log(y_)*y_+y(1,t_)+y(2,x_),y_):= exp(C_1*exp(-t)+1/2-(x+C_1)^2/4);
pdsolve(log(y_)*y_+y(1,t_)+b_*y(2,x_),y_):= if(isconstant(b), exp(C_1*exp(-t)+1/2-(x+C_1)^2/4/b) );

pdsolve(y(1,t_)+y(1,x_)+1/(x_+y_),y_):= W(exp(x) *C_1 *(t - x)) - x + 1;
pdsolve(y(1,t_)+y(1,x_)+c_/(x_+y_),y_):= if(isconstant(c), c*W(exp(x/c) *C_1 *(t - x)) - x+c );
pdsolve(y(1,t_)*x_-t_*y(1,x_)+b_*y_,y_):=if(isconstant(b), C_1*exp(b*atan(x/t))*(t*t+x*x) );
pdsolve(y(1,t_)*x_-t_*y(1,x_)+y_,y_):= C_1*exp(atan(x/t))*(t*t+x*x);
pdsolve(c_*y_+y(1,t_)+b_*y_*y(1,x_),y_):= if(isconstant(b,c), W(b*c*exp(C_1-c*t-c^2*x))/b/c );
pdsolve(c_*y_+y(1,t_)+y_*y(1,x_),y_):= if(isconstant(c), W(c*exp(C_1-c*t-c^2*x))/c );
pdsolve(y_+y(1,t_)+b_*y_*y(1,x_),y_):= if(isconstant(b), W(b*exp(C_1-t-x))/b );
pdsolve(y_+y(1,t_)+y_*y(1,x_),y_):= W(exp(C_1-t-x));

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

pdsolve(y(p_,t_)+y(m_,x_)+y(q_,x_),y_):= gsolution(-1,m,0,y,x,q)+gsolution(0,y,t,p)+C_3*(t^p/p!-x^m/m!);
pdsolve(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!) );
pdsolve(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))) );
pdsolve(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))) ));
pdsolve(y(p_,t_)+y_*y(q_,x_)+b_*y(m_,x_),y_):= if(q<>m and hasnot(b,x) and hasnot(b,y), C_1+C_2*(x*m!-b*t^p/p!) );
pdsolve(y(p_,t_)+y_*y(q_,x_)+b_*x_^n_*y(m_,x_),y_):= if(n<0 and q> m-n and hasnot(b,x) and hasnot(b,y), C_1+C_2*(int(x^(-n),x,m)-b*t^p/p!), dsolve(y*y(q,x)+b*x^n*y(m,x),y) );
pdsolve(y(p_,t_)+y_*y(q_,x_)+x_^n_*y(m_,x_),y_):= if(n<0 and q> m-n, C_1+C_2*(int(x^(-n),x,m)-t^p/p!), dsolve(y*y(q,x)+x^n*y(m,x),y) );
pdsolve(y(p_,t_)+y_*y(q_,x_)+y(m_,x_),y_):= if(m==q,-1-fallingfactorial(-p,p)*(C_1+t)^(-p)*(C_2+x)^m/m!,C_1+C_2*(x^m/m!-t^p/p!) );
pdsolve(y(p_,t_)+y(q_,x_)+y_*y(m_,x_),y_):= -fallingfactorial(-p,p)*(C_1+t)^(-p)*(C_2+x)^m/m!;
pdsolve(y(p_,t_)+y(q_,x_)+b_*y_*y(m_,x_),y_):= if(isconstant(b),-fallingfactorial(-p,p)*(C_1+t)^(-p)*(C_2+x)^m/(m-1)!/b,
	if(hasnot(b,x), (-1)^(1+p)*(C_2+x)^m/(m-1)!/(C_1+int(b,t,p)) ));
pdsolve(y(p_,t_)+a_*y(q_,x_)+y_*y(m_,x_),y_):= if(m==q,-a,0)-fallingfactorial(-p,p)*(C_1+t)^(-p)*(C_2+x)^m/m!;
pdsolve(y(p_,t_)+a_*y(q_,x_)+b_*y_*y(m_,x_),y_):= if(isconstant(b),if(m==q,-a/b,0)-fallingfactorial(-p,p)*(C_1+t)^(-p)*(C_2+x)^m/(m-1)!/b,
	if(hasnot(b,x), (-1)^(1+p)*(C_2+x)^m/(m-1)!/(C_1+int(b,t,p)) ));

pdsolve(y(p_,t_)+y(q_,x_)+b_*y(1,x_),y_):= if(hasnot(b,x) and hasnot(b,y), C_1*(x-b*t^p/p!)+C_3*exp((-b)^(1/(q-1))*x)+C_2 and C_2+C_1*mittag(p,(1/4)*t^p*b^2)*exp(-b*x/2) );
pdsolve(y(p_,t_)+a_*y(q_,x_)+y(1,x_),y_):= if(hasnot(a,x) and hasnot(a,y), C_1*(x-t^p/p!)+C_3*exp((-1/a)^(1/(q-1))*x)+C_2 and C_2+C_1*mittag(p,(1/4)*1/a*t^p)*exp(-x/2/a) );
pdsolve(y(p_,t_)+a_*y(q_,x_)+b_*y(1,x_),y_):= if(hasnot(a,b,x) and hasnot(a,b,y), C_1*(x-b*t^p/p!)+C_3*exp((-b/a)^(1/(q-1))*x)+C_2 and C_2+C_1*mittag(p,(1/4)*1/a*t^p*b^2)*exp(-b*x/2/a) );
pdsolve(y(p_,t_)+a_*y(q_,x_)+b_*x_^n_*y(1,x_),y_):= if(isconstant(a,b), gsolution(-b/a*x^n,1,0,y,x,q)+gsolution(0,y,t,p)
	+C_3*(t^p/p!-if(a+b<>0, 1/(a+b)*x^(1-n)/(1-n), 1/a*((-1/(1-n)^2)*x^(1-n)+1/(1-n)*log(x)*x^(1-n))))
	+if(p==1 and a==b, C_4*Ei(x^(1-n)/((1-n)*(a+b)*t)),0) );

pdsolve(y(p_,t_)+y(q_,x_)+c_*y(q_,x_),y_):= pde(y(p,t)+(1+c)*y(q,x),y);

pdsolve(y(1,t_)+y(3,x_)+asin(y_)*y(1,x_),y_):= sin(-1/(C_1+t)^2+(C_2+x)/(C_1+t));
pdsolve(y(1,t_)+y(3,x_)+aa_*asin(y_)*y(1,x_),y_):= if(isconstant(aa),sin(-1/aa/(C_1+aa*t)^2+(C_2+x)/(C_1+aa*t)) );
pdsolve(y(1,t_)+b_*y(3,x_)+asin(y_)*y(1,x_),y_):= if(isconstant(b),sin(-b/(C_1+t)^2+(C_2+x)/(C_1+t)) );
pdsolve(y(1,t_)+b_*y(3,x_)+aa_*asin(y_)*y(1,x_),y_):= if(isconstant(b,aa),sin(-b/aa/(C_1+aa*t)^2+(C_2+x)/(C_1+aa*t)) );
pdsolve(y(1,t_)+y(3,x_)+acos(y_)*y(1,x_),y_):= cos(-1/(C_1+t)^2+(C_2+x)/(C_1+t));
pdsolve(y(1,t_)+y(3,x_)+aa_*acos(y_)*y(1,x_),y_):= if(isconstant(aa),cos(-1/aa/(C_1+aa*t)^2+(C_2+x)/(C_1+aa*t)) );
pdsolve(y(1,t_)+b_*y(3,x_)+acos(y_)*y(1,x_),y_):= if(isconstant(b),cos(-b/(C_1+t)^2+(C_2+x)/(C_1+t)) );
pdsolve(y(1,t_)+b_*y(3,x_)+aa_*acos(y_)*y(1,x_),y_):= if(isconstant(b,aa),cos(-b/aa/(C_1+aa*t)^2+(C_2+x)/(C_1+aa*t)) );
pdsolve(y(1,t_)+y(3,x_)+asinh(y_)*y(1,x_),y_):= sinh(1/(C_1+t)^2+(C_2+x)/(C_1+t));
pdsolve(y(1,t_)+y(3,x_)+aa_*asinh(y_)*y(1,x_),y_):= if(isconstant(aa),sinh(1/aa/(C_1+aa*t)^2+(C_2+x)/(C_1+aa*t)) );
pdsolve(y(1,t_)+b_*y(3,x_)+asinh(y_)*y(1,x_),y_):= if(isconstant(b),sinh(b/(C_1+t)^2+(C_2+x)/(C_1+t)) );
pdsolve(y(1,t_)+b_*y(3,x_)+aa_*asinh(y_)*y(1,x_),y_):= if(isconstant(b,aa),sinh(b/aa/(C_1+aa*t)^2+(C_2+x)/(C_1+aa*t)) );
pdsolve(y(1,t_)+y(3,x_)+acosh(y_)*y(1,x_),y_):= cosh(1/(C_1+t)^2+(C_2+x)/(C_1+t));
pdsolve(y(1,t_)+y(3,x_)+aa_*acosh(y_)*y(1,x_),y_):= if(isconstant(aa),cosh(1/aa/(C_1+aa*t)^2+(C_2+x)/(C_1+aa*t)) );
pdsolve(y(1,t_)+b_*y(3,x_)+acosh(y_)*y(1,x_),y_):= if(isconstant(b),cosh(b/(C_1+t)^2+(C_2+x)/(C_1+t)) );
pdsolve(y(1,t_)+b_*y(3,x_)+aa_*acosh(y_)*y(1,x_),y_):= if(isconstant(b,aa),cosh(b/aa/(C_1+aa*t)^2+(C_2+x)/(C_1+aa*t)) );
pdsolve(y(1,t_)+y(3,x_)+log(y_)*y(1,x_),y_):= exp(1/(C_1+t)^2+(C_2+x)/(C_1+t));
pdsolve(y(1,t_)+y(3,x_)+c_*log(y_)*y(1,x_),y_):= if(isconstant(c),exp(1/c/(C_1+c*t)^2+(C_2+x)/(C_1+c*t)) );
pdsolve(y(1,t_)+b_*y(3,x_)+log(y_)*y(1,x_),y_):= if(isconstant(b), exp(b/(C_1+t)^2+(C_2+x)/(C_1+t)) );
pdsolve(y(1,t_)+b_*y(3,x_)+c_*log(y_)*y(1,x_),y_):= if(isconstant(b,c),exp(b/c/(C_1+c*t)^2+(C_2+x)/(C_1+c*t)) );
#pdsolve(y(1,t_)+y(3,x_)+c_*y_^2*y(1,x_),y_):= C_1+C_2^2/(sqrt(4C_1^2+C_2^2)*cosh(C_2*x-(c*C_1^2*C_2+C_2^3)*t+C_3)+2C_1);
#pdsolve(y(1,t_)+y(3,x_)+c_*y_^2*y(1,x_),y_):= if(isconstant(c), C_1-4C_1/(4C_1^2*(x-c*C_1^2*t)^2+1) );
#pdsolve(y(1,t_)+y(3,x_)+y_^2*y(1,x_),y_):= C_1-4C_1/(4C_1^2*(x-C_1^2*t)^2+1);
pdsolve(y(1,t_)+y(3,x_)+y_^2*y(1,x_),y_):=  sqrt(-6)* C_1 *tanh(2 C_1^3 *t + C_1 *x + C_2);
pdsolve(y(1,t_)+y(3,x_)+c_*y_^2*y(1,x_),y_):=  if(isconstant(c), sqrt(-6/c)* C_1 *tanh(2 C_1^3 *t + C_1*x + C_2) );
pdsolve(y(1,t_)+b_*y(3,x_)+y_^2*y(1,x_),y_):=  if(isconstant(b), sqrt(-6sgn(b))*abs(b)^(1/6)* C_1 *tanh(2 C_1^3 *t + C_1/cbrt(b) *x + C_2) );
pdsolve(y(1,t_)+b_*y(3,x_)+c_*y_^2*y(1,x_),y_):=  if(isconstant(b,c), sqrt(-6sgn(b)/c)*abs(b)^(1/6)* C_1 *tanh(2 C_1^3 *t + C_1/cbrt(b) *x + C_2) );
#pdsolve(y(1,t_)+y(3,x_)+c_*y_*y(1,x_),y_):= if(isconstant(c), -12/c* C_1^2/(C_2* t + C_1* x + C_3)^2 + C_1^2/c - 2C_2/C_1/c );
#pdsolve(y(1,t_)+y(3,x_)+y_*y(1,x_),y_):= -12 C_1^2/(C_2* t + C_1* x + C_3)^2 + C_1^2 - 2C_2/C_1 and (C_1+x)/(C_2+t);
#pdsolve(y(1,t_)+b_*y_*y(1,x)+c_*y(3,x),y_):=  if(c>0, 6c/b*sech(C_1+sqrt(1/2)/2*(x-c/2*t))^2, 6c/b*sec(C_1+sqrt(1/2)/2*(x-c/2*t))^2 );
#pdsolve(y(1,t_)+b_*y_^2*y(1,x)+c_*y(3,x),y_):= if(c>0, sqrt(6*c/b)*sech(C_1+x-c*t), sqrt(6*c/b)*sec(C_1+x-c*t) );
pdsolve(y(1,t_)+y_*y(1,x)+y(3,x),y_):=  -12 C_1^2 * tanh(C_2 *t + C_1 *x + C_3)**2 + 8 C_1^2- C_2/C_1;
pdsolve(y(1,t_)+y_*y(1,x)+d_*y(3,x),y_):=  if(inconstant(d), -12 C_1^2 *d* tanh(C_2 *t + C_1 *x + C_3)**2 + 8 C_1^2 *d - C_2/C_1 );
pdsolve(y(1,t_)+b_*y_*y(1,x)+y(3,x),y_):=  if(inconstant(b), -(12 C_1^3* tanh(C_2 *t + C_1 *x + C_3)**2 - 8 C_1^3 + C_2)/(b* C_1) );
pdsolve(y(1,t_)+b_*y_*y(1,x)+d_*y(3,x),y_):= if(inconstant(b,d), -(12 C_1^3 *d* tanh(C_2 *t + C_1 *x + C_3)**2 - 8 C_1^3 *d + C_2)/(b* C_1) );

pdsolve(y(1,t_)+b_*y(2,x_)+c_*y(1,x_),y_):= if(has(b,y), if(hasnot(c,x), int(1/int(1/b,y),y)=C_1-t-int(c,t)+x and C_2+C_1*(int(c,t)-x), 
	dsolve(b*y(2,x)+c*y(1,x),y) ),
	if(has(c,y),pde(y(1,t)+c*y(1,x),y)  ));
pdsolve(y(1,t_)+b_*y(2,x_)+y(1,x_),y_):= if(has(b,y), int(1/int(1/b,y),y)=C_1-2t+x and C_2+C_1*(t-x) );
pdsolve(y(1,t_)+a_*b_*y(2,x_)+y(1,x_),y_):= if(hasnot(a,x) and has(b,y), int(1/int(1/b,y),y)=C_1-int(a+1,t)+x and C_2+C_1*(t-x),
	if(hasnot(b,x) and has(a,y), int(1/int(1/a,y),y)=C_1-int(b+1,t)+x and C_2+C_1*(t-x) ));
pdsolve(y(1,t_)+a_*b_*y(2,x_)+c_*y(1,x_),y_):= if(hasnot(a,c,x) and has(b,y), int(1/int(1/b,y),y)=C_1-int(a+c,t)+x and C_2+C_1*(int(c,t)-x),
	if(hasnot(b,c,x) and has(a,y), int(1/int(1/a,y),y)=C_1-int(b+c,t)+x and C_2+C_1*(int(c,t)-x) ));
pdsolve(y(1,t_)+a_*x_*y(2,x_)+a_*y(1,x_),y_):= if(isconstant(a), C_1+C_2*log(x)+C_3*(t-1/a*x)+C_4*Ei(x/(a*t)) );
pdsolve(y(1,t_)+a_*y(2,x_)+b_*y_*y(1,x_),y_):= if(hasnot(a,b,x), 2a/b*C_1/(C_3+C_2*t+C_1*x)-C_2/C_1/b and 2a/b*C_1*tanh(C_3+C_2*t+C_1*x)-C_2/C_1/b );
pdsolve(y(1,t_)+a_*y(2,x_)+y_*y(1,x_),y_):= if(hasnot(a,x), 2a*C_1/(C_1*x+C_2*t+C_3)-C_2/C_1 and 2a*C_1*tanh(C_3+C_2*t+C_1*x)-C_2/C_1 );
pdsolve(y(1,t_)+y(2,x_)+b_*y_*y(1,x_),y_):= if(hasnot(b,x), 2C_1/b/(C_3+C_2*t+C_1*x)-C_2/C_1/b and 2/b*C_1*tanh(C_3+C_2*t+C_1*x)-C_2/C_1/b );
pdsolve(y(1,t_)+y(2,x_)+y_*y(1,x_),y_):= 2*C_1/(C_3+C_2*t+C_1*x)-C_2/C_1 and 2*C_1*tanh(C_3+C_2*t+C_1*x)-C_2/C_1;


#4;
pdsolve(a_+c_*y_^(-2)+y(1,t_)^2+y(1,x_)^2,y_):= if(isconstant(a,c), i*sqrt(c/a+a*(t+C_1)^2+a*(x+C_2)^2) ); 
pdsolve(a_+y_^(-2)+y(1,t_)^2+y(1,x_)^2,y_):= if(isconstant(a), i*sqrt(1/a+a*(t+C_1)^2+a*(x+C_2)^2) ); 

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

pdsolve(y(p_,t_)+y(1,x_)^2+b_*y(q_,x_)+y_,y_):=   -(C_1+x+exp(-t)/((2)^p-1)/(-1)^p)*exp(-t) and -(C_1-x+exp(-t)/((2)^p-1)/(-1)^p)*exp(-t);
pdsolve(y(p_,t_)+y(1,x_)^2+b_*y(q_,x_)+c_*y_,y_):=if(isconstant(c), -(C_1+x+exp(-c*t)/((2)^p-1)/(-c)^p)*exp(-c*t)/a and -(C_1-x+exp(-c*t)/((2)^p-1)/(-c)^p)*exp(-c*t)/a );
pdsolve(y(p_,t_)+a_*y(1,x_)^2+b_*y(q_,x_)+y_,y_):=if(isconstant(a), -(C_1+x+exp(-t)/((2)^p-1)/(-1)^p)*exp(-t)/a and -(C_1-x+exp(-t)/((2)^p-1)/(-1)^p)*exp(-t)/a );
pdsolve(y(p_,t_)+a_*y(1,x_)^2+y(q_,x_)+c_*y_,y_):=if(isconstant(c), -(C_1+x+exp(-c*t)/((2)^p-1)/(-c)^p)*exp(-c*t)/a and -(C_1-x+exp(-c*t)/((2)^p-1)/(-c)^p)*exp(-c*t)/a );
pdsolve(y(p_,t_)+a_*y(1,x_)^2+b_*y(q_,x_)+c_*y_,y_):=if(isconstant(c), -(C_1+x+exp(-c*t)/((2)^p-1)/(-c)^p)*exp(-c*t)/a and -(C_1-x+exp(-c*t)/((2)^p-1)/(-c)^p)*exp(-c*t)/a );

pdsolve(a_*y(1,t_)^2+y_*y(2,t_)+a_*y(1,x_)^2+y_*y(2,x_),y_):= if(isconstant(a), replace(dsolve(a*y(1,x)^2+y*y(2,x),y),x,t+x) );
pdsolve(a_*y(1,t_)^2+b_*y_*y(2,t_)+c_*y(1,x_)^2+d_*y_*y(2,x_),y_):= if(a*d==b*c, replace(dsolve(a*y(1,x)^2+b*y*y(2,x),y),x,t+x) );

#pdsolve(y(p_,t_)+y(q_,x_)+b_*y_+d_,y_):= if(hasnot(d,y), if(p<0 or q<0 and isconstant(b), C_1*exp((-b/2)^(1/p)*t+(-b/2)^(1/q)*x),
	if(p<=1 and q<=1,C_2-t^p/p!+x^q/q!,1)*if(hasnot(b,x),gsolution(-b,y,t,p),0)
	+if(q<=1,C_2+t^p/p!-x^q/q!,1)*if(hasnot(b,t),gsolution(-b,y,x,q),0)
	+psolution(1,q,b,d,y,t,x,p) )); 
#pdsolve(y(p_,t_)+a_*y(q_,x_)+y_+d_,y_):= if(hasnot(a,d,y), if(p<0 or q<0 and isconstant(a), C_1*exp((-1/2)^(1/p)*t+(-1/2/a)^(1/q)*x),
	if(p<=1 and q<=1,C_2+if(has(a,t),int(a,t,p),t^p/p!)-if(hasnot(a,t),int(1/a,x,q),x^q/q!),1)*gsolution(-1,y,t,p)
	+if(q<=1,C_2+t^p/p!-int(1/a,x,q),1)*if(hasnot(a,t),gsolution(-1/a,y,x,q),0)
	+psolution(a,q,1,d,y,t,x,p) )); 
#pdsolve(y(p_,t_)+a_*y(q_,x_)+b_*y_+d_,y_):=if(hasnot(a,b,d,y), if(p<0 or q<0 and isconstant(a,b), C_1*exp((-b/2)^(1/p)*t+(-b/2/a)^(1/q)*x),
	if(p<=1 and q<=1,C_2-if(has(a,t),int(a,t,p),t^p/p!)+if(hasnot(a,t),int(1/a,x,q),x^q/q!),1)*if(hasnot(b,x),gsolution(-b,y,t,p),0)
	+if(q<=1,C_2+t^p/p!-int(1/a,x,q),1)*if(hasnot(a,b,t),gsolution(-b/a,y,x,q),0)
	+psolution(a,q,b,d,y,t,x,p) )); 
pdsolve(y(p_,t_)+y(q_,x_)+y_+d_,y_):= if(hasnot(d,y), pdsolve(y(p,t)+y(q,x)+y,y)+psolution(1,q,1,d,y,t,x,p) ); 
pdsolve(y(p_,t_)+a_*y(q_,x_)+y_+d_,y_):=if(hasnot(a,d,y), pdsolve(y(p,t)+a*y(q,x)+y,y)+psolution(a,q,1,d,y,t,x,p) ); 
pdsolve(y(p_,t_)+y(q_,x_)+b_*y_+d_,y_):=if(hasnot(b,d,y), pdsolve(y(p,t)+y(q,x)+b*y,y)+psolution(1,q,b,d,y,t,x,p) ); 
pdsolve(y(p_,t_)+a_*y(q_,x_)+b_*y_+d_,y_):=if(hasnot(a,b,d,y), pdsolve(y(p,t)+a*y(q,x)+b*y,y)+psolution(a,q,b,d,y,t,x,p) ); 
pdsolve(a_*y(p_,t_)+b_*y(q_,x_)+c_*y_+d_,y_):=if(hasnot(a,b,c,d,y), pdsolve(a*y(p,t)+b*y(q,x)+c*y,y)+psolution(a,p,b,q,c,d,y,x,t) ); 
	
pdsolve(a_*y_+b_*y_+y(p_,t_)+y(q_,x_),y_):=if(hasnot(a,x),gsolution(-a,y,t,p),1)*if(hasnot(b,t),gsolution(-b,y,x,q),gsolution(-b,y,t,p))*(t^p/p!-x^q/q!);
pdsolve(a_*y_+b_*y_+y(p_,t_)+d_*y(q_,x_),y_):=if(hasnot(a,x),gsolution(-a,y,t,p),1)*if(hasnot(b,t),gsolution(-b/d,y,x,q),gsolution(-b,y,t,p))*if(hasnot(d,t),t^p/p!-int(1/d,x,q),1);

pdsolve(y(2,t_)-y(2,x_)+c_*y_^m_+d_*y_^n_,y_):= if(2m-1==n, (-1/4*c*(1-m)^2*((t+C_1)^2-(x+C_2)^2)-d/c/m)^(1/(1-m)) );
#pdsolve(y(2,t_)-y(2,x_)+c_*y_+d_*y_^n_,y_):= (2d/c/(1+n)*sinh(1/2*sqrt(-c)*(1-n)*(sinh(C_1)*x+cosh(C_1)*t)+C_2)^2)^(1/(1-n));

pdsolve(y(1,t_)+y(2,x_)+exp(y_)+d_,y_):= if(isconstant(d), -2log(sqrt(-1/d)+C_1*exp(sqrt(d/2)*x+d*t/2)) );
pdsolve(y(1,t_)+y(2,x_)+b_*exp(y_)+d_,y_):= if(isconstant(b,d), -2log(sqrt(-b/d)+C_1*exp(sqrt(d/2)*x+d*t/2)) );
pdsolve(y(1,t_)+a_*y(2,x_)+exp(y_)+d_,y_):= if(isconstant(a,d), -2log(sqrt(-1/d)+C_1*exp(sqrt(d/a/2)*x+d*t/2)) );
pdsolve(y(1,t_)+a_*y(2,x_)+b_*exp(y_)+d_,y_):= if(isconstant(a,b,d), -2log(sqrt(-b/d)+C_1*exp(sqrt(d/a/2)*x+d*t/2)) );
pdsolve(y(1,t_)+y(2,x_)+exp(c_+y_)+d_,y_):= if(isconstant(d) and d(c,x,2)==0, -c-2log(sqrt(-1/d)+C_1*exp(sqrt(d/2)*x+d*t/2)) );
pdsolve(y(1,t_)+y(2,x_)+b_*exp(c_+y_)+d_,y_):= if(isconstant(b,d) and d(c,x,2)==0, -c-2log(sqrt(-b/d)+C_1*exp(sqrt(d/2)*x+d*t/2)) );
pdsolve(y(1,t_)+a_*y(2,x_)+exp(c_+y_)+d_,y_):= if(isconstant(a,d) and d(c,x,2)==0,-c-2log(sqrt(-1/d)+C_1*exp(sqrt(d/a/2)*x+d*t/2)) );
pdsolve(y(1,t_)+a_*y(2,x_)+b_*exp(c_+y_)+d_,y_):= if(isconstant(a,b,d) and d(c,x,2)==0, -c-2log(sqrt(-b/d)+C_1*exp(sqrt(d/a/2)*x+d*t/2)) );
pdsolve(y(1,t_)+b_*y(2,x_)-y_^2+y_^3,y_):= if(isconstant(b),1+1/(exp(1/sqrt(-2b)*x+1/2t)*C_1-1) );
pdsolve(y(1,t_)+b_*y(2,x_)+y_^2-y_^3,y_):= if(isconstant(b),1+1/(exp(1/sqrt(2b)*x-1/2t)*C_1-1) );
pdsolve(y(1,t_)+b_*y(2,x_)+c_*y_^2+d_*y_^3,y_):= if(c+d==0, 1+1/(exp(sqrt(c/b/2)*x+d/2t)*C_1-1) );
pdsolve(y(1,t_)+b_*y(2,x_)+y_+d_*y_^n_,y_):= if(hasnot(b,d,t) and hasnot(b,d,x), (sqrt(-d)+C_1*exp(-(1-n)*(3+n)/2/(1+n)*t+sqrt(1/b*(1-n)^2/2/(1+n))*x))^(2/(1-n)) );
pdsolve(y(1,t_)+b_*y(2,x_)+c_*y_+y_^n_,y_):= if(hasnot(b,c,t) and hasnot(b,c,x), (sqrt(-1/c)+C_1*exp(-c*(1-n)*(3+n)/2/(1+n)*t+sqrt(c/b*(1-n)^2/2/(1+n))*x))^(2/(1-n)) );
pdsolve(y(1,t_)+b_*y(2,x_)+c_*y_+d_*y_^n_,y_):= if(hasnot(b,c,d,t) and hasnot(b,c,d,x), (sqrt(-d/c)+C_1*exp(-c*(1-n)*(3+n)/2/(1+n)*t+sqrt(c/b*(1-n)^2/2/(1+n))*x))^(2/(1-n)) );
pdsolve(y(1,t_)+b_*y(2,x_)+c_*y_+y_^2,y_):= if(isconstant(b,c), -c/4* (tanh(C_1 + 5c* t/12 + sqrt(c/6/b)*x/2) - 1)^2 );
pdsolve(y(1,t_)+b_*y(2,x_)+c_*y_+c_*y_^3,y_):= if(isconstant(b,c), 1/2i*tanh(C_1 + 3c* t/4 + sqrt(2c/b)/4 x)-1/2i );
pdsolve(y(1,t_)+b_*y(2,x_)+c_*y_+y_^3,y_):= if(isconstant(b,c), sqrt(-c)/2*(tanh(C_1 + 3c* t/4 + sqrt(2c/b)/4*x) - 1) );
pdsolve(y(1,t_)+b_*y(2,x_)+y_+y_^3,y_):= if(isconstant(b), 1/2i*tanh(C_1 + 3 t/4 + sqrt(2/b)/4 x)-1/2i );
pdsolve(y(1,t_)+b_*y(2,x_)+y_+y_^2,y_):= if(isconstant(b),-1/4*(tanh(C_1 + 5 t/12 + x/(2sqrt(6b))) - 1)^2 );
pdsolve(y(1,t_)+y(2,x_)+c_*y_^2+d_*y_^3,y_):= if(c+d==0, 1+1/(exp(sqrt(c/2)*x+d/2t)*C_1-1) );
pdsolve(y(1,t_)+y(2,x_)+y_^2+y_^3,y_):= -1/2 - 1/2tanh(C_1 + t/4 + sqrt(-2)/4*x);
pdsolve(y(1,t_)+y(2,x_)+y_+d_*y_^n_,y_):= if(hasnot(d,t) and hasnot(d,x), (sqrt(-d)+C_1*exp(-(1-n)*(3+n)/2/(1+n)*t+sqrt((1-n)^2/2/(1+n))*x))^(2/(1-n)) );
pdsolve(y(1,t_)+y(2,x_)+c_*y_+y_^n_,y_):= if(hasnot(c,t) and hasnot(c,x), (sqrt(-1/c)+C_1*exp(c*(n-1)*(3+n)/2/(1+n)*t+sqrt(c*(1-n)^2/2/(1+n))*x))^(2/(1-n)) );
pdsolve(y(1,t_)+y(2,x_)+c_*y_+d_*y_^n_,y_):= if(hasnot(c,d,t) and hasnot(c,d,x), (sqrt(-d/c)+C_1*exp(c*(n-1)*(3+n)/2/(1+n)*t+sqrt(c*(1-n)^2/2/(1+n))*x))^(2/(1-n)) );
pdsolve(y(1,t_)+y(2,x_)+c_*y_+y_^3,y_):= if(isconstant(c), sqrt(-c)/2*(tanh(C_1 + 3c* t/4 + sqrt(2c)/4*x) - 1) );
pdsolve(y(1,t_)+y(2,x_)+c_*y_+c_*y_^3,y_):= if(isconstant(c), 1/2i*tanh(C_1 + 3c* t/4 + sqrt(2c)/4 x)-1/2i );
pdsolve(y(1,t_)+y(2,x_)+y_+y_^3,y_):= 1/2i*tanh(C_1 + 3 t/4 + sqrt(2)/4 x)-1/2i;
pdsolve(y(1,t_)+y(2,x_)+y_+y_^2,y_):=-1/4*(tanh(C_1 + 5 t/12 + x/(2 sqrt(6))) - 1)^2;
pdsolve(y(1,t_)+y(2,x_)+c_*y_+y_^2,y_):= if(isconstant(c), -c/4* (tanh(C_1 + 5c* t/12 + sqrt(c)*x/(2 sqrt(6))) - 1)^2 );

pdsolve(y(p_,t_)+b_*y(p1_,x_)+c_*y(p2_,u_)+d_,y_):= if(hasnot(b,c,d,y), psolution(b,p1,c,p2,0,d,y,t,x,u,p)+C_1+C_2*(2*t^p/p!-x^p1/p1!/b-u^p2/p2!/c) );
pdsolve(y(p_,t_)+y(n_,x_)+y(m_,x_)+d_,y_):= if(hasnot(d,y), C_3*(1/p!*t^p-1/m!*x^m)+gsolution(-1,m,0,y,x,n)
	+psolution(1,n,1,m,0,d,y,t,x,p) );
pdsolve(y(p_,t_)+a_*y(n_,x_)+y(m_,x_)+d_,y_):= if(hasnot(a,d,y), C_3*(1/p!*t^p-1/m!*x^m)+gsolution(-1/a,m,0,y,x,n)
	+psolution(a,n,1,m,0,d,y,t,x,p) );
pdsolve(y(p_,t_)+y(n_,x_)+b_*y(m_,x_)+d_,y_):= if(hasnot(b,d,y), pdsolve(y(p,t)+y(n,x)+b*y(m,x),y)
	+psolution(1,n,b,m,0,d,y,t,x,p) );
pdsolve(y(p_,t_)+a_*y(n_,x_)+b_*y(m_,x_)+d_,y_):= if(hasnot(a,b,d,y), pdsolve(y(p,t)+a*y(n,x)+b*y(m,x),y)
	+psolution(a,n,b,m,0,d,y,t,x,p) );
pdsolve(f_+y(1,t_)+y_*y(1,x_)+y_*y(1,u_),y_):= C_1+f*(t+u+x)+y+y^2=0;

pdsolve(y(p_,t_)+b_*y(p1_,x_)+c_*y(p2_,u_)+y_,y_):= gsolution(-1,y,t,p)+gsolution(-1/b,y,x,p1)+gsolution(-1/c,y,u,p2);
pdsolve(y(p_,t_)+b_*y(p1_,x_)+c_*y(p2_,u_)+g_*y_,y_):= gsolution(-g,y,t,p)+gsolution(-g/b,y,x,p1)+gsolution(-g/c,y,u,p2);

pdsolve(y(p_,t_)+y(n_,x_)+y(m_,x_)+y_,y_):= if(m>0, gsolution(-1,y,t,p),0)+if(m>0 or isinteger(m),gsolution(-1,m,-1,y,x,n),0)+C_3*exp(x)*E(p,-3t^p);
pdsolve(y(p_,t_)+y(n_,x_)+y(m_,x_)+c_*y_,y_):= if(m>0 or isinteger(m) and hasnot(c,t),gsolution(-1,m,-c,y,x,n),0)
	+if(m>0 and hasnot(c,x),gsolution(-c,y,t,p),0)
	+if(isconstant(c) and p==1,C_3*exp((-2-c)*t+x), if(hasnot(c,x),gsolution(-c,y,t,p)*exp((-1)^(1/(m-n))*x),0 ));
pdsolve(y(p_,t_)+y(n_,x_)+b_*y(m_,x_)+y_,y_):= if(hasnot(b,t),if(m>0,gsolution(-1,y,t,p),0)+if(m>0 or isinteger(m),gsolution(-b,m,-1,y,x,n),0)+if(isconstant(b),C_3*E(p,-(2+b)*t^p)*exp(x),0),
	if(p<=1 and has(b,t),C_3*exp(-int(b,t,p)*(-1)^(1/n)+(-1)^(1/n)*x),0) );
pdsolve(y(p_,t_)+y(n_,x_)+b_*y(m_,x_)+c_*y_,y_):= if(hasnot(b,c,t),if(m>0,gsolution(-c,y,t,p),0)+if(m>0 or isinteger(m),gsolution(-b,m,-c,y,x,n),0)+if(isconstant(b,c),C_3*E(p,-(1+b+c)*t^p)*exp(x),0),
	if(hasnot(c,x),gsolution(-c,y,t,p)*(if(m>0,C_2,0)+if(isconstant(b),exp((-b)^(1/(n-m))*x),0)),0)
	+if(p<=1 and has(b,t),C_3*exp(int(b,t,p)*(c)^m-int(c,x)),0) );
pdsolve(y(p_,t_)+a_*y(n_,x_)+y(m_,x_)+y_,y_):= if(hasnot(a,t),if(m>0,gsolution(-1,y,t,p),0)+if(m>0 or isinteger(m),gsolution(-1/a,m,-1/a,y,x,n),0)+if(isconstant(a),C_3*E(p,-(a+2)*t^p)*exp(x),0),
	gsolution(-1,y,t,p)*(if(m>0,C_2,0)+if(isconstant(a),exp((-a)^(1/(m-n))*x),0))
	+if(p<=1 and has(a,t),C_3*exp((-1)^(m+n)*int(a,t,p)-x),0) );
pdsolve(y(p_,t_)+a_*y(n_,x_)+y(m_,x_)+c_*y_,y_):= if(hasnot(a,c,t),if(m>0,gsolution(-c,y,t,p),0)+if(m>0 or isinteger(m),gsolution(-1/a,m,-c/a,y,x,n),0)+if(isconstant(a,c),C_3*E(p,-(a+1+c)*t^p)*exp(x),0),
	if(hasnot(c,x),gsolution(-c,y,t,p)*(if(m>0,C_2,0)+if(isconstant(a),exp((-a)^(1/(m-n))*x),0)),0)
	+if(p<=1 and has(a,t),C_3*exp((-1)^(m+n)*int(a,t,p)*(c)^n-int(c,x)),0) );
pdsolve(y(p_,t_)+a_*y(n_,x_)+b_*y(m_,x_)+y_,y_):= if(hasnot(a,b,t),if(m>0,gsolution(-1,y,t,p),0)+if(m>0 or isinteger(m),gsolution(-b/a,m,-1/a,y,x,n),0)+if(isconstant(a,b),C_3*E(p,-(a+b+1)*t^p)*exp(x),0),
	if(p<=1 and has(a,t),C_3*exp((-1)^(m+n)*int(a,t,p)*(1/b)^n-int(1/b,x)),0)
	+if(p<=1 and has(b,t),C_3*exp(int(b,t,p)*(1/a)^m-int(1/a,x)),0) );
pdsolve(y(p_,t_)+a_*y(n_,x_)+b_*y(m_,x_)+c_*y_,y_):= if(hasnot(a,b,c,t),if(m>0,gsolution(-c,y,t,p),0)+if(m>0 or isinteger(m),gsolution(-b/a,m,-c/a,y,x,n),0)+if(isconstant(a,b,c),C_3*E(p,-(a+b+c)*t^p)*exp(x),0),
	if(hasnot(c,x),gsolution(-c,y,t,p)*(if(m>0,C_2,0)+if(isconstant(a,b),exp((-b/a)^(1/(n-m))*x),0)),0)
	+if(p<=1 and has(a,t),C_3*exp((-1)^(m+n)*int(a,t,p)*(c/b)^n-int(c/b,x)),0)
	+if(p<=1 and has(b,t),gsolution(-c,y,t,p)*(int(b,t,p)*(c/a)^m-int(c/a,x)),0) );

pdsolve(c_*y_+y(1,t_)+y(2,t_)+y(2,x_),y_):=if(isconstant(c),C_1*cos(1/2*sqrt(4c-1)*x)*exp((-1/2)*t)+C_2*exp((-1/2)*t)*sin(1/2*sqrt(4c-1)*x));
pdsolve(c_*y_+b_*y(1,t_)+y(2,t_)+y(2,x_),y_):=if(isconstant(b,c),C_1*cos(1/2*sqrt(4c-b^2)*x)*exp((-1/2)*b*t)+C_2*exp((-1/2)*b*t)*sin(1/2*sqrt(4c-b^2)*x));
pdsolve(c_*y_+b_*y(1,t_)+y(2,t_)+a_*y(2,x_),y_):=if(isconstant(a,b,c),C_1*cos(1/2*sqrt(4c/a-b^2/a)*x)*exp((-1/2)*b*t)+C_2*exp((-1/2)*b*t)*sin(1/2*sqrt(4c/a-b^2/a)*x));
pdsolve(-y_^4+y(1,t_)+y_^2*y(1,x_)+y_*y(2,x_),y_):= tanh(C_1-t+x) and tan(C_1-t-x);
pdsolve(-y_^4+y(1,t_)-y_^2*y(1,x_)+y_*y(2,x_),y_):= tanh(C_1-t-x) and tan(C_1-t+x);
pdsolve(y_^4+y(1,t_)+y_^2*y(1,x_)-y_*y(2,x_),y_):= tanh(C_1+t-x) and tan(C_1+t+x);
pdsolve(y_^4+y(1,t_)-y_^2*y(1,x_)-y_*y(2,x_),y_):= tanh(C_1+t+x) and tan(C_1+t-x);

pdsolve(y(p_,t_)+b_*y(q_,x_)+c_*y(q_,u_)+d_*y(q_,v_),y_):= C_1+if(isconstant(b,c,d), C_3*(3/p!*t^p-u^q/c/q!-x^q/b/q!-v^q/d/q!+x+u+v),0)
	+if(q==2,C_2*exp((b+c+d)^(1/p)*t)*sin(x)*sin(u)*sin(v),0 );
pdsolve(y(p_,t_)+b_*y(1,x_)+c_*y(1,u_)+d_*y(1,v_),y_):= C_1+C_2*(3*t^p/p!-u/c-x/b-v/d);
pdsolve(y(1,t_)+y_*y(1,x_)+y_*y(1,u_)+y_*y(1,v_),y_):= (C_1+u+v+x)/(C_2+3t);

pdsolve(y(m_,t_)+y(m_,x_)+y(p_,x_)+d_*y(q_,x_),y_):=  replace(dsolve(-1/d*y(p,x),y,x,q),C_1,C_1-t);
pdsolve(y(m_,t_)+y(m_,x_)+c_*y(p_,x_)+y(q_,x_),y_):=  replace(dsolve(-c*y(p,x),y,x,q),C_1,C_1-t);
pdsolve(y(m_,t_)+y(m_,x_)+c_*y(p_,x_)+d_*y(q_,x_),y_):=  replace(dsolve(-c/d*y(p,x),y,x,q),C_1,C_1-t);
pdsolve(y(m_,t_)+b_*y(m_,x_)+y(p_,x_)+d_*y(q_,x_),y_):= if(isconstant(b), replace(dsolve(-1/d*y(p,x),y,x,q),C_1,C_1-b*t) );
pdsolve(y(m_,t_)+b_*y(m_,x_)+c_*y(p_,x_)+y(q_,x_),y_):= if(isconstant(b), replace(dsolve(-c*y(p,x),y,x,q),C_1,C_1-b*t) );
pdsolve(y(m_,t_)+b_*y(m_,x_)+c_*y(p_,x_)+d_*y(q_,x_),y_):= if(isconstant(b), replace(dsolve(-c/d*y(p,x),y,x,q),C_1,C_1-b*t) );

pdsolve(y(m_,t_)+y(1,x_)+b_*y_^2*y(1,x_)+c_*y(3,x_),y_):=if(m<=1, if(c>0, sqrt(6*(c-1)/b)*sech(C_1+sqrt(abs(1-1/c))*(x-c*t**m/m!)), sqrt(6*(c-1)/b)*sec(C_1+sqrt(abs(1-1/c))*(x-c*t**m/m!)) ));
pdsolve(y(m_,t_)+a_*y(1,x_)+y_^2*y(1,x_)+c_*y(3,x_),y_):=if(m<=1, if(c>0, sqrt(6*(c-a))*sech(C_1+sqrt(abs(1-a/c))*(x-c*t**m/m!)), sqrt(6*(c-a))*sec(C_1+sqrt(abs(1-a/c))*(x-c*t**m/m!)) ));
pdsolve(y(m_,t_)+a_*y(1,x_)+b_*y_^2*y(1,x_)+c_*y(3,x_),y_):=if(m<=1, if(c>0, sqrt(6*(c-a)/b)*sech(C_1+sqrt(abs(1-a/c))*(x-c*t**m/m!)), sqrt(6*(c-a)/b)*sec(C_1+sqrt(abs(1-a/c))*(x-c*t**m/m!)) ));
pdsolve(y(m_,t_)+y(1,x_)+y_**2*y(1,x_)+y(3,x_),y_):=if(m<=1, i*sqrt(6)* C_1 *tanh((2 C_1^2 - 1) *C_1* t**m/m! + C_1 *x + C_2) );
pdsolve(y(m_,t_)+y(1,x_)-y_**2*y(1,x_)+y(3,x_),y_):=if(m<=1,   sqrt(6)* C_1 *tanh((2 C_1^2 - 1) *C_1* t**m/m! + C_1 *x + C_2) );
pdsolve(y(m_,t_)+y(1,x_)+y_*y(1,x_)+y(3,x_),y_):=if(m<=1, -12 C_1^2 *tanh(C_2* t**m/m! + C_1* x + C_3)^2 + 8 C_1^2 - 1 - C_2/C_1 );
pdsolve(y(m_,t_)+y(1,x_)+b_*y_*y(1,x_)+d_*y(3,x_),y_):=if(m<=1, -(C_1 + 12 C_1^3 *d* tanh(C_2 *t**m/m! + C_1 *x + C_3)**2 - 8 C_1^3 *d + C_2)/(b* C_1) );
pdsolve(y(m_,t_)+a_*y(1,x_)+y_*y(1,x_)+d_*y(3,x_),y_):=if(m<=1, -(a* C_1 + 12 C_1^3 *d* tanh(C_2 *t**m/m! + C_1 *x + C_3)**2 - 8 C_1^3 *d + C_2)/C_1 );
pdsolve(y(m_,t_)+a_*y(1,x_)+y_*y(1,x_)+y(3,x_),y_):=if(m<=1, -(a* C_1 + 12 C_1^3 * tanh(C_2 *t**m/m! + C_1 *x + C_3)**2 - 8 C_1^3  + C_2)/C_1 );
pdsolve(y(m_,t_)+a_*y(1,x_)+b_*y_*y(1,x_)+y(3,x_),y_):=if(m<=1, -(a* C_1 + 12 C_1^3 * tanh(C_2 *t**m/m! + C_1 *x + C_3)**2 - 8 C_1^3  + C_2)/(b* C_1) );
pdsolve(y(m_,t_)+a_*y(1,x_)+b_*y_*y(1,x_)+d_*y(3,x_),y_):=if(m<=1, -(a* C_1 + 12 C_1^3 *d* tanh(C_2 *t**m/m! + C_1 *x + C_3)**2 - 8 C_1^3 *d + C_2)/(b* C_1) );

pdsolve(y(m_,t_)+y_*y(1,x_)+y(2,x_)+y(3,x_),y_):= if(m<=1, 1/25*(-3 tanh(C_1 *t + C_2 + x/10)**2 + 6 tanh(C_1 *t + C_2 + x/10) - 250 C_1 + 3));
pdsolve(y(m_,t_)+a_*y_*y(1,x_)+b_*y(2,x_)-y(3,x_),y_):= if(m<=1, 1/25/a*(3b^2* tanh(C_1 *t**m/m! + C_2 + b*x/10)**2 + 6b^2* tanh(C_1 *t**m/m! + C_2 +b* x/10) - 250/b* C_1 - 3b^2) );
pdsolve(y(m_,t_)+y(1,x_)+y_*y(2,x_)+y(3,x_),y_):= if(m<=1, 3 C_1 *tanh((2C_1^2-1) *C_1 *t**m/m! + C_1 *x + C_2) );
pdsolve(y(m_,t_)+a_*y(1,x_)+b_*y_*y(2,x_)+c_*y(3,x_),y_):= if(m<=1, 3c/b*C_1 *tanh((2c*C_1^2-a) *C_1 *t**m/m! + C_1 *x + C_2) );


#5;
#pdsolve(y(p_,t_)+a_*y(q_,x_)^n_+b_*y(m_,x_)+c_*y_+f1_,y_):= if(isconstant(a,c), dsolve(y(q,x)^n+b/a*y(m,x)+c/a*y,y)+if(has(b,x),dsolve(-c*y,y,t,p),0)-f1/c );
#pdsolve(y(p_,t_)+y(q_,x_)^n_+b_*y(m_,x_)+c_*y_+f1_,y_):= if(isconstant(c,f1), dsolve(y(q,x)^n+b*y(m,x)+c*y,y)+if(has(b,x),dsolve(-c*y,y,t,p),0)-f1/c );

#pdsolve(y(1,t_)-y(2,x_)+a_*y_+b_*y_^2+y_^3,y_):= if(b+a+1==0, (C_1*exp(1/sqrt(2)*x+(0.5-a)*t)+a*C_2*exp(a/sqrt(2)*x+(a^2/2-a)*t))/(C_1*exp(1/sqrt(2)*x+(0.5-a)*t)+C_2*exp(a/sqrt(2)*x+(a^2/2-a)*t)+C_3) );
pdsolve(y(1,t_)-y(2,x_)+a_*y_+b_*y_^2+c_*y_^3,y_):= if(a+b+c==0, (C_1*exp(sqrt(2)/2*x+(1/2-a)*t))/(C_1*exp(sqrt(2)/2*x+(1/2-a)*t)+C_2) );
pdsolve(y(1,t_)-y(2,x_)+a_*y_+b_*y_^2+y_^3,y_):= if(b+a+1==0, (C_1*exp(sqrt(2)/2*x+(1/2-a)*t))/(C_1*exp(sqrt(2)/2*x+(1/2-a)*t)+C_2) );
pdsolve(y(1,t_)-y(2,x_)+y_+(-2)*y_^2+y_^3,y_):= exp(sqrt(2)/2 x-t/2)/(exp(sqrt(2)/2 x-t/2)+C_1);

pdsolve(y(p_,t_)+b_*y(p1_,x_)+c_*y(p2_,u_)+g_*y_+f1_,y_):= psolution(b,p1,c,p2,g,f1,y,t,x,u,p)
	+gsolution(-g,y,t,p)+gsolution(-g/b,y,x,p1)+gsolution(-g/c,y,u,p2);
	
pdsolve(y(p_,t_)+y(n_,x_)+y(m_,x_)+y_+d_,y_):= pdsolve(y(p,t)+y(n,x)+y(m,x)+y,y)+psolution(1,n,1,m,1,d,y,t,x,p);
pdsolve(y(p_,t_)+y(n_,x_)+y(m_,x_)+c_*y_+d_,y_):= pdsolve(y(p,t)+y(n,x)+y(m,x)+c*y,y)+psolution(1,n,1,m,c,d,y,t,x,p);
pdsolve(y(p_,t_)+y(n_,x_)+b_*y(m_,x_)+y_+d_,y_):= pdsolve(y(p,t)+y(n,x)+b*y(m,x)+y,y)+psolution(1,n,b,m,1,d,y,t,x,p);
pdsolve(y(p_,t_)+y(n_,x_)+b_*y(m_,x_)+c_*y_+d_,y_):= pdsolve(y(p,t)+y(n,x)+b*y(m,x)+c*y,y)+psolution(1,n,b,m,c,d,y,t,x,p);
pdsolve(y(p_,t_)+a_*y(n_,x_)+y(m_,x_)+y_+d_,y_):= pdsolve(y(p,t)+a*y(n,x)+y(m,x)+y,y)+psolution(a,n,1,m,1,d,y,t,x,p);
pdsolve(y(p_,t_)+a_*y(n_,x_)+y(m_,x_)+c_*y_+d_,y_):= pdsolve(y(p,t)+a*y(n,x)+y(m,x)+c*y,y)+psolution(a,n,1,m,c,d,y,t,x,p);
pdsolve(y(p_,t_)+a_*y(n_,x_)+b_*y(m_,x_)+y_+d_,y_):= pdsolve(y(p,t)+a*y(n,x)+b*y(m,x)+y,y)+psolution(a,n,b,m,1,d,y,t,x,p);
pdsolve(y(p_,t_)+a_*y(n_,x_)+b_*y(m_,x_)+c_*y_+d_,y_):= pdsolve(y(p,t)+a*y(n,x)+b*y(m,x)+c*y,y)+psolution(a,n,b,m,c,d,y,t,x,p);
	
pdsolve(a_-y_^4+y(1,t_)+y_^2*y(1,x_)+y_*y(2,x_),y_):= tanh(C_1-t+x) and tan(C_1-t-x);
pdsolve(a_-y_^4+y(1,t_)-y_^2*y(1,x_)+y_*y(2,x_),y_):= tanh(C_1-t-x) and tan(C_1-t+x);
pdsolve(a_+y_^4+y(1,t_)+y_^2*y(1,x_)-y_*y(2,x_),y_):= tanh(C_1+t-x) and tan(C_1+t+x);
pdsolve(a_+y_^4+y(1,t_)-y_^2*y(1,x_)-y_*y(2,x_),y_):= tanh(C_1+t+x) and tan(C_1+t-x);
	
pdsolve(y(p_,t_)+b_*y(q_,x_)+c_*y(q_,u_)+d_*y(q_,v_)+g_,y_):= if(hasnot(g,y), 
	psolution(b,q,c,q,d,q,0,g,y,t,x,u,v,p)+if(isconstant(b,c,d),C_3*(3/p!*t^p-u^q/c/q!-x^q/b/q!-v^q/d/q!),0) );
pdsolve(y(p_,t_)+b_*y(q_,x_)+c_*y(q_,u_)+d_*y(q_,v_)+y_,y_):= if(v<>x,gsolution(-1,y,t,p)+gsolution(-1/b,y,x,q)+gsolution(-1/c,y,u,q)+gsolution(-1/d,y,v,q) );
pdsolve(y(p_,t_)+b_*y(q_,x_)+c_*y(q_,u_)+d_*y(q_,v_)+g_*y_,y_):= if(v<>x,gsolution(-g,y,t,p)+gsolution(-g/b,y,x,q)+gsolution(-g/c,y,u,q)+gsolution(-g/d,y,v,q) );
pdsolve(y(p_,t_)+b_*y(1,x_)+c_*y(1,u_)+d_*y(1,v_)+g_*y_,y_):= if(v<>x,gsolution(-g,y,t,p)+gsolution(-g/b,y,x,1)+gsolution(-g/c,y,u,1)+gsolution(-g/d,y,v,1) );
pdsolve(y(q_,x_)+b_*y(p_,x_)+y(q2_,z_)+d_*y(p2_,z_)+f_*y_,y_):= dsolve(-b*y(p,x)-f*y,y,x,q)+dsolve(-b*y(p2,z)-f*y,y,z,q2);
pdsolve(y(p_,t_)+a_+b_*y_+c_*y(q_,x_)+y(q2_,x_),y_):=if(hasnot(a,b,c,t), dsolve(-a-b*y-c*y(q,x),y,x,q2)+C_3*E(p,-b*t^p) );
pdsolve(f_+y(1,t_)+y_*y(1,x_)+y_*y(1,u_)+y_*y(1,v_),y_):= C_1+f*(t+u+v+x)+y+3/2*y^2=0;

pdsolve(y(2,t_)+b_*y(2,x_)+c_*y(1,t_)+d_*y(1,x_)+g_*y_,y_):= gsolution(-c,-g,y,t,2)*(1+exp(-d/b*x))+gsolution(-d/b,-g/b,y,x,2)*(1+exp(-c*t));
pdsolve(y(2,x_)+b_*y(2,t_)+c_*y(1,x_)+d_*y(1,t_)+g_*y_,y_):= gsolution(-c,-g,y,x,2)*(1+exp(-d/b*t))+gsolution(-d/b,-g/b,y,t,2)*(1+exp(-c*x));
pdsolve(y(2,x_)+b_*y(2,t_)+c_*y(1,x_)+d_*y(1,t_)+y_,y_):= gsolution(-c,-1,y,x,2)*(1+exp(-d/b*t))+gsolution(-d/b,-1/b,y,t,2)*(1+exp(-c*x));
pdsolve(y(2,x_)+y(2,t_)+c_*y(1,x_)+d_*y(1,t_)+g_*y_,y_):= gsolution(-c,-g,y,x,2)*(1+exp(-d*t))+gsolution(-d,-g,y,t,2)*(1+exp(-c*x));
pdsolve(y(2,x_)+y(2,t_)+c_*y(1,x_)+d_*y(1,t_)+y_,y_):= gsolution(-c,-1,y,x,2)*(1+exp(-d*t))+gsolution(-d,-1,y,t,2)*(1+exp(-c*x));


#6;
pdsolve(y(q_,x_)+b_*y(p_,x_)+y(q2_,z_)+d_*y(p2_,z_)+f_*y_+g_,y_):= dsolve(-b*y(p,x)-f*y,y,x,q)+dsolve(-b*y(p2,z)-f*y,y,z,q2)-g/f;
pdsolve(y(p_,t_)+b_*y(q_,x_)+c_*y(q_,u_)+d_*y(q_,v_)+y_+f1_,y_):= psolution(b,q,c,q,d,q,1,f1,y,t,x,u,v,p)
	+gsolution(-1,y,t,p)+gsolution(-1/b,y,x,q)+gsolution(-1/c,y,u,q)+gsolution(-1/d,y,v,q);
pdsolve(y(p_,t_)+b_*y(q_,x_)+c_*y(q_,u_)+d_*y(q_,v_)+g_*y_+f1_,y_):= psolution(b,q,c,q,d,q,g,f1,y,t,x,u,v,p)
	+gsolution(-g,y,t,p)+gsolution(-g/b,y,x,q)+gsolution(-g/c,y,u,q)+gsolution(-g/d,y,v,q);

pdsolve(y(m_,t_)+y(q_,x_)+b_*y(p_,x_)+y(q2_,z_)+d_*y(p2_,z_)+f_*y_,y_):= 
	if(hasnot(f,g,t), dsolve(-b*y(p,x)-f*y,y,x,q)+dsolve(-d*y(p2,z)-f*y,y,z,q2)+C_3*E(m,-f*t^m) );
pdsolve(y(m_,t_)+a_*y(q_,x_)+b_*y(p_,x_)+c_*y(q2_,z_)+d_*y(p2_,z_)+f_*y_,y_):= 
	if(hasnot(f,g,t), dsolve(-b/a*y(p,x)-f/a*y,y,x,q)+dsolve(-d/c*y(p2,z)-f/c*y,y,z,q2)+C_3*E(m,-f*t^m) );

pdsolve(y_-y_**2-y_^3+y(1,t_)+y_*y(1,x_)+y(2,x_),y_):= tanh(C_1-t+x);


#7;
pdsolve(y(m_,t_)+y(q_,x_)+b_*y(p_,x_)+y(q2_,z_)+d_*y(p2_,z_)+f_*y_+g_,y_):= 
	if(hasnot(f,g,t), dsolve(-b*y(p,x)-f*y,y,x,q)+dsolve(-d*y(p2,z)-f*y,y,z,q2)+C_3*E(m,-f*t^m)-g/f );
pdsolve(y(m_,t_)+a_*y(q_,x_)+b_*y(p_,x_)+c_*y(q2_,z_)+d_*y(p2_,z_)+f_*y_+g_,y_):= 
	if(hasnot(f,g,t), dsolve(-b/a*y(p,x)-f/a*y,y,x,q)+dsolve(-d/c*y(p2,z)-f/c*y,y,z,q2)+C_3*E(m,-f*t^m)-g/f );


#=;
pdsolve(a_=b_,y_):=pdsolve(replace(toy(a-b),y(x),y),y);
pdsolve(y(p_,t_)=a_+b_*y_,y_):=if(hasnot(a,b,t), pdsolve(a+b*y)+C_3*exp(b*t),pdsolve(y(1,t)-a-b*y));
pdsolve(y(p_,t_)=a_+b_*y_+y(q_,x_),y_):=if(hasnot(a,b,t), dsolve(-a-b*y,y,x,q)+C_3*E(p,b*t^p) );
pdsolve(y(p_,t_)=a_+b_*y_+c_*y(q_,x_),y_):=if(hasnot(a,b,c,t), dsolve(-a/c-b/c*y,y,x,q)+C_3*E(p,b*t^p) );


#0;
pdsolve(y(p_,t_)+k1_*y(q_,x_)+a1_*y_+b1_*z,z(p_,t_)+k2_*z(q_,x_)+a2_*y_+b2_*z):=
	block(x1:=-(a1+b2)/2+sqrt((a1+b2)^2-4*a1*b2+4*a2*b1)/2,x2:=-(a1+b2)/2-sqrt((a1-b2)^2+4*a2*b1)/2,
	if(isconstant(k1,k2) and p>0 and q>0,(y1:=dsolve(x1*y,y,t,p)+dsolve(x1/k1*y,y,x,q),y2:=dsolve(x2*y,y,t,p)+dsolve(x2/k2*y,y,x,q)),
	(y1:=pde(y(p,t)+k1*y(q,x)-x1*y,y),y2:=pde(y(p,t)+k2*y(q,x)-x2*y,y)) ),
	y=((a1+x2)*y1-(a1+x1)*y2)/((x1-x2)*a2) and z=(y1-y2)/(x1-x2) );
pdsolve(y(p_,t_)+y(q_,x_)+a1_*y_+b1_*z,z(p_,t_)+z(q_,x_)+a2_*y_+b2_*z):=
	block(x1:=-(a1+b2)/2+sqrt((a1+b2)^2-4*a1*b2+4*a2*b1)/2,x2:=-(a1+b2)/2-sqrt((a1-b2)^2+4*a2*b1)/2,
	if(p>0 and q>0,(y1:=dsolve(x1*y,y,t,p)+dsolve(x1*y,y,x,q),y2:=dsolve(x2*y,y,t,p)+dsolve(x2*y,y,x,q)),
	(y1:=pde(y(p,t)+y(q,x)-x1*y,y),y2:=pde(y(p,t)+y(q,x)-x2*y,y)) ),
	y=(a1+x2)/a2/(x1-x2)*y1-(a1+x1)/a2/(x1-x2)*y2 and z=(y1-y2)/(x1-x2) );
pdsolve(y(p_,t_)+y(q_,x_)+y(q2_,x_)+a1_*y_+b1_*z,z(p_,t_)+z(q_,x_)+z(q2_,x_)+a2_*y_+b2_*z):=
	block(x1:=-(a1+b2)/2+sqrt((a1+b2)^2-4*a1*b2+4*a2*b1)/2,x2:=-(a1+b2)/2-sqrt((a1-b2)^2+4*a2*b1)/2,
	if(p>0 and q>0,(y1:=dsolve(x1*y,y,t,p)+dsolve(x1*y,y,x,q)+dsolve(x1*y,y,x,q2),y2:=dsolve(x2*y,y,t,p)+dsolve(x2*y,y,x,q)+dsolve(x2*y,y,x,q2)),
	(y1:=pde(y(p,t)+y(q,x)+y(q2,x)-x1*y,y),y2:=pde(y(p,t)+y(q,x)+y(q2,x)-x2*y,y)) ),
	y=(a1+x2)/a2/(x1-x2)*y1-(a1+x1)/a2/(x1-x2)*y2 and z=(y1-y2)/(x1-x2) );


pdsolve(a_):=pdsolve(replace(toy(a),y(x),y),y);
pdsolve(a_=b_):=pdsolve(toy(a-b),y);
pdsolve(y(p_,t)=b_):=if(p<1 and hasnot(b,t),block(f:=pdsolve(y(1,t)-toy(b),y),if(hasnot(f,C_1+t) and hasnot(f,(-C_1-t)),replace(f,t,t**p/p!),pdsolve(y(p,t)-toy(b),y))),
	pdsolve(y(p,t)-toy(b),y) );
pdsolve(ds(y_,t,p_)=b_):=if(p<1 and hasnot(b,t),block(f:=pdsolve(y(1,t)-toy(b),y),if(hasnot(f,C_1+t) and hasnot(f,(-C_1-t)),replace(f,t,t**p/p!),pdsolve(y(p,t)-toy(b),y))),
	pdsolve(y(p,t)-toy(b),y) );
pdsolve(a_=b_,c_=d_):=pdsolve(toy(a-b),toy(c-d));