commit 9622dc7e2fce6ae7349338ae2be289e42d72795e
author avm99963 <me@avm99963.com> Fri Apr 29 15:44:05 2022 +0200
committer avm99963 <me@avm99963.com> Fri Apr 29 15:53:16 2022 +0200
tree ea36e1e1eb2d0678357e0eb140abd9f67dc0e7b8
parent 72184a0f0ebdc29d56961fe8b307bb50afbb66d9

Afegir sessions 1, 2 del lab de fenòmens

A la sessió 1 falta completar l'exercici extra.

Change-Id: I76f9a74f7648b77c3c60ac91cc4a29477386f383

quad10/fenomens/lab/p2/energ.f90

diff --git a/quad10/fenomens/lab/p2/energ.f90 b/quad10/fenomens/lab/p2/energ.f90
new file mode 100644
index 0000000..4ee8508
--- /dev/null
+++ b/quad10/fenomens/lab/p2/energ.f90
@@ -0,0 +1,14 @@
+real*8 function ENERG(S, L, PBC)
+  integer*2 :: S(1:L, 1:L)
+  integer*4 :: I, J, L
+  integer*4 :: PBC(0:L+1)
+  real*8 :: ENE
+  ENE = 0.0d0
+  do I = 1, L
+    do J = 1, L
+      ENE = ENE - S(I, J)*S(PBC(I + 1), J) - S(I, J)*S(I, PBC(J + 1))
+    enddo
+  enddo
+  ENERG = ENE
+  return
+endfunction