(* Doppelspalt Beugungsbild. Es wird die Intensitaet dargestellt. *)

(* Initialisierung *)

Schirmabstand = 10;
Spaltbreite = 1;
k = 20;

ds = Spaltbreite/2;

(* Abstaende vom oberen und vom unteren Spalt *)
Da[y_] := Abs[Schirmabstand + I (y - ds)];
Db[y_] := Abs[Schirmabstand + I (y + ds)];

(* Summe der beiden Kugelwellen *)
Welle[y_]  := Exp[I k Da[y]]/Da[y] + Exp[I k Db[y]]/Db[y]

(* Darstellung der Intensitaet im Bereich -10 < y < 10 *)
Plot[Abs[Welle[y]] ,{y,-10,10}]
