
Re: Handle function implicitly accounting for independent variables
Posted:
Jan 7, 2013 2:08 PM


Dear Steven Lord,
I changed my approach a little bit.
At the moment I'm doing like this: 1. First I wrote a standard function phi_33_kai
function phi_33 = phi_33_kai(k1,k2,k3,beta,i) k = sqrt(k1(i).^2 + k2.^2 + k3.^2); k30 = k3 + k1(i).*beta(i); k0 = sqrt(k1(i).^2 + k2.^2 + k30.^2); E_k0 = 1.453.*(k0.^4./((1 + k0.^2).^(17/6))); phi_33 = E_k0./(4.*pi.*k.^4).*(k1(i).^2 + k2.^2); end
2. Then, I call that in the previous for loop as follows:
clc,clear all,close all k1 = (1E02:0.1:100); load beta.mat %% the procedure to calculate beta is not reported anymore (see first msg)
F_33 = zeros(1,numel(k1)); count = 0; for i = 1:numel(k1) count = count + 1; phi_33 = @(k2,k3) phi_33_kai(k1,k2,k3,beta,i); F_33(count) = quad2d(phi_33,100,100,100,100); end
Is this way comparable to yours or it has any relevant flaw?
Btw, would you mind to extend your approach with a whole piece of code: I was not able to completely understand it. I thank you in advance.
Kindest regards, Francesco Perrone

