Tak to wygląda w pseudokodzie, Przeanalizuj sobie i Przepisz do matlaba, to się Nauczysz sumować szeregi:
def double_fact(n):
m = 1
while n > 0:
m *= n
n -= 2
return m
def real_value(x):
return 1 / (sqrt(1 - ( 16 * x * x)))
def some_series(x, e):
s_sum = 1 # sum
n = 1 # comp index
k = 1 # comp index
while True:
tmp = ( double_fact(k) / (pow(2, n) * factorial(n))) * (pow(4 * x, k + 1))
s_sum += tmp
if abs(s_sum - real_value(x)) < e:
break
n += 1
k += 2
return s_sum
double_fact
: https://en.wikipedia.org/wiki/Double_factorial - to Twoje 1, 1 * 3, 1 * 3 * 5, ...
some_series
liczy wartośc dla x
z zadaną dokładnością e
(u Ciebie 0.01
); Puszczasz to w for
z jakimś skokiem w zakresie (-0.25; 0.25)
i pozamiatane:)