Voici un code en c++

Bonjour,
je veux que mon code s'arrête lorsque T[i+1]=T[i] et qu'il me renvoie juste le T[i+1] et là c'est pas ce qu'il fait (quand vous l'exécutez vous allez avoir en dernier des valeurs qui se rapprochent de 0.19 et moi je veux qu'en t[1] j'aurai 0.19 et pas en t[27 svp aidez moi]
voici mon code

#include <iostream>
#include <iomanip>
#include <cmath>

using namespace std;
double fonction (double M, double a, double b, double G)
{
float T[27],t[27];
T[0]=1;
t[0]=0;
int n;
for (int i=0;i<27;i++)
{
int j=0;
T[i+1]=sqrt(M/((a+(b*t[j])/2)+(G)*T[i]));
if (T[i]==T[i+1])
{
T[i+1]=T[i];
n=i+1;
break;
}
t[j+1]=t[j]+T[i];
cout<<"t["<<j+1<<"] ="<<setprecision (5) <<t[j+1]<<" "<<"T["<<i<<"] ="<<setprecision (5) <<T[i]<<endl;
j++;
}

system("pause");
}
int main()
{
double M,a,b,G;
cout<<"donnez M,a,b,G:"<<endl;
cin>>M>>a>>b>>G;
fonction (M,a,b,G);
}


merci.