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