# Mittag-Leffler function;

mittag(0,b_,x_) := 1/(1-x)/Gamma(b);
mittag(a_,1,x_) := mittag(a,x);
mittag(1,2,x_) := exp(x)-1;
mittag(1,3,x_) := exp(x)-1-x;
mittag(2,2,x_) := cosh(sqrt(x))-1;
mittag(2,2,x_^2) := cosh(x)-1;
mittag(2,2,-x_^2) := cos(x)-1;
mittag(a_,inf,x_) :=0;

mittag(a_,ln(a_,x_)):= x;
mittag(0,x_):= 1/(1-x);
mittag(1,x_):= exp(x);
mittag(2,x_^2):= cosh(x);
mittag(2,-x_^2):= cos(x);
#mittag(4,x_):=cos(x^(1/4))/2+cosh(x^(1/4))/2;

mittag(a_,0):=1;