DIMEN3

 !Program name: DIMEN3

real, dimension(10) :: A=(/ 1,2,3,4,5,6,7,8,9,10 /)

Real, dimension(10) :: B

where (A>5)

    B=1.

Elsewhere

    B=0

end where

write (*,*) A, B

stop

end

Comentários

Postar um comentário

Postagens mais visitadas deste blog

RECURSIVE

  !Program name: RECURSIVE !calculate the factorial of a number n! program recursive      integer :: i, f    character (len=1) :: answer    Print *, 'Input the value to be factorized: '    Read *, i ! Variable function:    f = factorial(i)    Print *, 'The value of factorial of', i, ' is: '    Print *, f 10   Print *, 'Do you want to calculate another factorial (Y/N)?'    Read *, answer    If (answer == 'Y') then         Print *, 'Input another value to be factorized: '         Read *, i         f = factorial(i)         Print *, 'The value of the factorial of ', i, ' is: '         Print *, f         Go to 10       else       print *, 'Ok. Goodbye.'    end if   end program recursive   ! computes the factorial of n (n!) recursive function factorial (n) result (fac) ! function result implicit none      ! dummy arguments    integer :: fac    integer, i
Leia mais

HYDROGEN-RADIAL

!--------------------------------------------------------------- !Program name: HYDROGEN-RADIAL !Radial equation of the hydrogen-like atom !Objective: plot the radial wave function for each n,l pair ! The equation is atomic units !--------------------------------------------------------------- !The Coulomb potential energy equation is: ! V(r) = -2Z/r ! The ceentrifugal potential equation is: ! VL=l(l+1)/r^2 ! The energy is: ! E=-2*Z/n^2 !--------------------------------------------------------------- ! Use of Numerov method to find the wave function ! Function fn from Numerov equation is ! Fn = 1 + K2*h^2/12 , where K2=-2(Veff - E) ! In the code, h is changed into dr !---------------------------------------------------------------   Program hydrogen   Implicit none   integer :: mesh, n, l, i, j, zeta, m, maxiter=100   Real*8 rmin, dr, zmesh, y_out_m   double precision :: e, fh12, norm, eps=1.0E-6    double precision, allocatable :: r(:), r2
Leia mais

RUNGE-KUTTA

  !Program name: RK !Runge-Kutta method for differential equation dy/dx=-xy !The exact solution: y=exp(-x^2/2) !Interval X=0 to X=3 where y(0)=1 !Y(n+1) <=> X+h   program Runge_Kutta real :: x, y, h, exact, error Real :: k1, k2, k3, k4 integer :: i f(x,y)=-x*y   !fourth order Runge-Kutta h=0.05 nstep=3./h  !60 steps y=1.   !y(0)=1 do i=0,nstep    x=h*i    k1=h*f(x,y)    k2=h*f(x+0.5*h,y+0.5*k1)    k3=h*f(x+0.5*h,y+0.5*k2)    k4=h*f(x+h,y+k3)    y=y+(1/6.*(k1+2.*k2+2.*k3+k4))    exact=exp((-(x+h)**2)/2)    error=exact-y    print *, i, x+h, y,exact, error !for h=0.05  y(1)=0.60653    error=-0.00000 !for h=0.05  y(3)=0.00111    error=-7.63E-08 end do Print *,   !Second order Runge Kutta h=0.05 nstep=3./h  !60 steps y=1.   !y(0)=1 do i=0,nstep    x=i*h    !It guarantees interval from 0 to 3    k=h*f(x,y)    y=y+h*f(x+0.5*h,y+0.5*k)    exact=exp((-(x+h)**2)/2)    error=exact-y    print *, i, x+h, y,exact
Leia mais