human
I am working on the file "4_5-Potentialtheorie_Quell-Panel-Verfahren.ipynb".
The first task is:
Unser Quell-Paneel macht schon (fast) was es soll. Es hat nur noch einen Schönheitsfehler: aus dem Paneel scheint entgegen der Strömungsrichtung ein Volumenstrom auszutreten. Wir möchten aber, dass das Paneel zur Anströmrichtung hin eine undurchlässige Wand simuliert. Die stromabgewandte Seite wird später auf der Innenseite des umströmten Körpers liegen, so dass die dort austretende Strömung keine Bedeutung hat).
Wir müssen die Stärke $\lambda$ also genau so wählen, dass der Staupunkt gerade auf dem Paneel zum liegen kommt. Da der pro Länge und Einheitstiefe austretende Volumenstrom ohne überlagerte Translationsströmung zu beiden Seiten des Quell-Panels gleich ist, tritt auf der der Strömung zugewandten Seite genau $\lambda/2$ aus. Damit der Staupunkt also auf dem Quell-Paneel liegt, muss dieser auf die Länge und Einheitstiefe bezogene Volumenstrom also gerade durch die Translationsströmung kompensiert werden:
$$\frac{\lambda}{2} = u_\infty$$
Wenn wir oben im Python-Code lamb = 1 setzen, sollten wir also das gewünschte Ergebnis bekommen. 
Überlagerung mehrerer Paneele zur Beschreibung eines umströmten Körpers
Zur Abbildung des umströmten Körpers mit Paneelen, müssen wir sicherstellen, dass die Geschwindigkeit am Ort der Paneele genau tangential verläuft. Anders ausgedrückt, muss die Geschwindigkeitskomponente normal zum Paneel $u_n$ gleich Null sein. Der Einfachheit halber fordern wir, dass $u_n$ nur im Mittelpunkt des Paneels Null sein muss.
Da die Normalenkomponente der Geschwindigkeit gerade die Ableitung der Potentialfunktion nach der Normalenrichtung ist, gilt:
$$u_n(x_m,y_m) = \frac{\partial \phi}{\partial n}\left(x_m, y_m\right) \stackrel{!}{=} 0$$
Und jetzt wird's etwas komplizierter. Denn die Potentialfunktion $\phi$ in dieser Gleichung ist nicht einfach die eines einzelnen Paneels, sondern die aller $N$ linear überlagerter Paneele und der überlagerten Translationsströmung:
$$\phi(x,y) = u_\infty x + v_\infty y + \sum_{i=1}^N...
Can you write Python program for it?

What is a the Python program you requested?
assistant
#Here is the Python program you requested.

def norm_int(panel_i, panel_j):
    # Berechnet das Integral in der Koeffizientenmatrix
    
    def integrand(s):
        delta_x = panel_j.mx - panel_i.x(s)
        delta_y = panel_j.my - panel_i.y(s)
        return ((delta_x*math.cos(panel_j.beta) 
                + delta_y*math.sin(panel_j.beta))
                / (delta_x*delta_x + delta_y*delta_y))
    
    return integrate.quad(integrand, 0.0, panel_i.laenge)[0]