So I have a system of equations that I need to iterate until convergence.
I've tried doing this with for loops, while loops, symbolicly, and lastly using approximation methods.
Maybe convergence doesnt mean what I think it means, maybe I'm just missing something.
Here is my code...
clear all; close all; clc;
T = 20000;
rho = .002377;
A = 2827;
v_tip = 650;
R = 30;
a = 6;
Theta_Tw = -10*pi/180;
gamma = 8.1;
Cd0 = 0.01;
Sigma = 0.085;
k = 1.15;
mu = linspace(0.05,0.3,26);
i_hp = [-10; 0; 10]*(pi/180);
phi_hp = 0;
beta1C = 0;
beta1S = 0;
Omega = v_tip/R;
CT = T/(rho*A*v_tip^2);
delta = 1-mu.^2+9*mu.^4/4;
lambda_i_initial = [sqrt(CT/2), CT./(2*mu(2:end))];
lambda = mu.*tan(i_hp)+lambda_i(1,:);
lambda_i = CT./(2*(mu.^2+lambda.^2));
I need to take the initial solution for lambda_i and insert that into the equation for lambda, and then take the lambda solution and insert it into lambda_i, then lambda_i into lambda and so on until convergence.
Any help would be appreciated.
