The code below is small part of a MATLAB code, translate the code below from MATLAB to Python!

global NB NBB NS NG NLB NTR NTRL NT NSHC NSVS NSHR NSH NREG global VSLACK TOLER PBASE VLMAX VLMIN ITMAX global BIND BSN BNAM PG QG PD QD V DEL QGMAX QGMIN VSH global FB TB NCKT YL ZL BL TAP RAT KV LEN TAPMAX TAPMIN TAPSTP global SNO SUS SUSMAX SUSMIN SUSSTP global YB BP BPP LBP PBP UBP LBPP PBPP UBPP global LAM global PGMAX PGMIN VGMAX VGMIN ap bp cp global ofp ifp % output and input file pointers global JBF JB global DELP DELQ lfep lfeq global MF MT X0 global SLFLOW SBUS SLOSS NSI global LNSI LNVI FBI TBI PAI PBI QAI QBI RL XL ZLM2 MF MT MS YS global LMIS LJAC U Pa Pb Qa Qb

LMIS(0*NT+1:1*NT) = U(FBI).^2 + 2*U(FBI).*(Pa.*RL+Qa.*XL-U(TBI)/2) + (Pa.^2+Qa.^2).*ZLM2; LMIS(1*NT+1:2*NT) = U(TBI).^2 + 2*U(TBI).*(Pb.*RL+Qb.*XL-U(FBI)/2) + (Pb.^2+Qb.^2).*ZLM2;

LMIS(2*NT+1:3*NT) = (U(FBI)+Pa.*RL+Qa.*XL).*tan(DEL(TBI)-DEL(FBI)) - (Pa.*XL-Qa.*RL) ; LMIS(3*NT+1:4*NT) = (U(TBI)+Pb.*RL+Qb.*XL).*tan(DEL(FBI)-DEL(TBI)) - (Pb.*XL-Qb.*RL) ;

LMIS(4*NT+1:4*NT+1*(NB-1)) = PG(LNSI) - PD(LNSI) + MF(LNSI,:)*Pa + MT(LNSI,:)*Pb - real(U(LNSI).*YS(LNSI));

LMIS(4*NT+(NB-1)+1:4*NT+(NB-1)+(NB-NG)) = QG(LNVI) - QD(LNVI) + MF(LNVI,:)*Qa + MT(LNVI,:)*Qb + imag(U(LNVI).*YS(LNVI));

global NB NBB NS NG NLB NTR NTRL NT NSHC NSVS NSHR NSH NREG global VSLACK TOLER PBASE VLMAX VLMIN ITMAX global BIND BSN BNAM PG QG PD OD V DEL QGMAX QGMIN VSH global FB TB NCKT YL ZL BL TAP RAT KV LEN TAPMAX TAPMIN TAPSTP global SNO SUS SUSMAX SUSMIN SUSSTP global YB BP BPP LBP PBP UBP LBPP PBPP UBPP global LAM global PGMAX PGMIN VGMAX VGMIN ap bp cp global ofp ifp $output and input file pointers global JBF JB global DELP DELQ lfep lfeg global MF MT XO global SLFLOW SBUS SLOSS NSI global LNSI LNVI FBI TBI PAI PBI QAI QBI RL XL ZLM2 MF MT MS YS global LMIS LJAC U Pa Pb Qa Qb LMIS (0*NT+1:1*NT) = U (FBI).^2 LMIS (1*NT+1:2*NT) U (TBI).^2 + 2*U (FBI). * (Pa. *RL+Qa. *XL-U (TBI)/2) + (Pa.^2+Qa.^2). *ZLM2; + 2*U (TBI). * (Pb. *RL+Qb. *XL-U (FBI)/2) + (Pb.^2+Qb.^2).*ZLM2; LMIS (2*NT+1:3*NT) = (U (FBI) +Pa. *RL+Qa. *XL). *tan (DEL (T :) -DEL (FBI)) - (Pa. *XL-Qa. *RL); LMIS (3*NT+1:4*NT) = (U (TBI) +Pb. *RL+Qb. *XL). *tan (DEL (FBI) -DEL (TBI)) - (Pb. *XL-Qb. *RL); LMIS (4*NT+1:4*NT+1* (NB-1)) = PG (LNSI) - PD (LNSI) + MF (LNSI, :) *Pa + MT (LNSI, :) *Pb - real (U (LNSI).*YS (LNSI)); LMIS (4*NT+ (NB-1) +1:4*NT+ (NB-1) + (NB-NG)) = QG (LNVI) QD (LNVI) + MF (LNVI, :) *Qa + MT (LNVI, :) *Qb + imag (U (LNVI). *YS (LNVI))