use matlab to solve this please:The net heat exchange by radiation from plate $1$ with

radius $b$ to plate $2$ with radius a that are separated by

distance $c$ is given by:

$q=b^2{F}_{1-2}(T_1^4-T_2^4)$

where $T_1$ and $T_2$ are the absolute temperatures of the

plates.

$=5\mathrm{.}669{10}^{-8}\frac{W}{m^2-K^4}$ is the Stefan$-$Boltzman constant, and $F_{1-2}$ is a shape factor

which, for the arrangement in the figure, is given by:

$F_{1-2}=\frac{1}{2}[Z-\sqrt[\phantom{2}]{Z^2-4x^2{Y}^2}]$

where $x=\frac{a}{c}$ and, $Y=\frac{c}{b},$ and $Z=1+(1+x^2)Y^2.$

Write a script file that calculates the heat exchange q$.$ For input the program asks

the user to enter values for $T_1,T_2,a,b,$ and $c.$ For output, the program prints $a$

summary of the geometry and temperatures and then prints the value of $q.$

As an example, use your script to calculate the results for $T_1=400K,T_2=600K,a=1m,$

$b=2m,$ and $c=0\mathrm{.}1,1,$ and $10m$ and displays the results as follows:use matlab:The net heat exchange by radiation from plate $1$ with

radius $b$ to plate $2$ with radius a that are separated by

distance $c$ is given by:

$q=b^2{F}_{1-2}(T_1^4-T_2^4)$

where $T_1$ and $T_2$ are the absolute temperatures of the

plates.

$=5\mathrm{.}669{10}^{-8}\frac{W}{m^2-K^4}$ is the Stefan$-$Boltzman constant, and $F_{1-2}$ is a shape factor

which, for the arrangement in the figure, is given by:

$F_{1-2}=\frac{1}{2}[Z-\sqrt[\phantom{2}]{Z^2-4x^2{Y}^2}]$

where $x=\frac{a}{c}$ and, $Y=\frac{c}{b},$ and $Z=1+(1+x^2)Y^2.$

Write a script file that calculates the heat exchange q$.$ For input the program asks

the user to enter values for $T_1,T_2,a,b,$ and $c.$ For output, the program prints $a$

summary of the geometry and temperatures and then prints the value of $q.$

As an example, use your script to calculate the results for $T_1=400K,T_2=600K,a=1m,$

$b=2m,$ and $c=0\mathrm{.}1,1,$ and $10m$ and displays the results as follows: