Настроение сейчас - Бооринг

function su1(n,a,b)
grstart
setcolor('r')
sy=0
sg=0
for x=a:(b-a)/10:b
y=((log(x)-(2/9))*pr1(n,x))-4
sy=sy+y
sg=sg+1
ca=sy/sg
punct(x,y)
end
setcolor('none')
punct(a,ca)
for x=a:(b-a)/10:b
setcolor('g')
punct(x,ca)
end
kvadrat

-------------------------------------------------------------------

function pr1=pr1(n,x)
pr1=1
for k=1:n
pr1=pr1*(abs(x-k)*cos(((k*x)^(1/3))/2))
end


-------------------------------------------------------------------

function kvadrat=kvadrat
if abs(su1)<0,01
punct(-1;su1)
setcolor('b')
end

setcolor('none')
vector(0.2,0.2)
setcolor('g')
vector(0,-0.4)
vector(-0.4,0)
vector(0,0.4)
vector(0.4,0)
setcolor('none')
vector(-0.2,-0.2)
setcolor('b')
end



------------------------------------------------------------------------

И стало светло!!!!!!


[700x525]