$\%$ close all

clear, clc

$\%\%$

a $=2$; $\%$ the radius for observation

b $=0\mathrm{.}2$; $\%$ the radius of the drill hole

Nr $=50$; $\%$ Nr: Number of r$-$direction grid, you can change it

Np $=50$; $\%$ Nt: Number of phi$-$direction grid, you can change it

J$0=1$;

N $=520$; $\%$ Summation number, you can change it

$\%\%$ Initialize Matrix Allocation

Jr$=$

Jp$=$

Jx$=$

Jy$=$

r$=$

phi$=$

$[$R$,$PHI$]=$

$\%\%$ Main loop $/$ Calculation

$\%\%$ Plot

$\%$ Use quiver function

quiver$($X$,$Y$,$Jx$,$Jy$)$;

Given the skeleton code above how would I solve the question given in the image. Please note a$=2,$ b$=0\mathrm{.}2$ and J$0=1$ and that the code should be suitable for MATLAB

Find current distribution on metal plate with hole

Initial condition: Metal plate

Modified condition: hole

Infinite metal plate

Hole radius $0\mathrm{.}2$ m

Observation radius $2$ m

Origin $=$ hole origin

Infinite Plate

J$0=1$