Why Do I Find Wrong Root?
Info
This question is closed. Reopen it to edit or answer.
Show older comments
My .m file is
function [ ] = deneme ( )
syms x;
a=6.22*10^-5
b=1/[2.2093e+03/(x*8.6400e-05)+2.6876e+03*0.35*10^-3/(237*9.8800e-05)+2.7910e+03/(0.915*266.15*5.1746e-04)]
c=b/a
d=1-exp(-c)
e=65.8316
f=13.6797
g=d*a*(e-f)
h=g/(2*8.6400e-05)
j=2.46
k=204
l=3.9044
m=0.0104
n=0.0936
o=1.8*10^-3
p=196.7505
r=0.2760830447e-3-0.3578069421e-5*j+2.882405796*10^(-8)*j^2-2.421711925*10^(-10)*j^3+9.472074611*10^(-13)*j^4-2.340978515*10^(-14)*j^5+9.704212241*10^(-16)*j^6-4.908365400*10^(-18)*j^7-2.350097753*10^(-19)*j^8+3.564616924*10^(-21)*j^9-1.471394557*10^(-23)*j^10
s=k*(1-m)*1.8*10^-3/r
t=0.23e-1*s^.8*l^.4*n/o
u=(h/(k*p))
v=solve( x-t*230*u^0.5, x)
end
In system, all symbols is constant except for x. the true roots must be 1997 and 0 but after the calculations i find
v=1.8276429115700991713526650261042e-40 - 1.7442364916156394530764545262007e-48*i roots.
I am waiting for your solutions.
Sincere thanks.
1 Comment
Walter Roberson
on 27 Mar 2013
Answers (0)
This question is closed.
Community Treasure Hunt
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!