- double precision function prdscl(A,B,n)\r
- integer i,n\r
- double precision A(n), B(n), S\r
-C\r
- S = 0.0D0\r
- do i=1, n\r
- S = S + A(i) * B(i)\r
- enddo\r
-C\r
- prdscl = S\r
- return\r
- end\r
+ double precision function prdscl(A,B,n)
+ integer i,n
+ double precision A(n), B(n), S
+C
+ S = 0.0D0
+ do i=1, n
+ S = S + A(i) * B(i)
+ enddo
+C
+ prdscl = S
+ return
+ end