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.

Réponses (0)

Cette question est clôturée.

Community Treasure Hunt

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

Start Hunting!

Translated by