ЛР1 Программирование (1.1-1.2 В21) - IV

LabWork1
📑 IT Online Mini Campus. Task 1.1. Task 1.2

IV. Программа  SageMath (Python)

var('x'); A,B,N=.55,1,10; X1,X2=1000,6000

def F(f,x,a,b): 

    try: return round(f(x),4)

    except: return NaN

    if a>=b: return NaN

def create_table(f,a,b,n):

    d=(b-a)/(n-1); x0=randint(X1,X2)*.1^4    

    t=[[round(x,4),F(f,x,a,b)] for x in srange(a,b+d,d)]

    pretty_print(html(latex(f)))

    pretty_print(html('x=%.4f ↦ f(x)=%s'%(x0,F(f,x0,a,b))))

    show(table([['x','f(x)']]+t))

    return t

def create_plottable(f,a,b,n): 

    t=create_table(f,a,b,n)

    list_plot(t).show(gridlines=True,figsize=4)

f21(x)=(sin(x*pi/2)+x^(1/5))/(sqrt(abs(cos(pi*x)+1))*exp(sqrt(x)))

create_plottable(f21,A,B,N)

import sympy; from sympy import *

sympy.init_printing(use_unicode=True)

x=Symbol('x')

y=(sin(x*pi/2.)+x^(1/5.))/(sqrt(abs(cos(pi*x)+1.))*exp(1.)^sqrt(x));

a,b,n=.55,1,10; d=(b-a)/(n-1)

print('x   -> f(x)'); print(15*'-')

for t in range(n):

    x0=a+t*d

    if cos(pi*x0)!=-1 and x0>0 and a<b:

        pprint('%.2f, %.4f'%(x0,y.evalf(subs={x:x0}))); print(15*'-')

    else:

        pprint('%.2f, %s'%(x0,'NaN'))

plot(y,(x,.55,1),ylim=(0,30));