Interpoláció és extrapoláció polinomokra
A segédprogramot Farkas Erika írta.
A matematikai háttér:
Az segédprogrammal könnyedén elvégezhetjük az extra- illetve interpolációs eljárást, azaz ismert értékekbõl az ismeretlenekre való következtetést.
A szubrutin használata:
A program meghívása:
CALL POLINT(xa,ya,n,x,y,dy) sor begépelésével történik.
A zárójelben szereplõ értékeket a fõprogramban
is deklarálni kell.
Az n és x változókat a fõprogramban
be kell kérni az ismert pontok koordinátáival együtt.
A bekért pontok x koordinátáit az xa nevû
vektorba, y koordinátáit az ya vektorba kell eltárolni.
Ezután lehet a subroutine-t meghivni. Végül a fõprogramban
y értékét kell kiiratni.
A szubrutin meghívása:
CALL POLINT(xa,ya,n,x,y,dy)
Az egyes változók típusa a következõ:
n: egész szám : az ismert pontok száma
x: valós szám : amelyik pontban keressük az értéket
(a pont x koordinátája)
y: valós szám : a keresett érték (a
pont y koordinátája)
xa, ya: n elemû valós vektorok; a program az xa
vektorban tárolja rendre a megadott pontok x koordinátáját,
ya-ban pedig értelemszerûen az y koordinátákat
dy: hibabecslés, valós szám
Példa a segédprogram használatára:
program teszt
integer n,i
real xa(20),ya(20),x,y,dy
CALL POLINT(xa,ya,n,x,y,dy)
write (*,*)"A keresett ertek: y=",y
stop
end
Próbafuttatás:
Az ismert pontok szama 2
Az 1. pont x koordinátája
5
Az 1. pont ykoordinátája
6
Az 2. pont x koordinátája
6
Az 2. pont y koordinátája
8
Melyik pontban keres az erteketi? x= 99
A keresett érték y=194
Innen töltheted le a segédprogramot: polint.for