x = 0:50;
y = @(x) exp(-x./25).*sqrt((25+x)/(25-x));
int_y = integral(y,0,50)
x = 0:50;
y = @(x) exp(-x./25).*sqrt((25+x)/(25-x));
int_y = cumtrapz(x,y(x));
figure
plot(x, real(int_y), x,imag(int_y))
grid
legend('Real','Imaginary', 'Location', 'best')
.