Changes between Version 3 and Version 4 of fortran

Timestamp:

Oct 17, 2016, 8:59:10 PM (8 years ago)

Author:

Leon Kos

Comment:

Reposition and page outline

fortran

@@ +1 @@
+[[PageOutline]]
 = FORTRAN 95 examples =
 
@@ -33 +34 @@
 end
 }}}
+
+== Console2/Console2/Console2.f90 ==
+{{{
+#!fortran
+real x,dx,s
+!open(1,file='y1.dat', status='unknown')
+a=0
+b=1.57
+n=100
+dx=(b-a)/n
+S=0
+do i=0,n
+    if(mod(i,2)==0) then 
+       S=S+dx/3*2*f(i*dx+a)
+   else 
+       S=S+dx/3*4*f(i*dx+a)
+   end if
+end do
+ 
+print *, 'Sim = ', S 
+pause
+
+contains
+
+function f(x) 
+f=sin(x)*sin(2*x)*sin(3*x)
+end function f
+end
+}}}
+
+== Console3/Console3/Console3.f90 ==
+{{{
+#!fortran
+real x,dx, a, b
+open(1,file='y2.dat', status='unknown')
+a=-0.5
+b=0.5
+n=100
+dx=(b-a)/n
+do i=1,n
+    write(1,*) i*dx + a, 1/sqrt(1-(i*dx+a)**2)
+end do
+end
+}}}
+
+== Console4/Console4/Console4.f90 ==
+{{{
+#!fortran
+real x,dx, a, b
+open(1,file='y2.dat', status='unknown')
+a=0
+b=1.14
+n=100
+dx=(b-a)/n
+do i=1,n
+    write(1,*) i*dx + a, sin(i*dx+a)*sin(2*i*dx+a)*sin(3*i*dx+a)
+end do
+end
+}}}
+
+== Console5/Console5/Console5.f90 ==
+{{{
+#!fortran
+program Simp
+
+integer n
+real y1, y2, y3, x, h, a, s
+n=100
+h=2/n
+do i=1,n
+   
+        x(i)=-1+h(i-1)
+
+        y1(i)=f1(x(i))
+        y2(i)=f2(x(i),h)
+        y3(i)=f3(x(i),h)
+        
+         a=h/6*(y1+4*y2+y3)
+         s=sum(a)
+end do
+
+contains
+real function f1(x)
+    real x
+    f1=x*(5-4*x)**(-0.5)
+  end function
+real function f2(x,h)
+    real x, h
+    f2=(x+h/2)*(5-4*(x+h/2))**(-0.5)
+  end function
+real function f3(x,h)
+    real x, h
+    f3=(x+h)*(5-4*(x+h))**(-0.5)
+
+ end function
+print *, s
+end program
+}}}
+
+== Console6/Console6/Console6.f90 ==
+{{{
+#!fortran
+program kolco
+
+integer :: n=1000
+real :: r1=1.0, r2=2.0, q=3.0, d=5.0, eq=-4.8e-10, m=9.1e-28 ! ������� ��������� � ���-�������
+real ::  dt, t, x, vx
+
+t=1.0e-9
+dt=t/float(n)
+vx=0
+x=2.5
+
+
+open (1,file='kolco',status='unknown')
+
+do i=1,n 
+   E=q*x/(x**2+r1**2)**(1.5)- q*(5-x)/((5-x)**2+r2**2)**(1.5)
+   vx=vx+(eq*E/m)*dt
+   x=x+vx*dt
+   write(1,*) vx, x, i*dt
+end do
+
+end program
+}}}
+
+== Console7/Console7/Console7.f90 ==
+{{{
+#!fortran
+program kolco
+
+integer :: n=10000
+real(8) :: r1=1.0, r2=2.0, q=3.0, d=5.0, eq=4.8e-10, m=9.1e-28 ! ������� ��������� � ���-�������
+real(8) ::  dt, t, x, vx
+
+t=1.0e-6
+dt=t/float(n)
+vx=0
+x=2.5
+
+
+open (1,file='kolco',status='unknown')
+
+do i=1,n 
+   E=q*x/sqrt((x**2+r1**2)**3)- q*(5-x)/sqrt(((5-x)**2+r2**2)**3)
+   vx=vx+(eq*E/m)*dt
+   x=x+vx*dt
+   write(1,*) x, i*dt, vx
+end do
+
+end program
+}}}
+
+== Console8/Console8/Console8.f90 ==
+{{{
+#!fortran
+external fct, out
+real aux(8,2)
+common /a/ dt, k, n, A, Vm, om
+real :: pr(5) = (/0.,2*3.14159*10., 2*3.14159/50, .0001, .0/)
+real :: du(2) = (/0.5, 0.5/)
+integer :: nd = 2
+real :: dt=2*3.14159265
+real :: t
+A = 2.0
+om = 1.0
+Vm = A*om
+u(1) = Vm/Vm
+u(2) = 0
+k = 1
+n=1
+
+open (unit=2, file = 'osc_rk0.dat' ,status='unknown')
+write(*,*) 'x=' , u(1), 'v=', u(2)
+pause
+call rkgs(pr, u, du, nd, ih, fct, out, aux)
+write(2,*) 'ih=', ih
+stop
+end
+!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
+
+subroutine fct (t, u, du)
+real u(2), du(2)
+common /a/ dt, k, n, A, Vm, om
+du(1) = -u(2)/(1+0.1*sin(om*dt))
+du(2) = u(1)
+return
+end
+
+
+subroutine out (t, u, du, ih, nd, pr)
+real u(2), du(2), pr(5)
+common /a/ dt, k, n, A, Vm, om
+
+if(t.ge.k*(dt/50.)) then
+write (2,1) t/om, u(1)*Vm, u(2)*A
+k=k+1
+else
+end if
+
+1 format (f8.3, f8.3, f8.3)
+return
+end
+
+
+
+}}}
+
+== Console9/Console9/Console9.f90 ==
+{{{
+#!fortran
+!  Console9.f90 
+!
+!  FUNCTIONS:
+!  Console9 - Entry point of console application.
+!
+
+!****************************************************************************
+!
+!  PROGRAM: Console9
+!
+!  PURPOSE:  Entry point for the console application.
+!
+!****************************************************************************
+
+    program Console9
+
+    implicit none
+
+    ! Variables
+
+    ! Body of Console9
+    print *, 'Hello World'
+
+    end program Console9
+
+}}}
+
 
 == Console10/Console10/Console10.f90 ==
