Roman Putanowicz Wiki en:teaching:subjects:it:labs

en:teaching:subjects:it:labs:lab1

Anonymous (anonymous@undisclosed.example.com) — 2017-10-02T15:54:45+0200

author=Roman Putanowicz backgroundtext= title=Lab 1: Introduction to GNU/Linux operating system showbuttons=off \input{/home/prac/putanowr/dokuwiki/data/media/wiki/latex/it_preamb.tex} \renewcommand{\chaptername}{Lab Handouts} \setcounter{chapter}{0} \begin{comment}

en:teaching:subjects:it:labs:lab2

Anonymous (anonymous@undisclosed.example.com) — 2017-10-02T15:54:45+0200

author=Roman Putanowicz backgroundtext= title=Lab 2: Getting started with Octave showbuttons=off \input{/home/prac/putanowr/dokuwiki/data/media/wiki/latex/it_preamb.tex} \renewcommand{\chaptername}{Lab Handouts} \setcounter{chapter}{1} \begin{comment} \begin{equation*}\text{.}\quad\quad \sin(a)^{\sin(b)+1} \end{equation*}$\rho_s$

en:teaching:subjects:it:labs:lab3

Anonymous (anonymous@undisclosed.example.com) — 2017-10-02T15:54:45+0200

author=Roman Putanowicz title=Lab 3: Defining and using functions backgroundtext= showbuttons=off \input{/home/prac/putanowr/dokuwiki/data/media/wiki/latex/it_preamb.tex} \setcounter{chapter}{2} \begin{comment} \begin{equation*} f(x) = \frac{1}{1+(x\sin(x))^2} \end{equation*}$\alpha$

en:teaching:subjects:it:labs:lab4

Anonymous (anonymous@undisclosed.example.com) — 2017-10-02T15:54:45+0200

