Υπολογισμός των επι μέρους όρων.Το υποπρόγραμμα: SUBROUTINE SUBSUM

subroutine subsum ( m, n, p, f )

implicit double precision (a-h,o-z)

save

parameter(mx=8000)

dimension p(n), f(m)

common/xydata/x(mx),y(mx),mm

data itric /0/

if(itric.eq.0) then

itric = 1

if(m.gt.mx) stop

mm=m

if( (n/3)*3 .ne. n) stop

nodes = n/3

open(1,file='xy.data',status='old')

do 10 i=1,m

read(1,*)x(i),y(i)

10 continue

close(1)

endif

do 1 i = 1,m

call ann(p,nodes,x(i),per)

f(i)= per-y(i)

1 continue

end

Τα σημεία x-y είναι στο αρχείο: xy.data

Η διαδικασία αυτή εκτελείται μόνο την πρώτη φορά που καλείται η ρουτίνα SUBSUM

Προηγούμενη διαφάνεια

Επόμενη διαφάνεια

Μετάβαση στην πρώτη διαφάνεια

Προβολή έκδοσης γραφικών