The figure below shows a rectangular domain with prescribed temperature

and heat flux

conditions on its edges.

The theoretical temperature distribution in this plate obtained by solving the steady$-$state heat equation is

In the expression above, note that

is the hyperbolic sine function and you can compute it with numpy.sinh

In practice, this summation is truncated to a finite number of terms in order to perform the calculation; including more terms increases the accuracy of the calculation.

Write a function calc$\_$Texact that takes as inputs the x and y coordinates to calculate the temperature and the number of terms nterms to compute in the summation. x and y are scalars $($floats$)$ and nterms is an integer.

In your code below, the following function will be graded: