clc; clear;
%Eingabe Kabeldaten
R1_str = 0.0991;  %[Ohm/km]
L1_str = 0.287e-3; %[mH/km]

%Eingabe Leitungslänge
l_K1 = 300; %[m]

%Widerstände berechnen
X_K1 = L1_str*(l_K1/1000)*2*pi*50; %[Ohm]
R_K1 = R1_str*(l_K1/1000);  %[Ohm]
S_bez = 0.01; %MVA
U_bez = 230e-3; %kV
Z_bez = U_bez*U_bez/(S_bez); %[Ohm]

%Modell erzeugen
mpc = erstesModell;

%Modell übernimmt Zahlenwerte aus echten Werten
mpc.branch(1,3) = R_K1 /Z_bez;
mpc.branch(1,4) = X_K1 /Z_bez;
mpc.bus(1,10) = U_bez;
mpc.bus(2,10) = U_bez;
mpc.baseMVA = S_bez;
mpc.bus(2,3) = 3000;

%Modellberechnung starten und auf Struktur Regime ausgeben
regime = runpf(mpc);