Bisekcija - iskanje ničel

 

Naloga:
Napišite program, ki poišče ničle funkcije.

program bisekcija
implicit none
real :: a,b,c,f
write(*,*) 'vnesi levo krajisce intervala na katerem je nicla, a=?'
read(*,*)a
write(*,*) 'vnesi desno krajisce intervala na katerem je nicla, b=?'
read(*,*)b
if (f(a)*f(b)>0) then
   write(*,*) 'funkcija v a in b ni razlicno predznacena, ne bom nadaljeval'
else
 do while (abs(a-b).gt.1e-5)  !1e-5 je natancnost na katero bomo dolocili nicle
  c=(a+b)/2.
  write(*,'(3f10.6,3f10.2,f10.6)')a,b,c,f(a),f(b),f(c),abs(a-b)
  if (f(a)*f(c).lt.0) then
   b=c
  else
   a=c
  endif
 enddo
  write(*,'(3f10.6,3f10.2,f10.6)')a,b,c,f(a),f(b),f(c),abs(a-b)
  write(*,*)'nicla je=',c
endif
end program bisekcija

function f(x)
implicit none
real::x,f
! tu napisemo funkcijo, katere niclo hocemo poiskati
f=0.2*x**5 - 4.02*x**4 - 2.009*x**3 - 4.099*x**2 - 19.0023*x + 50
end function f

Tekst sedaj shranimo kot datoteko bisekcija.f90 in prevedemo.

Tu je sedaj še print screen postopka prevajanja. V kolikor vam kakšna podrobnost ni jasna, mi sporočite.



Nazaj na osnovno stran za numerične metode za kemike.

 

Fakulteta za kemijo in kemijsko tehnologijo

Katedra za fizikalno kemijo

Univerza v Ljubljani

Večna pot 113, K3.009

1000 Ljubljana

Slovenija

tel: +386 1 479 8540

e-pošta: tomaz pika urbic at fkkt.uni-lj.si

Govorilne ure so po dogovoru.

Prešernove nagrade