Surface Plot for Basic Reproductive Number of Ross , Macdonald Malaria model

I am trying to replicate a plot for basic reproductive number

here is my code

[mu2,a] =meshgrid(0.01:0.001:0.03,0.005:0.01:0.5);

b =0.5;

c=0.5;

m=20;

r=0.01;

% Ro = (m.*a.^2.*b.*c)./(r.*mu2);

surf(mu2,a,Ro)

The problem is that the graph I get is not exactly what is in the paper, especially the range of a,mu2 and R0, there seem to be a mistake on the R0 formula , I don't know what is wrong. Please help

Answer 1

[Anmol Dhiman]

You can use the following code

[mu2,a] =meshgrid(0.1:0.01:0.5,0.01:0.001:0.03);
b =0.5;
c=0.5;
m=20;
r=0.01;
Ro = (m.*a.^2.*b.*c)./(r.*mu2);
surf(mu2,a,Ro)
hold on
z2 =0*mu2 + 1
surf(mu2,a,z2,'red')

I got this figure using the above code:

[Chukwuma Ugoala]

Thank you anmol,

I have implement the correction you made on the code but I still get this error attached below

[Anmol Dhiman]

Hi update the surf command to below line

surf(mu2,a,z2,'MarkerFaceColor','red')

[Chukwuma Ugoala]

Thank you very much