gsolution1(a_,y_,x_,4) := if(hasnot(a,x), exp(a^(1/4)*x) );
gsolution2(a_,y_,x_,4) := if(hasnot(a,x), exp(-a^(1/4)*x) );
#gsolution3(a_,y_,x_,4) := if(hasnot(a,x), exp((-1)^(3/4)*(a)^(1/4) *x)  );
#gsolution4(a_,y_,x_,4) := if(hasnot(a,x), exp(-(-1)^(3/4)*(a)^(1/4)*x)  );
gsolution3(a_,y_,x_,4) := if(hasnot(a,x), cos(a^(1/4)*x)  );
gsolution4(a_,y_,x_,4) := if(hasnot(a,x), sin(a^(1/4)*x)  );

gsolution1(b_,y_,x_,3) := if(hasnot(b,x), exp(b^(1/3)* x) );
gsolution2(b_,y_,x_,3) := if(hasnot(b,x), exp(-(-1)^(1/3)*b^(1/3)* x) );
gsolution3(b_,y_,x_,3) := if(hasnot(b,x), exp((-1)^(2/3)*b^(1/3)* x) );

gsolution1(b_,y_,x_,2) := if(hasnot(b,x), exp(b^(1/2)* x) );
gsolution2(b_,y_,x_,2) := if(hasnot(b,x), exp(-b^(1/2)* x) );