@@ -919 +1158 @@
 }}}
 
-== Console2/Console2/Console2.f90 ==
-{{{
-#!fortran
-real x,dx,s
-!open(1,file='y1.dat', status='unknown')
-a=0
-b=1.57
-n=100
-dx=(b-a)/n
-S=0
-do i=0,n
-    if(mod(i,2)==0) then 
-       S=S+dx/3*2*f(i*dx+a)
-   else 
-       S=S+dx/3*4*f(i*dx+a)
-   end if
-end do
- 
-print *, 'Sim = ', S 
-pause
-
-contains
-
-function f(x) 
-f=sin(x)*sin(2*x)*sin(3*x)
-end function f
-end
-}}}
 
 == Console20/Console20/Console20.f90 ==
@@ -1731 +1942 @@
 }}}
 
-== Console3/Console3/Console3.f90 ==
-{{{
-#!fortran
-real x,dx, a, b
-open(1,file='y2.dat', status='unknown')
-a=-0.5
-b=0.5
-n=100
-dx=(b-a)/n
-do i=1,n
-    write(1,*) i*dx + a, 1/sqrt(1-(i*dx+a)**2)
-end do
-end
-}}}
 
 == Console30/Console30/Console30.f90 ==
@@ -2553 +2750 @@
 }}}
 
