1. Function, switch statement. Write a function that takes two square matrices of the same size and an operator name as inputs and apply the operator on the two matrices. The operator name can be ‘+’, ‘-‘, ‘*’, and ‘/’. For example, when the function opName input is '+', the function calculates A+B; when the function opName input is '-', the function calculates A-B.

Function C = ApplyOp(A, B, opName)

Switch opName

Case ‘+’

C=…;

…

In your main script, create the following matrices and call your function with all operators.

A=421 B = 340

020 715

109 095

(2) Function. You work for an engineering design firm that has been contracted to evaluate the concentration of pollutants. A tall chimney releases (20 kg/s) of smoke into a steady uniform wind of speed (u( m/s)). An expression for the decay of the peak pollutant concentration C(ppm) with distance downwind is given as:

C(x,y,z)= M / 4πρDx * e^−((y^2+z^2)u / 4Dx)

Assuming that the diffusion coefficient D to be 0.01 m2/s, and the density of air to be

ρ=1.225 kg/m^3

• Write a function which gets the location(x,y,z) and the wind speed, u, as inputs and calculates the pollutant concentration.
• Call your function to find how pollution concentration at location (x=3000, y=0, z=10) changes with wind speed 0<u<12 m/s. Plot the pollution concentration vs wind speed.
• question (3) Loops and flow control. Make function called loopTest(N) that loops through the values 1 through N and for each number n it should display ‘n is divisible by 2’, ‘n is divisible by 3’, ‘n is divisible by 2 AND 3’ or ‘n is NOT divisible by 2 or 3’. Use a for loop, the function mod or rem to figure out if a number is divisible by 2 or 3, and num2str to convert each number to a string for displaying. You can use any combination of if, else, and elseif.