author=Roman Putanowicz backgroundtext= title=Lab 4: Control flow : loops and conditional statements showbuttons=off \input{/home/prac/putanowr/dokuwiki/data/media/wiki/latex/it_preamb.tex} \setcounter{chapter}{3} \begin{comment} $ \sum$$\prod $\begin{equation*} \begin{split} a_i &= sin(2\pi i + 1) \\ b_i &= \frac{1}{i+1} \end{split} \end{equation*}\begin{equation*} \sum_{j=1}^{j=N}\sum_{i…

en:teaching:subjects:it:labs:lab5

Anonymous (anonymous@undisclosed.example.com) — 2017-10-02T15:54:45+0200

author=Roman Putanowicz backgroundtext= title=Lab 5: More on control flow; Octave versus Matlab showbuttons=off \input{/home/prac/putanowr/dokuwiki/data/media/wiki/latex/it_preamb.tex} \setcounter{chapter}{4} \begin{comment} $x_0$\begin{equation*} x_{1+n} = \frac{1}{2}\left( x_n + \frac{A}{x_n}\right) \end{equation*}$\sqrt{A}$$x_n$\begin{equation*} \left|\sqrt{A} - x_n\right| < \varepsilon \en…

en:teaching:subjects:it:labs:lab6

Anonymous (anonymous@undisclosed.example.com) — 2017-10-02T15:54:45+0200

author=Roman Putanowicz backgroundtext= title=Lab 6: Solving problems; more on plotting showbuttons=off \input{/home/prac/putanowr/dokuwiki/data/media/wiki/latex/it_preamb.tex} \setcounter{chapter}{5} \begin{comment} \begin{equation*} \mathbf{r}(t)= \begin{cases} x(t) = 2*cos(t)\\ y(t) = 2*sin(t) \end{cases} \quad \text{for} \quad t\in [0,2\pi] \end{equation*}

en:teaching:subjects:it:labs:lab7

Anonymous (anonymous@undisclosed.example.com) — 2017-10-02T15:54:45+0200

~~UNDERCONSTRUCTION~~ author=Roman Putanowicz backgroundtext= title=Lab 7: Solving problems; operating on files showbuttons=off \input{/home/prac/putanowr/dokuwiki/data/media/wiki/latex/it_preamb.tex} \begin{comment}

en:teaching:subjects:it:labs:lab8

Anonymous (anonymous@undisclosed.example.com) — 2017-10-02T15:54:45+0200

~~UNDERCONSTRUCTION~~ author=Roman Putanowicz backgroundtext= title=Lab 1: Introduction to GNU/Linux operating system showbuttons=off \input{/home/prac/putanowr/dokuwiki/data/media/wiki/latex/it_preamb.tex} \begin{comment}

en:teaching:subjects:it:labs:sol_1_1

Anonymous (anonymous@undisclosed.example.com) — 2017-10-02T15:54:45+0200

back Solution to exercise 1.1 Issue the following shell commands: mkdir projectA mkdir projectA/data projectA/results mkdir -p projectA/src/include mkdir projectA/src/lib In order to create the file 'caseA.' use your favourite editor, e.g: gedit projectA/data/caseA.dat

en:teaching:subjects:it:labs:sol_1_2

Anonymous (anonymous@undisclosed.example.com) — 2017-10-02T15:54:45+0200

back Solution to exercise 1.12 Issue the following commands: cd projectA cp data/caseA.data data/caseB.dat

en:teaching:subjects:it:labs:sol_1_3

Anonymous (anonymous@undisclosed.example.com) — 2017-10-02T15:54:45+0200

back Solution to exercise 1.3 The data file can be downloaded from here: . To unpack it, type: gunzip months.txt.gz There are several solutions to this exercise. One can observe that the months which names end with 'er' are the last four months and can use

en:teaching:subjects:it:labs:sol_2_1

Anonymous (anonymous@undisclosed.example.com) — 2017-10-02T15:54:45+0200

back Solution to exercise 2.1 Run Octave: octave -q Crate variables a and b by assigning them a value: a = 0.7 b = 1.0 Enter the expression. Mind operators priorities and be careful with the braces: sin(a)^(sin(b)+1) Again assign values to variables

en:teaching:subjects:it:labs:sol_2_1_1

Anonymous (anonymous@undisclosed.example.com) — 2017-10-02T15:54:45+0200

back Solution to exercise 2.1.1 Run Octave: octave -q Crate variables a and b by assigning them a value: a = 0.7 b = 1.0 Enter the expression. Mind operators priorities and be careful with the braces: sin(a)^(sin(b)+1) Again assign values to variables

en:teaching:subjects:it:labs:sol_2_1_2

Anonymous (anonymous@undisclosed.example.com) — 2017-10-02T15:54:45+0200

back Solution to exercise 2.1.2 Dot and vector products of two vectors u and v can be caluclated using Octave built-in functions dot(u,v) and cross(u,v) respectively. disp(“---------------------------------------------------------------”

en:teaching:subjects:it:labs:sol_2_1_3

Anonymous (anonymous@undisclosed.example.com) — 2017-10-02T15:54:45+0200

back Solution to exercise 2.1.3 The key to the solution is to calculate the last element of the sequence and then use operator : (colon) to generate the sequence. a = input(“Give sequence initial term : ”) d = input(“Give sequence difference :

en:teaching:subjects:it:labs:sol_2_3

Anonymous (anonymous@undisclosed.example.com) — 2017-10-02T15:54:45+0200

back Solution to exercise 2.3

en:teaching:subjects:it:labs:sol_3_1_1

Anonymous (anonymous@undisclosed.example.com) — 2017-10-02T15:54:45+0200

back Solution to exercise 3.1.1 The implementation is straightforward : the difference of the area of a circle sector and a triangle. function A = csegarea(R, alpha) A = 1/2 * R^2 * alpha - 1/2*R^2*sin(alpha); endfunction However less operations are done if we slightly rearange the terms:

en:teaching:subjects:it:labs:sol_3_1_2

Anonymous (anonymous@undisclosed.example.com) — 2017-10-02T15:54:45+0200

back Solution to exercise 3.1.2 In order to normalize a vector we have to calculate its length. It can be using the function norm(). function y = normalize(x) L = norm(x); y = x/L; endfunction Variables x and y are vector variables. Please note how vector and scalar operands are handled for division operator on line 3.

en:teaching:subjects:it:labs:sol_3_1_3

Anonymous (anonymous@undisclosed.example.com) — 2017-10-02T15:54:45+0200

author=Roman Putanowicz title=Solution to exercise 3.1.3 back Solution to exercise 3.1.3 There are three main issues in this problem: * element-by-element operations on vectors * controlling plot annotations * saving plot in a file

en:teaching:subjects:it:labs:sol_4_1_1

Anonymous (anonymous@undisclosed.example.com) — 2017-10-02T15:54:45+0200

author=Roman Putanowicz title=Solution to exercise 4.1.1 back Solution to exercise 4.1.1 Firstly we should note that our expression can be transformed to an equivalent form: \begin{equation*} \sum_{j=1}^{j=N}b_j\left(\sum_{i=1}^{i=j}a_i\right) \end{equation*} Then we should also note that there is no need to calculate the inner sum from scratch for each new

en:teaching:subjects:it:labs:sol_4_1_2

Anonymous (anonymous@undisclosed.example.com) — 2017-10-02T15:54:45+0200

author=Roman Putanowicz title=Solution to exercise 4.1.2 back Solution to exercise 4.1.2 The problem can be solved with a simple application of “while” loop: L = input(“Give L : ”); Rd = L; steps = 0; while Rd > 0.1*L;

en:teaching:subjects:it:labs:sol_5_1_1

Anonymous (anonymous@undisclosed.example.com) — 2017-10-02T15:54:45+0200

author=Roman Putanowicz title=Solution to exercise 5.1.1 back Solution to exercise 5.1.1 N = 20; x = floor(N*rand()); printf(“Guess an integer number in range <0,%d>\n”,N); while 1 g = input("Your gueess : "); if x == g disp("Congratulations, this is the right guess."); exit(0); endif

en:teaching:subjects:it:labs:sol_5_1_2

Anonymous (anonymous@undisclosed.example.com) — 2017-10-02T15:54:45+0200

author=Roman Putanowicz title=Solution to exercise 5.1.2 back Solution to exercise 5.1.2 A = input(“Give A : ”); x = input(“Give x0 : ”); epsilon = input(“Give epsilon : ”); n = 0; error = abs(sqrt(A) - x); while ( error >= epsilon)

en:teaching:subjects:it:labs:sol_5_1_3

Anonymous (anonymous@undisclosed.example.com) — 2017-10-02T15:54:45+0200

author=Roman Putanowicz title=Solution to exercise 5.1.3 back Solution to exercise 5.1.3 R = input(“Give R : ”); function p = polyarea(R,n) p = n * 1/2 * R^2 * sin(2*pi/n); endfunction epsilon = 0.005; trueArea = pi*R^2;

en:teaching:subjects:it:labs:sol_6_1_1

Anonymous (anonymous@undisclosed.example.com) — 2017-10-02T15:54:45+0200

author=Roman Putanowicz title=Solution to exercise 6.1.1 back Solution to exercise 6.1.1 h = figure(); axis([-5,5,-5,5]); xy = []; # array of vertex coordinates stp = 0; quitchar = [-1, toascii(“Q”)]; closechar = [toascii(

en:teaching:subjects:it:labs:sol_6_1_2

Anonymous (anonymous@undisclosed.example.com) — 2017-10-02T15:54:45+0200

author=Roman Putanowicz title=Solution to exercise 6.1.2 back Solution to exercise 6.1.2 The $N$ segments that partition the parameter space are defined by values \begin{equation*} \{t_1, t_2,\ldots,t_{N+1}\}\enspace . \end{equation*} For each value of $t_{i}$ we can calculate the corresponding radius vector \begin{equation*} \mathbf{r}_i = \mathbf{r}(t_i) = [x(t_i), y(t_i)] \enspace . \end{equation*}\begin{equation*} \Delta\mathbf{r}_i = \mathbf{r}_{i+1} - \ma…

en:teaching:subjects:it:labs:sol_6_1_3

Anonymous (anonymous@undisclosed.example.com) — 2017-10-02T15:54:45+0200

author=Roman Putanowicz title=Solution to exercise 6.1.3 back Solution to exercise 6.1.3 function rotatecolor(R,a,n) fi = linspace(0,2*pi-2*pi/n,n); x = R*cos(fi); y = R*sin(fi); h = zeros(n); for i=1:n h(i) = square(x(i),y(i),a); endfor drawnow() circulate(h, n);

en:teaching:subjects:it:labs:sol_6_3

Anonymous (anonymous@undisclosed.example.com) — 2017-10-02T15:54:45+0200

function rotatecolor(R,a,n) fi = linspace(0,2*pi-2*pi/n,n); x = R*cos(fi); y = R*sin(fi); h = zeros(n); for i=1:n h(i) = square(x(i),y(i),a); endfor drawnow() circulate(h, n); endfunction function circulate(h, N) n = length(h); i = 1; k = 0; while k < N c = get(h(i), "facecolor"); set(h(i),"facecolor","red"); drawnow() fname = sprintf("frame%03d.png", i); print(fname, "-dpng","-S300,300"); set(h(i),"…

en:teaching:subjects:it:labs:solp1p1

Anonymous (anonymous@undisclosed.example.com) — 2017-10-02T15:54:45+0200