radiation pattern of horn antenna

hi everybody

i am trying to do this integral in matlab please can anyone help me with that

this is the code that i did and it seems working but i am not getting a good result

thanks

x = linspace(-A/2,A/2,numberMultipathComponent);

fun= @(x) (cos((pi.*x)./A) .* exp(-j*k*(sqrt(R0+(x.^2))-(x.*sin(thetarad)))));

fethetahor = integral(fun,(-A/2),(A/2),'ArrayValued',1);

[David Goodmanson]

Hi khaled, first thing is, R0 needs to be squared in your function

[khaled whi]

thanks for your answer David

indeed, R0 needs to be squared