I don’t have your function to plot, so I created one. The rest of the code is (I hope) relatively straightforward to understand. I leave it to you to insert the correct structure of the lightning arrestor structure, although you may not need it since the mesh plot will obscure it. Change it as necessary to create the result you want.
The Code —
f = @(b,r) b(1).*exp(-b(2).*r); % Function Of Radius Only
a = linspace(0, 2*pi, 50); % Angle
r = linspace(0, 100, 50); % Radius
b = [30 0.05]; % Parameter Vector
[A,R] = meshgrid(a, r); % Creat Mesh Matrices
fr = f(b,R); % Evaluate Function Of ‘R’ Matrix
[xm,ym,zm] = pol2cart(A,R,fr); % Convert To Cartesian
figure(1)
mesh(xm, ym, zm)
grid on
The Plot —
