channelflow.org

This page is read only. You can view the source, but not change it. Ask your administrator if you think this is wrong.
====== Math 445 lab 5: programming basics ======

**Problem 1:** Write a **for** loop that will print statements of the form

<code>
1 times pi is 3.1416
2 times pi is 6.2832
</code>

for the numbers 1 through 10.

----

**Problem 2:** If $F$ is a temperature in Farenheit, then $C = \frac{5}{9}(F-32)$ is the same temperature in Celsius. Write a function ''farenheit2celsius'' that takes a Celsius temperature as input, converts it to Celsius, prints a statement of the form ''20 Farenheit is -6.6667 Celsius.'', and returns the Celsius value as its output. 

----

**Problem 3:** Write a **for** loop that use the above function to print a list of statements 

<code>
-10 Farenheit is -23.3333 Celsius.
0 Farenheit is -17.7778 Celsius.
</code>

for Farenheit temperatures from -10 to 100, in steps of 10.

----

**Problem 4:** Write a function ''average'' that takes an input vector $x$ and computes the average (mean) value of its components, according to the formula

\begin{equation*}
\text{average}(x) = \frac{1}{N} \sum_{i=1}^N x_i
\end{equation*}

Test your function by comparing its output to the output of the built-in Matlab function **mean** on a random vector of length 100 whose elements are uniformly distributed between 0 and 1.

----

**Problem 5:** The formula for matrix-vector multiplication $y = Ax$ is 

<latex>
y_i = \sum_{j=1}^n A_{ij} x_j
</latex>

In this formula, $A$ is an $m \times n$ matrix, $x$ is an $n$-dimensional column vector, and $y$ is an $m$-dimesional column vector.

Write a function ''matvecmult'' that takes a matrix $A$ and a vector $x$ as arguments, computes $y = Ax$ according to that formula, and returns the vector $y$. Compare your ''matvecmult'' to Matlab's built-in matrix-vector multiplication operator ''*'' on a random $4 \times 4$ matrix and a random 4d column vector.

...to be continued...
channelflow.org

User Tools

Site Tools

Page Tools
