public double CalculateReturns(int age, int retireAge, int deathAge, double mean, double stdDeviation, double amount, TaxStatus taxType, FilingStatus status, double income, List <Riders> Riders)
{
double amountWithFees = amount;
if (Riders.Contains(Models.Riders.DeathBenefit))//death benefit rider
{
amountWithFees -= amountWithFees * .005;
}
double[] account = new double[deathAge - retireAge];
int count = 0;
double temp = amountWithFees;
double principle = amountWithFees;
for (int j = age; j < deathAge; j++)
{
double rate = mean;
if (j == retireAge)
{
double assetAtRetire = temp;
}
if (j < retireAge)
{
temp = temp * Math.Pow(1 + rate, 1);
}
if (j >= retireAge)
{
account[count] = CalcWithdrawal(mean, temp, deathAge - j + 1, taxType, status, principle / (deathAge - retireAge));
double Taxamount = TaxHelper.CalcTaxes(TaxHelper.CalcWithdrawalAmount(rate, temp, deathAge - retireAge), taxType, status, principle);
temp -= CalcWithdrawal(mean, temp, deathAge - j + 1, taxType, status, principle / (deathAge - retireAge));
temp = temp * Math.Pow(1 + rate, 1);
count++;
}
}
double averageWithdrawal = 0;
for (int i = 0; i < deathAge - retireAge; i++)//calculates average withdrawal during retirement
{
averageWithdrawal += account[i];
}
averageWithdrawal = averageWithdrawal / (deathAge - retireAge);
return(averageWithdrawal);
}
public double CalcWithdrawal(double rate, double presentValue, int yearsWithdrawing, TaxStatus taxType, FilingStatus status, double principle)
{
return(TaxHelper.CalcTaxedWithdrawals(rate, presentValue, yearsWithdrawing, taxType, status, principle));
}
public double[] CalculateReturns(int age, int retireAge, int deathAge, double mean, double stdDeviation, double amount, TaxStatus taxType, FilingStatus status, double income, List <Riders> Riders)
{//same as deferred variable except lumpsum instead of continuous payments
double amountWithFees = amount;
Boolean isGMWB;
if (Riders.Contains(Models.Riders.GMWB))//Checks to see if GMWB is a rider
{
isGMWB = true;
amountWithFees -= amountWithFees * .005;
}
else
{
isGMWB = false;
}
Boolean isGMAB;
if (Riders.Contains(Models.Riders.GMAB))//checks to see if GMAB is a rider
{
isGMAB = true;
amountWithFees -= amountWithFees * .005;
}
else
{
isGMAB = false;
}
if (Riders.Contains(Models.Riders.DeathBenefit))
{
amountWithFees -= amountWithFees * .005;
}
double[] MedianAverageWithdrawal = new double[4000];
for (int i = 0; i < 4000; i++)
{
// double[] account = new double[deathAge-retireAge];
int count = 0;
double temp = amountWithFees;
double principle = amountWithFees;
double minWithdrawal = 0;
double withdrawalAmount = 0;
for (int j = age; j < deathAge; j++)
{
Random rand = new Random();
double rate = mean + stdDeviation * (rand.NextDouble() * (6) - 3);
if (j == retireAge)
{
if (isGMAB && (principle > temp))
{
temp = principle;
}
minWithdrawal = principle / (deathAge - retireAge);// CalcWithdrawal(mean, principle, deathAge - j + 1, taxType, status, principle);
}
if (j < retireAge)
{
temp = temp * Math.Pow(1 + rate, 1);
}
else
{
double withdrawal = CalcWithdrawal(mean, temp, deathAge - j + 1, taxType, status, principle / (deathAge - retireAge));
double taxes = TaxHelper.CalcTaxes(TaxHelper.CalcWithdrawalAmount(rate, temp, deathAge - retireAge + 1), taxType, status, principle);
temp -= withdrawal;
withdrawal = withdrawal - withdrawal * .03;
if (isGMWB)
{
if (withdrawal < minWithdrawal)
{
withdrawal = minWithdrawal;
}
}
withdrawalAmount += withdrawal;
temp = temp * Math.Pow(1 + rate, 1);
count++;
}
}
withdrawalAmount = withdrawalAmount / (deathAge - retireAge); //calculates average withdrawal
MedianAverageWithdrawal[i] = withdrawalAmount; //stores the average withdrawal for this trial
}
// trials.Add(MedianAverageWithdrawal);//adds an array of the averages to the end of the lsit, will be taken out later and used
return(MedianAverageWithdrawal);
}