Табулирование трансцендентных функций

Форум, посвященный общим вопросам, связанным с MATLAB

Модератор: Admin

Belypled
Пользователь
Сообщения: 1
Зарегистрирован: Вт май 15, 2018 10:21 pm

Табулирование трансцендентных функций

Сообщение Belypled » Вт май 15, 2018 10:49 pm

Добрый вечер всем! Уже два дня сижу над задачей, никак не получается решить, помогите пожалуйста.
Дан интегральный косинус:
Изображение
y=0.577215

Нужно протабулировать Ci(x) на отрезке [a, b] с шагом h и точностью eps, основываясь на ряде Тейлора, предварительно вычислив его:
Изображение

где а=0.4, b=4, h=0.2, eps=Изображение

Вот что у меня получилось, не судите строго, только начала работать в матлабе :D

Код: Выделить всё

function q = aa(a, b, h, eps)
q=zeros(1, (b-a)/h);
y=0.577215
v=1
n=1
t=0
I=0
for x=a:h:b
    for n=n:1:inf
    I=(((-1)^n)*(x^(2*n)))/((2*n)*(factorial(2*n)))
    t=t+I   
    end
    q(1, v)=t+log(x)+y
    v=v+1
end
hold on
plot(a:h:b, q(1, :))
end


График строится, но с графиком первой функции не сходится. Буду рада, если предложете свой вариант или укажете на ошибку в моем коде :D