Why Do I Find Wrong Root?

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.

Answers (0)

This question is closed.

Asked:

on 27 Mar 2013

Community Treasure Hunt

Find the treasures in MATLAB Central and discover how the community can help you!

Start Hunting!