Kindly help me integrate this function.

I am trying to integrate the following expression using 'quad' fuction for all the values of vector 'h' , kindly help me with the code as the vector 'h' is being problamatic to me. Code starts below:

h=1:0.02:10;

f= @(x) ((sin(x)).^2)./(((sin(x)).^2)+(h));

% a, the lower limit

a= 0 ; % b, the upper limit

b= 34 ; answer=quad(f,a,b);

%%% Code ends

[Star Strider]

I thought we solved this in http://www.mathworks.com/matlabcentral/answers/129173-kindly-help-me-integrate-this-function.

Answer 1

[Eugene]

Depends if you're trying a double integral or just evaluating the integral for various values for h. If its NOT a double integral then:

h=1:0.02:10;
f= @(x,c) sin(x).^2./(sin(x).^2+c);
% a, the lower limit
a= 0 ; % b, the upper limit
b= 34 ; 
for i=1:length(h);
  answer(i)=quad(@(x)f(x,h(i)),a,b);
end
%%%Code ends