Introduktion till vetenskapliga ber¨ akningar II. Ovning 4, ht 2009 ¨
1. F¨oljande Fortran-program inneh˚aller ett antal fel. F¨ors¨ok hitta dem.
PROGRAM test1
! Detta program inneh˚aller n˚agra fel, f¨ors¨ok finna dem
PRINT *, Skriv ett tal READ * talet
PRINT *, "Talet som du skrev var ", talet
STOP
END PROGRAM test1
F¨ors¨ok kompilera programmet och se om det fungerar efter korrigering!
2. Vilken utskrift f˚ar man fr˚an f¨oljande program? (fundera f¨orst innan du provar!)
PROGRAM test2 IMPLICIT NONE REAL :: a,b,c,d,e INTEGER :: f,g,h a = 1.5
b = 2.
c = a + b f = a + b d = a/b g = a/b e = d*a+b h = g*a+b PRINT *, c,d,e PRINT *, f,g,h STOP
END Program test2
3. Sv¨angningstiden f¨or en matematisk pendel ¨ar T = 2πq
l
g, d¨ar tyngdkraftsaccelera- tionen g = 9.81 m/s2. Skriv ett program, som ber¨aknar sv¨angningstiden f¨or en pendel av given l¨angd.
4. Skriv ett program, som l¨aser in komponenterna av en tredimensionell vektor aaa =aa xııııı + y + zkkk och ber¨kk aknar dess l¨angd |aa| =aaa px2 + y2 + z2.
5. I Fibonacci-serien ¨ar varje tal summan av de tv˚a f¨oreg˚aende talen: 1, 1, 2, 3, 5, 8, . . ..
Skriv ett program som genererar ett givet antal tal i serien, och skriver ut dem.
6. Ber¨akna ett n¨armev¨arde f¨or summan av f¨oljande alternerande serie:
∞
X
k=1
(−1)k−1
k .
Pr¨ova olika summeringsordning, och f¨ors¨ok ocks˚a summera positiva och negativa termer skilt f¨or sig. Summans exakta v¨arde ¨ar ln 2 ≈ 0.69314718 . . ..