#!/usr/bin/env python3

from scipy.optimize import root

# gegebene Spannungen
U1 = 19
U3 =  9

# R1 kann beliebig vorgegeben werden
R1 =  1

def f(x):
  # gesuchter Verktor x enthält die Elemente U2 und R2
  U2, R2 = x
  
  # Summe der Ströme in Knoten 1
  s1 = (U1 - U2) / R1 + (U3 - U2) / R1 - U2 / R2

  # Summe der Ströme in Knoten 2
  s2 = (U2 - U3) / R1 - U3 / R2

  # Beide Summen sollen am Ende 0 sein
  return [s1, s2]

# Startwerte
U2s = 0
R2s = 1

# Numerische Nullstellensuche
res = root(f, [U2s, R2s])

if res.success:
  print('Lösung:', res.x)
else:
  print('Keine Lösung gefunden')