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?
Kshitij Singh answered .
2025-11-20
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);