ΣΕΜΦΕ forum

Κώδικας: Επιλογή όλων

 program Newton_Raphson
     
      implicit none
      integer P,i
      parameter(P=50)
      real  x(P)
      real x0,e,ei
      real f,fd
       

      print *,"# Enter x0:"
      read(5,*)x0
      print *,"# Enter the error:"
      read(5,*)e
       
      if (e .le. 0.0) stop 'illegal value of error'
      i=2
      x(1)=x0 

      open(unit=20,file="newton_r2.dat")
      write(20,*)1,x(1)

      do while(i .le. P)
      x(i)=x(i-1)-f(x(i-1))/fd(x(i-1))
      ei=abs(x(i)-x(i-1))/x(i)
      print *,i,x(i)
      write(20,*)i,x(i),ei
      if (ei .le. e) stop 'Procedure completed successfully'
      i=i+1
      enddo
      print *,'# Method failed after P=',P,'iterations'
      close(20)
      end
     
      real function f(x)
      implicit none
      real x
      f=sqrt(2.0)*x*(1-x)-x
      return
      end
     
      real function fd(x)
      implicit none
      real x
      fd=sqrt(2.0)+2*sqrt(2.0)*x-1
      return
      end
      

Κώδικας: Επιλογή όλων

      program kinisi
      implicit none
      real*8 t,tf,dt,x0,y0,w1,w2,vx,vy,ax,ay,x,y
c------INPUT
      print *, 'enter x0, y0, w1, w2, tf, dt'
      read(5,*) x0, y0, w1, w2,tf,dt
      print *, 'ok malaka'
c------INITIALIZE
      x=x0
      y=y0
      t=0
      open(unit=11,file='kinisi.dat')
c------COMPUTE
      do while (t .le. tf)
         x = x0*cos(w1*t)
         vx=(x0/w1)*sin(w1*t)
         ax=(-x0/(w1*w1))*cos(w1*t)
         y=y0*sin(w2*t)
         vy=(-y0/w2)*cos(w2*t)
         ay=(-y0/(w2*w2))*sin(w2*t)                        
         t=t+dt
         write(11,*)t,x,y,vx,vy,ax,ay
      enddo
      close(11)
      end