Hi everyone.
I am writing again because my code is not working properly. I need your help
because I can not find the error.
I have to solve the set of equations presented in the pic for a certain time
interval:
<http://octave.1599824.n4.nabble.com/file/n4663954/cylinder.png>
This is the part of my code that should solve these equations:
/f=@(x)(x*(besselj(1, x))/(besselj(0, x))-Bi);
zeta=[1:N]';
C=zeta;
for n=1:N
zeta(n)=fsolve(f, ((n-1)*pi+pi/4));
C(n)=(2/zeta(n))*(besselj(1, zeta(n)))/((besselj(0,
zeta(n)))^2+(besselj(1,
zeta(n)))^2);
end
i=0;
theta=zeta;
theta_star=Fo*0;
T=Fo*0;
for t=0:t_step:tmax
i=i+1;
for n=1:N
theta(n)=C(n)*exp(-(zeta(n))^2*Fo(i));
theta_star(i)=theta_star(i)+theta(n);
end
T(i)=theta_star(i)*(Ti-Tf)+Tf;
end/
Where /t/ is the time variable and /Fo(i)=/:
/t=[0:t_step:tmax];
Fo=(Alpha.*t/Lc^2)';/
with /Alpha/ known.
The resulting curve /T=f(t)/ apparently is wrong. Does anybody find anything
wrong?
--
View this message in context:
http://octave.1599824.n4.nabble.com/Transcendental-equation-tp4663840p4663954.html
Sent from the Octave - General mailing list archive at Nabble.com.