Please someone help me correct my code. I want to plot a 2D graph for distinct alpha in one single 2D plot. See the attached photo for the guidience

alpha=[0.7, 0.8,0.9,1];
P=1;
t=0.5;
s=[-5:0.2:5];
A=(1-s.^2)+((P-4).*(t.^alpha))./(gamma(alpha+1));
Plot(s,A,'-r');
xlabel('x')
ylabel('u')
legend()

Answer 1

[Walter Roberson]

alpha=[0.7, 0.8,0.9,1];
P=1;
t=0.5;
s = (-5:0.2:5).';
A=(1-s.^2)+((P-4).*(t.^alpha))./(gamma(alpha+1));
plot(s,A);
xlabel('x')
ylabel('u')
legend( "\alpha = " + alpha)

Answer 2

[Omorodion Solomon]

Thank you so much. Lastly, Please how can I apply same to plot this?

[Walter Roberson]

For finite cases:

m = (1:2:2*n-3).';
urt = r + sum(prod(m.^2)./r.^(2*n-1) .* t.^(alpha .* n) ./ gamma(alpha .* n + 1),1)

for infinite upper bound you would need to use symsum() and symprod()