function [lamm_pip,A0disk, lammbase,A0proj] = Make_lam_vect(lampole, lamtip, INPUTS, j) % This function takes the starting values contained in the vector lampole % and integrates ove a biconcave disk contour to obtain the corresponding % values for lamda m and Ao at the base of the pipet. vecsize = length(lampole); % Initialize arrays lamm_pip = zeros(1,vecsize); A0disk = zeros(1,vecsize); % Taum_tip = zeros(1,npts); lammbase = zeros(1,vecsize); A0proj = zeros(1,vecsize); % Taum_proj = zeros(1,npts); dims = INPUTS{3}; Rp = dims{3}; for i = 1:vecsize % First we calculate the integration end values for the disk % disp(num2str(lampole(i))); [res, ~] = do_outer_disk_from_pole_fitd(INPUTS,j,lampole(i)); r0tip = get_r0_disk(res(2,end), INPUTS); lamf_pip = Rp/r0tip; lamm_pip(i) = get_lamm_of_Taum_fitd(res(3,end),lamf_pip,INPUTS); A0disk(i) = res(2,end); % Taum_tip(i) = res(3,end); % Then we calculate the end values for the projection. [resproj, ~] = do_projection_disk_fitd(INPUTS,j,lamtip(i)); A0proj(i) = resproj(2,end); % Taum_proj(i) = resproj(3,end); r0base = get_r0_disk(A0proj(i),INPUTS); lamfbase = Rp/r0base; lammbase(i) = get_lamm_of_Taum_fitd(resproj(3,end),lamfbase,INPUTS); end