How do I write the following summation for power series in code ?

I need to solve for the above power series and plot the solutions from the above formula against radius.

'a' is the coefficeint in the power series and 'p' is a dimensionless radial coordinate, n= 0 ,1,2, ..., infinity.'q' is a pressure term.

Also, I have found the values of 'a' for a0 , a2 , a4 , so how do I code the above formula such that the specific values of 'a' are called as well as the summation is done?

[John D'Errico]

DON'T create a list of variables named a0 a2 a4, etc. Instead, learn to use vectors.

But once you have done that, surely you know how to write a loop and form a sum?

[EMIL RAJAN VARGHESE]

Thanks this was great help

Answer 1

[Arya Chandan Reddy]

Hi, I understand that you are trying to compute the above function for a specific value of 'p' . Assuming that you have coefficients stored in an array in the format

a = [a0 a2 a4 so on..];

You can try to sum first 100 terms or till the available coefficients.

Here is the code for it

for i = 1:size(a)
    Wp = Wp + a(i)*(1 - p^(2*i +2));
end
Wp = Wp * q^(2/3);

Hope it helps