SharedMATHPHYS / bessel functions.sagewsOpen in CoCalc
%auto
# This cell automatically evaluates on startup -- or run it manually if it didn't evaluate.
# Here, it initializes the Jupyter kernel with the specified name and sets it as the default mode for this worksheet.
%default_mode jupyter_kernel

from numpy import arange
def gpl(func, n=1, k=0, b=20, max=None, min=None):
G = Graphics()
for i in range(8):
G += plot(func(n*(i + k), x), (x, 0, b), hue=sin(i/10), aspect_ratio='automatic', ymax=max, ymin=min)
return G

Error in lines 1-1 Traceback (most recent call last): File "/cocalc/lib/python2.7/site-packages/smc_sagews/sage_server.py", line 1013, in execute exec compile(block+'\n', '', 'single') in namespace, locals File "", line 1, in <module> File "/cocalc/lib/python2.7/site-packages/smc_sagews/sage_jupyter.py", line 32, in __call__ raise ValueError("You may not run Sage kernels from a Sage worksheet.\nInstead use the sage_select command in a Terminal to\nswitch to a different version of Sage, then restart your project.") ValueError: You may not run Sage kernels from a Sage worksheet. Instead use the sage_select command in a Terminal to switch to a different version of Sage, then restart your project.

# ЦЕЛОГО ПОЛОЖИТЕЛЬНОГО порядка

func = bessel_J
gpl(func)


# ЦЕЛОГО ОТРИЦАТЕЛЬНОГО порядка

gpl(func, n=-1)


# ПОЛОВИННОГО ПОЛОЖИТЕЛЬНОГО порядка

gpl(func, k=1/2)


# ПОЛОВИННОГО ОТРИЦАТЕЛЬНОГО порядка

gpl(func, n=-1, k=1/2, min=-50, max=50)


################################

# ЦЕЛОГО ПОЛОЖИТЕЛЬНОГО порядка

func = bessel_Y
gpl(func, min=-10)


# ЦЕЛОГО ОТРИЦАТЕЛЬНОГО порядка

gpl(func, n=-1, max=10, min=-10)


# ПОЛОВИННОГО ПОЛОЖИТЕЛЬНОГО порядка

gpl(func, k=1/2, min=-10)


# ПОЛОВИННОГО ОТРИЦАТЕЛЬНОГО порядка

gpl(func, n=-1, k=1/2)


# ЦЕЛОГО ПОЛОЖИТЕЛЬНОГО порядка

func = bessel_I
gpl(func)


# ЦЕЛОГО ОТРИЦАТЕЛЬНОГО порядка

gpl(func, n=-1)


# ПОЛОВИННОГО ПОЛОЖИТЕЛЬНОГО порядка

gpl(func, k=1/2)


# ПОЛОВИННОГО ОТРИЦАТЕЛЬНОГО порядка

gpl(func, n=-1, k=1/2, b=5, min=-50, max=50)


# ЦЕЛОГО ПОЛОЖИТЕЛЬНОГО порядка

func = bessel_K
gpl(func, max=10)


# ЦЕЛОГО ОТРИЦАТЕЛЬНОГО порядка

gpl(func, n=-1, max=10)


# ПОЛОВИННОГО ПОЛОЖИТЕЛЬНОГО порядка

gpl(func, k=1/2, max=10)


# ПОЛОВИННОГО ОТРИЦАТЕЛЬНОГО порядка

gpl(func, n=-1, k=1/2, max=50)


# ЦЕЛОГО ПОЛОЖИТЕЛЬНОГО порядка

func = spherical_bessel_J
gpl(func)


# ЦЕЛОГО ОТРИЦАТЕЛЬНОГО порядка

gpl(func, n=-1, max=10, min=-10)


# ПОЛОВИННОГО ПОЛОЖИТЕЛЬНОГО порядка

gpl(func, k=1/2)


# ПОЛОВИННОГО ОТРИЦАТЕЛЬНОГО порядка

gpl(func, n=-1, k=1/2, max=5)


# ЦЕЛОГО ПОЛОЖИТЕЛЬНОГО порядка

func = spherical_bessel_Y
gpl(func, min=-10)


# ЦЕЛОГО ОТРИЦАТЕЛЬНОГО порядка

gpl(func, n=-1, min=-2)


# ПОЛОВИННОГО ПОЛОЖИТЕЛЬНОГО порядка

gpl(func, k=1/2, min=-1)


# ПОЛОВИННОГО ОТРИЦАТЕЛЬНОГО порядка

gpl(func, n=-1, k=1/2, max=5, min=-5)