Diffusion Problems code in Fortran

! 0<=x<=1, 0<=y<=1 ! with boundary conditions phi=xy !* implicit none integer, parameter :: nx=25, ny=25 integer :: i,j,iterations real :: x,y,z,dx,dy,omega real, dimension(1:nx,1:ny) :: AE,AW,AN,AS,AP,SP,SU real, dimension(0:nx+1,0:ny+1) :: phi !*dx=1./float(nx); dy=1./float(ny) ! initialize SU and SP to interior cell values of zeroSU=0; SP=0phi=0.0 ! initialize AE,AW,AN,AS,AP over all cells to! the values over the …

Diffusion Problems code in Fortran Read More »