R Диференциальные уравнения

Вот таким способом я решил уравнение y''-5y'+6y=2e^x

   eq1 <- function(t, y,param) {
return(list(c(
y[2],
5*y[2] - 6*y[1]+ 2*exp(t)    
)))
 }
#################################################
require("ggplot2") # для построения графиков
require("rasterVis")
require("fields")
require("deSolve") # решение дифф. уравнений с начальными условиями
require("bvpSolve") # решение дифф. уравнений с краевыми условиями
require("dplyr") # манипуляции с данными
#
y.start <- c(y1 = 1, y2 = 1) # начальные условия
t <- seq(0, 5, by = 0.01)# ось Х - ов
#
sol <- ode(y = y.start, times = t, func = eq1)
sol <- data.frame(sol)
head(sol)
#
equ <- exp(sol$time)
 ggplot(sol) +
 geom_line(aes(time, y1), size = 1) + labs(x = "t",y = expression(paste(y[1])), 
 title = "Решение ОДУ с начальными условиями") +
 geom_line(aes(time,equ),col="red",lty=2,lwd=1.1)

Есть идеи у кого как решать такое уравнение y'=sinx*cosx В нем нет константных значений и триганометрия. Не до конца представляю себе реализацию но возможно ли каким либо образом переобразовать данный код для решения данного уравнения?