hi there. can some fix error in this code.

function G = FI_I(I)

%segma_N=10.^-7;
a = 2;     %Alpha
b = 2;     %Beta
c = 0.2;   %Gamma
d = 0.8;   %Big omega prime
k = 1;
z = 4;      %Propagation Distance
%segma_db = 8;  %Atmospheric attenuation coefficient
segma = (8*0.115129255);

%p_db = 0:1:30; %power
%p = 10.^(p_db/10);
%R = 0.5;    %Optoelectronic conversion factor
s = 0.3;    %Zero boresight 

I1 = exp(-segma*z)*0.11512925; %Atmospheric attenuation

Aperture_radius = 0.1;   %Aperture_radius
segma_s = 0.2; %Jitter standard deviation
omega_z = 2.5;   %Beam Width

v = (sqrt(pi)*Aperture_radius)/(sqrt(2)*omega_z);

omega_zeq = sqrt((omega_z.^2*(sqrt(pi)*erf(v)))/(2*v*exp(-v.^2))); 

g = omega_zeq/(2*segma_s);

A0 = (erf(v)).^2;

A = ((2*a.^(a/2))/(c.^(1+a/2)*gamma(a)))*(((c*b)/(c*b+d)).^(b+a/2));

ak = (nchoosek(b-1,k-1))*(((c*b+d).^(1-k/2))/gamma(k))*((d/c).^(k-1))*((a/b).^(k/2));

%PDF for a non-zero boresight pointing error
part_1 = (((g.^2)*A*exp(-s.^2/(2*segma_s.^2)))./((A0*I1).^(g.^2))).*(I.^(g.^2-1));
Q = integral(@integration,2,10);
G = (part_1)*ak*Q;

end

[Walter Roberson]

integration is not defined in your call to integral()

[rana zain]

https://drive.google.com/drive/folders/1FjJXWrscCk4JVfjeoAsQ-YwJ0OApYK0P?usp=sharing

I have written three matlab file :

1. integration.m file for the integration inside equation 14 .
2. FI_I.m file for the equation 14 and call the integration function (integration.m) inside it.
3. non_zero_o.m file for finally write a code for plotting the equation 16 where I tried to call the function for equation 14

I am trying to plot a figure similar the given Fig.1 for non-zero boresight pointing error (details is given on the pdf) using the equation 14 and 16. And for Zero-boresight pointing error using the equation 13 and 16.

N.B: In zero-boresight (equation 13) there is a function named Meijer for which I am stuck to write the code.