2

I am using a Matlab function called besselh which gives an exact solution to Theoderon's lift deficiency function for a two-dimensional airfoil executing a simple harmonic motion in incompressible flow with a given range of k.

I am trying to plot the real part of this function with the given range of k. I feel I am missing something really simple in my code. Perhaps, I might be using the for loop wrong by including the "plot" function within the loop?

for k=10^-10:0.2:1.5 %range of k
H=(besselh(1,2,k))/(besselh(1,2,k)+i*besselh(0,2,k)) %Bessel function
plot(k,real(H))
%plot(k,imag(H))
end

Below is a picture of how the plot should look like (solid line).

Any help would be appreciated, thank you in advance!

This is how the plot should look like

Jesus
  • 179
  • 2
  • 2
  • 13

1 Answers1

2

You are trying to plot each point individually. Just calculate the function across the whole range for k and then plot:

k = 10^-10:0.2:1.5; %range of k
% notice the ./ means element-wise vector division
H = (besselh(1,2,k)) ./ (besselh(1,2,k) + i * besselh(0,2,k)); %Bessel function
plot(k,real(H))
StaticBeagle
  • 5,070
  • 2
  • 23
  • 34