-== Console4/Console4/Console4.f90 ==
-{{{
-#!fortran
-real x,dx, a, b
-open(1,file='y2.dat', status='unknown')
-a=0
-b=1.14
-n=100
-dx=(b-a)/n
-do i=1,n
-    write(1,*) i*dx + a, sin(i*dx+a)*sin(2*i*dx+a)*sin(3*i*dx+a)
-end do
-end
-}}}
-
-== Console5/Console5/Console5.f90 ==
-{{{
-#!fortran
-program Simp
-
-integer n
-real y1, y2, y3, x, h, a, s
-n=100
-h=2/n
-do i=1,n
-   
-        x(i)=-1+h(i-1)
-
-        y1(i)=f1(x(i))
-        y2(i)=f2(x(i),h)
-        y3(i)=f3(x(i),h)
-        
-         a=h/6*(y1+4*y2+y3)
-         s=sum(a)
-end do
-
-contains
-real function f1(x)
-    real x
-    f1=x*(5-4*x)**(-0.5)
-  end function
-real function f2(x,h)
-    real x, h
-    f2=(x+h/2)*(5-4*(x+h/2))**(-0.5)
-  end function
-real function f3(x,h)
-    real x, h
-    f3=(x+h)*(5-4*(x+h))**(-0.5)
-
- end function
-print *, s
-end program
-}}}
-
-== Console6/Console6/Console6.f90 ==
-{{{
-#!fortran
-program kolco
-
-integer :: n=1000
-real :: r1=1.0, r2=2.0, q=3.0, d=5.0, eq=-4.8e-10, m=9.1e-28 ! ������� ��������� � ���-�������
-real ::  dt, t, x, vx
-
-t=1.0e-9
-dt=t/float(n)
-vx=0
-x=2.5
-
-
-open (1,file='kolco',status='unknown')
-
-do i=1,n 
-   E=q*x/(x**2+r1**2)**(1.5)- q*(5-x)/((5-x)**2+r2**2)**(1.5)
-   vx=vx+(eq*E/m)*dt
-   x=x+vx*dt
-   write(1,*) vx, x, i*dt
-end do
-
-end program
-}}}
-
-== Console7/Console7/Console7.f90 ==
-{{{
-#!fortran
-program kolco
-
-integer :: n=10000
-real(8) :: r1=1.0, r2=2.0, q=3.0, d=5.0, eq=4.8e-10, m=9.1e-28 ! ������� ��������� � ���-�������
-real(8) ::  dt, t, x, vx
-
-t=1.0e-6
-dt=t/float(n)
-vx=0
-x=2.5
-
-
-open (1,file='kolco',status='unknown')
-
-do i=1,n 
-   E=q*x/sqrt((x**2+r1**2)**3)- q*(5-x)/sqrt(((5-x)**2+r2**2)**3)
-   vx=vx+(eq*E/m)*dt
-   x=x+vx*dt
-   write(1,*) x, i*dt, vx
-end do
-
-end program
-}}}
-
-== Console8/Console8/Console8.f90 ==
-{{{
-#!fortran
-external fct, out
-real aux(8,2)
-common /a/ dt, k, n, A, Vm, om
-real :: pr(5) = (/0.,2*3.14159*10., 2*3.14159/50, .0001, .0/)
-real :: du(2) = (/0.5, 0.5/)
-integer :: nd = 2
-real :: dt=2*3.14159265
-real :: t
-A = 2.0
-om = 1.0
-Vm = A*om
-u(1) = Vm/Vm
-u(2) = 0
-k = 1
-n=1
-
-open (unit=2, file = 'osc_rk0.dat' ,status='unknown')
-write(*,*) 'x=' , u(1), 'v=', u(2)
-pause
-call rkgs(pr, u, du, nd, ih, fct, out, aux)
-write(2,*) 'ih=', ih
-stop
-end
-!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
-
-subroutine fct (t, u, du)
-real u(2), du(2)
-common /a/ dt, k, n, A, Vm, om
-du(1) = -u(2)/(1+0.1*sin(om*dt))
-du(2) = u(1)
-return
-end
-
-
-subroutine out (t, u, du, ih, nd, pr)
-real u(2), du(2), pr(5)
-common /a/ dt, k, n, A, Vm, om
-
-if(t.ge.k*(dt/50.)) then
-write (2,1) t/om, u(1)*Vm, u(2)*A
-k=k+1
-else
-end if
-
-1 format (f8.3, f8.3, f8.3)
-return
-end
-
-
-
-}}}
-
-== Console9/Console9/Console9.f90 ==
-{{{
-#!fortran
-!  Console9.f90 
-!
-!  FUNCTIONS:
-!  Console9 - Entry point of console application.
-!
-
-!****************************************************************************
-!
-!  PROGRAM: Console9
-!
-!  PURPOSE:  Entry point for the console application.
-!
-!****************************************************************************
-
-    program Console9
-
-    implicit none
-
-    ! Variables
-
-    ! Body of Console9
-    print *, 'Hello World'
-
-    end program Console9
-
-}}}
-
-== Eyler/Eyler/Eyler.f90 ==
-{{{
-#!fortran
-real v,x,dt,t,dv,Fm
-open(2,file='y2.dat', status='unknown')
-q1=1.0e-9
-q2=1.0e-9
-r1=0.1
-r2=0.2
-d=0.5
-t0=d/v0
-n=100
-pi=3.141593
-E0=0.0
-v0=0.0
-dt=t0/n
-x=x/d
-write(*,*) dt
-open(1,file='xv.dat', status='unknown')
-Fm=q1/me/(vo**2)
-print*, Fm
-do k=1,n
-    t=k*dt
-    v=v+Fm*sin(2.0*pi*x)*dt
-    x=x+v*dt
-    
-write(1,*) t*t0/1.0e-9, v*v0, x*d, sin(2.0*pi*x)
-if(x.ge.1.0) then
-write(1,*) t*t0/1.0e-9, v*vo, x*d
-stop
-else
-end if
-end do
-end
-}}}
 
 == uskarenie-protona/Console36/Console36.f90 ==
@@ -3056 +3025 @@
 
 }}}
+
+== Eyler/Eyler/Eyler.f90 ==
+{{{
+#!fortran
+real v,x,dt,t,dv,Fm
+open(2,file='y2.dat', status='unknown')
+q1=1.0e-9
+q2=1.0e-9
+r1=0.1
+r2=0.2
+d=0.5
+t0=d/v0
+n=100
+pi=3.141593
+E0=0.0
+v0=0.0
+dt=t0/n
+x=x/d
+write(*,*) dt
+open(1,file='xv.dat', status='unknown')
+Fm=q1/me/(vo**2)
+print*, Fm
+do k=1,n
+    t=k*dt
+    v=v+Fm*sin(2.0*pi*x)*dt
+    x=x+v*dt
+    
+write(1,*) t*t0/1.0e-9, v*v0, x*d, sin(2.0*pi*x)
+if(x.ge.1.0) then
+write(1,*) t*t0/1.0e-9, v*vo, x*d
+stop
+else
+end if
+end do
+end
+}}}
+
+