Tak to uz tak nejak funguje. Jen mam poslendi problem vsechny hodnoty ktere iteracni metoda vypise jsou rovny 0. nevite nekde co by to mohlo zpusobit?
kdyztak prikladam skript:
program Bifurkacni_diagram
implicit none
integer :: i
real :: dr, xn, r, x0, r0, xnn, rmax, k
x0 = 0.5
r0 = 0.1
dr = 0.01
xn = x0
r = r0
rmax = 5
OPEN (1, FILE = 'data.dat', ACCESS = 'APPEND',STATUS = 'NEW')
do while ( r < rmax)
do i = 1, 1000
xnn = r*xn*(1-xn)
xn = xnn
if ( i > 800) then
WRITE(1,*) "r = ",r," xnn = ",xnn," iterace =",i
endif
enddo
r = r + dr
enddo
close(1)
end program
dikz moc
Offline