PROGRAM rc1
      CALL start(t,R,V,C,charge,dt,ncalc)
      CALL Euler(t,R,V,C,charge,dt,ncalc)
      STOP
      END
      SUBROUTINE start(t,R,V,C,charge,dt,ncalc)
*	pocetno vreme
      t=0.0
*	pocetno naelektrisanje (C)
      charge=0.0
*	otpor (Om)
      R=1000.
*	napon (V)
      V0=1.
*     kapacitet (F)
      C=1.0E-6  
*	korak vremena (s)
      dt=0.001
*	trajanje (s)
      tmax=0.02
*	ukupan broj koraka
      ncalc=tmax/dt+0.5
      RETURN
      END
      SUBROUTINE Euler(t,R,V,C,charge,dt,ncalc)
      DO 10  icalc=1,ncalc	         
*	pocetak intervala
      w=2*3.14*100
      V=V0*cos(w*t)
      slope=V/R-charge/(R*C)
      charge=charge+slope*dt	
*	vreme	
      t=t+dt
	Vr=R*slope	
      CALL output(t,charge,Vr,V)	    
   10	CONTINUE
      RETURN
      END
      SUBROUTINE output(t,charge,Vr,V)
*	stampanje rezultata	
	WRITE(9,*)t,charge,Vr,V
	RETURN
      END