I am trying to solve the functions below numerically using MATLAB
This is my code where I have used "integral" command to solve for the second integral, while I used Reimann sum to solve for the first. Take for example c=1/128.
x= -1/128:0.0001:1/128; for j=1:length(x) den= @(y) pi^2.* sqrt(1- (x(j)+y).^2) .* sqrt(1- y.^2); fun= @(y) 1./den(y); ymin= max(-0.999,-0.999-x(j)); ymax= min(0.999,0.999-x(j)); pdfX(j)=integral(fun,ymin,ymax); end sum1=0; for j=1:length(omegat) func1(j)=pdfX(j)*0.0001; sum1=sum1+func1(j); end F=sum1;
Do you think it is correct? Are there any other ways to make it more accurate?
Thanks
Best Answer