#!/usr/bin/env python3fromsympyimport symbols, oo, pprint, summation, Limit, plot
print('10.')
n, k = symbols('n, k', integer=True)
an =1/ (n **2+1)
bn = n +2for i, o inenumerate([an, bn]):
pprint(o)
l = Limit(n * (o / o.subs({n: n +1}) -1), n, oo)
s = summation(o, (n, 0, oo))
for s in [l, l.doit(), s]:
pprint(s)
print()
print()
p = plot(*[sum([o.subs({n: n0})
for n0 inrange(m)])
for m inrange(10)],
legend=True,
show=False)
colors = ['red', 'green', 'blue', 'brown', 'orange',
'purple', 'pink', 'gray', 'skyblue', 'yellow']
for o, color inzip(p, colors):
o.line_color = color
p.show()
p.save(f'sample9_{i}.png')
0 コメント:
コメントを投稿