Write matlab code given information below.It is finite difference method for partial diffential equation in a semicurcilar boundary.

Figure 1 shows a semicircle of radius r=1.2 meters. Find the temperature distribution over the surface by solving V u=0 numerically in polar coordinate Use Ar=10 cm and 40= n/8. u=80 .C r=1.2 m u=0 'C Figure 1: Semicircle plate with zero temperature at its bottom and 80 C at its circular boundary u(i, j) = u(i + 1, j) + u(i - 1,j) lu(i + 1,j) - u(i - 1,j) 2 2 (Ar) 2 + (Ar) 2 + 2Ar (r;) 2 (40 ) 2 ri + 1 u(i,j + 1) + u(i,j -1), (r;)2 (40)2How can Apply? A single cell + Withj- Uij + 1=0 ( Ar ) ? 2 Ar ( 612 ) 2 In this case, i=0, 1,2,3, 12, \\quad j=0,1,2, 3,...8 i=0,1,2,3 412,j=0,1,2,3,...8 , the point u(0,0) u(0,0) is the center of the semicircle. The boundary conditions are as follows, for \\quad 1=0,11, \\quad u(i, 0)=0, \\quad j=0,....8, \\quad u(12,j) =80 for i=0,...11,u(i,0)=0,j=0,.....8,u(12,j)=80 for \\quad j=0 \\quad \\theta =0, \\quad for \\quad j=8 \\quad \\theta =\\pi