Can u help to write the code scilab to solve the transcendantal equation by Newton Raphson?
JOSE KANTOLE
jkantole at gmail.com
Tue Dec 28 01:39:27 CET 2010
Hello
I am fairly new to using Scilab but need to solve a transcendental
equation. I would like to write the code to solve this equation
.The value sought is F0
Given data: emin=0.688;
NTU= 1.488 ;
B=0.334;
p=-ln(1-emin);
x=( NTU-p) *2/B;
y=exp(x) ;
F0=(y-1)/(1-y*(1-emin));
If F0>10; or F0<0, then
y=exp(2*x/p);
This is the first approximation F0=sqrt(y/(1-emin))
My problem is to solve the problem using Raphson-Newton method
F1=(1-emin)*F0;
p0=ln(1+F0);
p20=p0*p0;
G0= p0*(1+17*p20/450)/(1+p20/100)+p20/4
p=ln(1+ F1);
p2=p*p;
G1= p*(1+17*p2/450)/(1+p2/450)+p2/4;
x=-ln (1-emin)+B/2*(G0-G1)-NTU;
y=B/2*(p0/F0-(1-emin)*p/(1+F1)); this is the derivative of x
Fon=F0-x/y ( Newton-Raphson)
Hopefully someone can see what is needed to get things working. Thanks
in advance all.
More information about the users
mailing list