Advertisement
Not a member of Pastebin yet?
Sign Up,
it unlocks many cool features!
- ###Определение пользовательской функции f(x):
- f <- function(x) # Заголовок функции
- {(exp(x)-exp(-x))/2} # Тело функции
- ###Построение графика функции f(x):
- dx <- -0.001 # Задание дельта-отступа от точки разрыва
- ## Задание диапазонов точек независимой переменной x
- # Значение аргумента x на левой ветви графика
- x1<-seq(-1000,length = 100001)
- # Значение аргумента x на правой ветви графика
- x2 <- seq(-1000,length = 100001)
- y1 <- f(x1) #Значение функции f(x) на левой части графика
- y2 <- f(x2) # Значение функции f(x) на правой части графика
- ## Вычерчивание левой части графика оранжевым цветом
- plot(x1,y1,type = "l", ylab = "y", xlab = "x",
- xlim=c(-5,8),ylim=c(-30,30),# Интервалы вывода графика
- col="orange", lwd =2, lty = "solid",
- main = "Исследование функции f(x)") #Заголовок графика
- ##Вычерчивание правой части графика желтым цветом
- lines(x2,y2,type = "l",col = "yellow",lwd=2)
- ##Нанесение на график осей,вертикальной ассимптоты и линий сетки
- # Задание точек для построения на графике вертикальных линий сетки
- p<-seq(-2,8,0.5)# Шаг сетки 0.5
- #вывод штриховых вертикальных линий сетки красного цвета
- abline (v=p,col = "red", lty = "dashed")
- #Задание точек для построения на графике горизонтальных линий сетки
- p<-seq(-30,30,10)# Шаг сетки 5
- #вывод штриховых вертикальных линий сетки зеленого цвета
- abline (v=p,col = "green", lty = "dashed")
- #нанесение на график осей ox и oy серого цвета
- abline (h=0,v=0,col="gray50" ,lty = "solid")
- ###Нахождение нулей функции f(x)
- ## Задание точности определения нулей фукнции
- dex<-.0000001
- #Интервал поиска 1-го(левого) корня
- m<--2.2;n<-1.6
- ## Вычисление и вывод значений функции на границах левого интерала
- z1<-f(m);z2<-f(n)
- z1;z2 ## Знаки разные
- ## Сохранение полной информации о найденном 1-ом корне
- vv1 <- uniroot(f,c(m,n),tol=dex)
- vv1 ## Вывод полной информации о найденном 1-ом корне
- # интервала поиска 2-го(правого) корня
- m1<- -1.4;n1<-0.2
- ## Вычисление и вывод значений функции на границах правого интерала
- z3<-f(m);z4<-f(n)
- z3;z4
- ## Сохранение полной информации о найденном 2-ом корне
- vv2 <- uniroot(f,c(m1,n1),tol=dex)
- vv2 ## Вывод полной информации о найденном 2-ом корне
- ## Выделение нулей функции на графике
- points(c(vv1[1],vv2[1]),c(vv2[2],vv2[2]),
- col ="red",pch=21,bg="black")
Advertisement
Add Comment
Please, Sign In to add comment
Advertisement