y(x):=sqrt(1+cos(4*x));

'integrate(y(x),x,0,%pi/4);

integrate(y(x),x,0,%pi/4);

float(%);