MeijerG function

meijerG function

syms x s;
b=1:1:10;
y1=gamma(b-s);
b1=11:1:20;
y2=gamma(1-b1+s);
for i=1:1:9
z1=y1(i).*y1(i+1);
end
for i1=1:1:9
z2=y2(i1).*y2(i1+1);
end
a=1:1:10;
y3=gamma(a-s);
a1=11:1:20;
y4=gamma(1-a1+s);
for j=1:1:9
z3=y3(j).*y3(j+1);
end
for j1=1:1:9
z4=y4(j1).*y4(j+1);
end
%t1=((y1-s).*(y4+s));
%t2=((y2+s).*(y3-s));
t1=y1.*y4;
t2=y2.*y3;
t3=t1./t2;
t4=(t3).*(x.^s);
t5=int(t4,s);
ezplot(t5(1));
grid;