diff --git a/CMakeLists.txt b/CMakeLists.txt
index 96047eefe..386febb30 100644
--- a/CMakeLists.txt
+++ b/CMakeLists.txt
@@ -6,7 +6,7 @@
 cmake_minimum_required(VERSION 3.23)
 project(
   BSMPT
-  VERSION 3.0.2
+  VERSION 3.0.3
   LANGUAGES C CXX
   DESCRIPTION
     "BSMPT - Beyond the Standard Model Phase Transitions : A C++ package for the computation of the EWPT in BSM models"
diff --git a/README.md b/README.md
index be9a69e7d..9881d473f 100644
--- a/README.md
+++ b/README.md
@@ -4,7 +4,7 @@ SPDX-FileCopyrightText: 2021 Philipp Basler, Margarete Mühlleitner and Jonas M
 SPDX-License-Identifier: GPL-3.0-or-later
 -->
 
-Program: BSMPT version 3.0.2
+Program: BSMPT version 3.0.3
 
 Released by: Philipp Basler, Lisa Biermann, Margarete Mühlleitner, Jonas Müller, Rui Santos and João Viana
 
diff --git a/tools/ModelGeneration/sympy/SM.py b/tools/ModelGeneration/sympy/SM.py
index f6e47c52b..b6e305003 100644
--- a/tools/ModelGeneration/sympy/SM.py
+++ b/tools/ModelGeneration/sympy/SM.py
@@ -13,7 +13,7 @@
 sigma0 = Matrix([[1,0],[0,1]])
 sigma1 = Matrix([[0,1],[1,0]])
 sigma2 = Matrix([[0,-I],[I,0]])
-sigma3 = Matrix([[1,0],[0,-1]]) 
+sigma3 = Matrix([[1,0],[0,-1]])
 m_electron = symbols('C_MassElectron',real=True)
 m_mu = symbols('C_MassMu',real=True)
 m_tau = symbols('C_MassTau',real=True)
@@ -69,7 +69,7 @@
 ER = symbols('eR muR tauR', real=True)
 EL = symbols('eL muL tauL', real=True)
 
-ye = sqrt(2)*m_electron/v 
+ye = sqrt(2)*m_electron/v
 ymu = sqrt(2)*m_mu/v
 ytau = sqrt(2)*m_tau/v
 
@@ -94,9 +94,9 @@
 DR = symbols('dR sR bR', real=True)
 QuarkBase = UL + DL + UR + DR
 
-yb = sqrt(2)*m_bottom/v 
-yc = sqrt(2)*m_charm/v 
-yd = sqrt(2)*m_down/v 
+yb = sqrt(2)*m_bottom/v
+yc = sqrt(2)*m_charm/v
+yd = sqrt(2)*m_down/v
 ys = sqrt(2)*m_strange/v
 yt = sqrt(2)*m_top/v
 yu = sqrt(2)*m_up/v
@@ -109,10 +109,10 @@
 URVector = Matrix([[x] for x in UR])
 DRVector = Matrix([[x] for x in DR])
 
-VQuark = ULVector.transpose() * VCKM * DownCoupling * DRVector * phi[0] 
-VQuark+= DLVector.transpose() * DownCoupling * DRVector * phi[0] 
+VQuark = ULVector.transpose() * VCKM * DownCoupling * DRVector * phi[0]
+VQuark+= DLVector.transpose() * DownCoupling * DRVector * phi[1]
 VQuark+= ULVector.transpose() * UpCoupling * URVector * phi[1].conjugate()
-VQuark+= -DLVector.transpose() * Dagger(VCKM)*UpCoupling*URVector*phi[1].conjugate()
+VQuark+= -DLVector.transpose() * Dagger(VCKM)*UpCoupling*URVector*phi[0].conjugate()
 VQuark = simplify(VQuark[0,0])
 
 
@@ -152,6 +152,6 @@
 
     if method == 'treeSimpl':
         toyModel.printTreeSimplified()
-    
+
     if method == 'CTSimpl':
-        toyModel.printVCTSimplified()
\ No newline at end of file
+        toyModel.printVCTSimplified()