Why Do I Find Wrong Root?
Infos
Cette question est clôturée. Rouvrir pour modifier ou répondre.
Afficher commentaires plus anciens
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 commentaire
Walter Roberson
le 27 Mar 2013
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!