Advertisement
Guest User

Untitled

a guest
Mar 25th, 2019
116
0
Never
Not a member of Pastebin yet? Sign Up, it unlocks many cool features!
text 3.42 KB | None | 0 0
  1. ###Определение пользовательской функции f(x):
  2. f <- function(x) # Заголовок функции
  3. {(exp(x)-exp(-x))/2} # Тело функции
  4. ###Построение графика функции f(x):
  5. dx <- -0.001 # Задание дельта-отступа от точки разрыва
  6. ## Задание диапазонов точек независимой переменной x
  7. # Значение аргумента x на левой ветви графика
  8. x1<-seq(-1000,length = 100001)
  9. # Значение аргумента x на правой ветви графика
  10. x2 <- seq(-1000,length = 100001)
  11. y1 <- f(x1) #Значение функции f(x) на левой части графика
  12. y2 <- f(x2) # Значение функции f(x) на правой части графика
  13. ## Вычерчивание левой части графика оранжевым цветом
  14. plot(x1,y1,type = "l", ylab = "y", xlab = "x",
  15. xlim=c(-5,8),ylim=c(-30,30),# Интервалы вывода графика
  16. col="orange", lwd =2, lty = "solid",
  17. main = "Исследование функции f(x)") #Заголовок графика
  18. ##Вычерчивание правой части графика желтым цветом
  19. lines(x2,y2,type = "l",col = "yellow",lwd=2)
  20. ##Нанесение на график осей,вертикальной ассимптоты и линий сетки
  21. # Задание точек для построения на графике вертикальных линий сетки
  22. p<-seq(-2,8,0.5)# Шаг сетки 0.5
  23. #вывод штриховых вертикальных линий сетки красного цвета
  24. abline (v=p,col = "red", lty = "dashed")
  25. #Задание точек для построения на графике горизонтальных линий сетки
  26. p<-seq(-30,30,10)# Шаг сетки 5
  27. #вывод штриховых вертикальных линий сетки зеленого цвета
  28. abline (v=p,col = "green", lty = "dashed")
  29. #нанесение на график осей ox и oy серого цвета
  30. abline (h=0,v=0,col="gray50" ,lty = "solid")
  31. ###Нахождение нулей функции f(x)
  32. ## Задание точности определения нулей фукнции
  33. dex<-.0000001
  34. #Интервал поиска 1-го(левого) корня
  35. m<--2.2;n<-1.6
  36. ## Вычисление и вывод значений функции на границах левого интерала
  37. z1<-f(m);z2<-f(n)
  38. z1;z2 ## Знаки разные
  39. ## Сохранение полной информации о найденном 1-ом корне
  40. vv1 <- uniroot(f,c(m,n),tol=dex)
  41. vv1 ## Вывод полной информации о найденном 1-ом корне
  42. # интервала поиска 2-го(правого) корня
  43. m1<- -1.4;n1<-0.2
  44. ## Вычисление и вывод значений функции на границах правого интерала
  45. z3<-f(m);z4<-f(n)
  46. z3;z4
  47. ## Сохранение полной информации о найденном 2-ом корне
  48. vv2 <- uniroot(f,c(m1,n1),tol=dex)
  49. vv2 ## Вывод полной информации о найденном 2-ом корне
  50. ## Выделение нулей функции на графике
  51. points(c(vv1[1],vv2[1]),c(vv2[2],vv2[2]),
  52. col ="red",pch=21,bg="black")
Advertisement
Add Comment
Please, Sign In to add comment
Advertisement