commit 45a8a462b0f64b2a96796a25ae064ffe252e225d
author avm99963 <me@avm99963.com> Sat Jun 04 12:41:03 2022 +0200
committer Adrià Vilanova Martínez <me@avm99963.com> Sat Jun 04 18:47:37 2022 +0200
tree 6631fcdda187b70505decaceefaeb7934d217c95
parent 802583e7b41e076c94a8279782af52dbd5df2f55

Afegir sessió 4 del lab de fenomens

Change-Id: I68a9c725acd6cde3fd43413fe673e5ec7794f7e6

quad10/fenomens/lab/p4/energ.f90

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