Возможно, вы заметили, что y = sin(x)*cos(x)^2
неотрицательно в интервале 0<=x<=pi/2
, поэтому ваша цель может быть достигнута через интеграл sin(x)*cos(x)^2
за указанный интервал.
Возможно, вы можете попробовать integrate
:
s <- integrate(function(x) sin(x)*cos(x)**2,0,pi/2)$value
такой, что
> s
[1] 0.3333333