---------------
#  Aufgabe 3  #
---------------


# der Buchstabe "t" ist schon verbraucht, 

t         # ist Funktion fuer Transponieren

# der Buchstabe "f" ist noch frei:

f         # Error: object 'f' not found

# Wir definieren eine Funktion f:

f = function( z )
{
  result = exp(-z^2)
  return(result)
}

f         # ok, ist jetzt eine Funktion
f(1i)     # 2.718282+0i


# Wir wollen die Cauchy'sche Integralformel
# ueberpruefen: die Kurve C ist ein Kreis 
# mit Radius r um z0:
# Die Funktion soll zwei Werte zurueckgeben, 
# das Integral und f(z0):

TestCauchyInt = function( f , z0 , r )
{
  n = 10000
  phi = seq(from=0,to=2*pi,length=n)
  z = z0 + r*exp(1i*phi)
  dz = diff(z)
  z = z[-1]
  integral = sum( f(z)/(z-z0) * dz ) * 1/(2*pi*1i)
  exact = f(z0)
  result = c(integral,exact)
  names(result) = c("Integral","f(z0)")
  return(result)
}

# ein Funktionsaufruf:
TestCauchyInt(f,z0=1i,r=2)