Kundur's 1994 Classic 10-Bar System Doesn't Converge in Matpower

[lmrklz]

I am trying to simulate the classical 10-bus system presented by Kundur (1994) using runpf with the newtonpf solver, but the power flow does not converge within 10 iterations, regardless of the parameter changes I apply. However, I have already simulated this same system in other software tools, where it converged without any issues.
Therefore, I would like to ask if anyone could provide some guidance or suggestions regarding possible reasons for this behavior. I am sharing the code used to model the system below, in case there is any mistake or inconsistency that I may have overlooked.

function mpc = case10kundur
%CASE10KUNDUR Power flow data for custom 10-bus Kundur system (Level 1)
% System fully renumbered to buses 1..10.
% Includes 5 parallel lines between buses 5 and 6.
%
% MATPOWER Case Format : Version 2

%%----- Power Flow Data -----%%
mpc.version = '2';

%% System MVA base
mpc.baseMVA = 100;

%%----- Bus Data -----%%
% bus_i type Pd Qd Gs Bs area Vm Va baseKV zone Vmax Vmin
mpc.bus = [
1 3 0 0 0 0 1 1.00 0 345 1 1.1 0.7; % slack
2 2 0 0 0 0 1 1.00 0 345 1 1.1 0.7; % PV
3 2 0 0 0 0 1 1.00 0 345 1 1.1 0.7; % PV
4 1 0 0 0 0 1 1.00 0 345 1 1.1 0.7;
5 1 0 0 0 0 1 1.00 0 345 1 1.1 0.7;
6 1 0 0 0 7.63 1 1.00 0 345 1 1.1 0.7; % shunt 763 MVAr
7 1 3271 1015 0 6.00 1 1.00 0 345 1 1.1 0.7; % load + shunt 600 MVAr
8 1 0 0 0 17.10 1 1.00 0 345 1 1.1 0.7; % shunt 1710 MVAr
9 1 0 0 0 0 1 1.00 0 345 1 1.1 0.7;
10 1 3384 971 0 0 1 1.00 0 345 1 1.1 0.7; % load
];

%%----- Generator Data -----%%
% bus Pg Qg Qmax Qmin Vg mBase status Pmax Pmin
mpc.gen = [
1 3981 0 9999 -9999 0.9800 100 1 9000 0; % G1
2 1736 0 9999 -9999 0.9646 100 1 9000 0; % G2
3 1154 0 9999 -9999 1.0400 100 1 9000 0; % G3
];

%%----- Branch Data (lines + transformers) -----%%
% fbus tbus r x b rateA rateB rateC ratio angle status angmin angmax

mpc.branch = [

%%==================== LINES ====================%%

%% Line 4–5
4   5   0.0000   0.0040   0        250   250   250   0      0   1   -360   360;

%% 5 parallel lines between 5–6
5   6   0.0015   0.0288   0.2346   250   250   250   0      0   1   -360   360;
5   6   0.0015   0.0288   0.2346   250   250   250   0      0   1   -360   360;
5   6   0.0015   0.0288   0.2346   250   250   250   0      0   1   -360   360;
5   6   0.0015   0.0288   0.2346   250   250   250   0      0   1   -360   360;
5   6   0.0015   0.0288   0.2346   250   250   250   0      0   1   -360   360;
%% Line 8–9
8   9   0.0010   0.0030   0        250   250   250   0      0   1   -360   360;

%%================ TRANSFORMERS ================%%

%% T1: 1–4
1   4   0        0.0020   0        250   250   250   0.8857  0   1   -360   360;

%% T2: 2–5
2   5   0        0.0045   0        250   250   250   0.8857  0   1   -360   360;

%% T3: 3–6
3   6   0        0.0125   0        250   250   250   0.9024  0   1   -360   360;

%% T4: 6–7
6   7   0        0.0030   0        250   250   250   1.0664  0   1   -360   360;

%% T5: 6–8
6   8   0        0.0026   0        250   250   250   1.0800  0   1   -360   360;

%% T6: 9–10 (tap variable – level 1)
9  10   0        0.0010   0        250   250   250   0.9750  0   1   -360   360;

];