Is it possible to solve the equation for x?

I am trying to find the value for x and all parameters values are given.

any suggestion or help will be great. Thank you

clear all
syms x
a(1) = 0.89;
a(2) = 6.7;
a(3) = 26.12;
a(4) = 0.07;
a(5) = 78.05;
a(6) = 34.01;
a(7) = 0.33;
a(8) = 22.49;
a(9) = -0.06;
T = 50;
fun = a(1)./a(2).*exp(-(x-a(3))./a(2)).*exp(-exp(-(x-a(3))./a(2)))+a(4).*a(5).*a(6).^a(5).*x.^(-1-a(5)).*1./exp((a(6)./x).^a(5))+(1-a(4)-a(1))./a(7).*(1-((x-a(8))./a(7)).*a(9)).^(1./a(9)-1).*1./exp((1-((x-a(8))./a(7)).*a(9)).^(1./a(9)))-(1-1./T) == 0;
sol = fsolve(fun)

Answer 1

[Stephan]

syms x
a(1) = 0.89;
a(2) = 6.7;
a(3) = 26.12;
a(4) = 0.07;
a(5) = 78.05;
a(6) = 34.01;
a(7) = 0.33;
a(8) = 22.49;
a(9) = -0.06;
T = 50;
fun = a(1)./a(2).*exp(-(x-a(3))./a(2)).*exp(-exp(-(x-a(3))./a(2)))+a(4).*a(5).*a(6).^a(5).*x.^(-1-a(5)).*1./exp((a(6)./x).^a(5))+(1-a(4)-a(1))./a(7).*(1-((x-a(8))./a(7)).*a(9)).^(1./a(9)-1).*1./exp((1-((x-a(8))./a(7)).*a(9)).^(1./a(9)))-(1-1./T) == 0;
sol = solve(fun,x)

[hassan tahir]

Thank you for your prompt response. But it doesn't work.

It is showing error in sol = solve(fun,x)

[Stephan]

try:

syms x
a(1) = 0.89;
a(2) = 6.7;
a(3) = 26.12;
a(4) = 0.07;
a(5) = 78.05;
a(6) = 34.01;
a(7) = 0.33;
a(8) = 22.49;
a(9) = -0.06;
T = 50;
fun = a(1)./a(2).*exp(-(x-a(3))./a(2)).*exp(-exp(-(x-a(3))./a(2)))+a(4).*a(5).*a(6).^a(5).*x.^(-1-a(5)).*1./exp((a(6)./x).^a(5))+(1-a(4)-a(1))./a(7).*(1-((x-a(8))./a(7)).*a(9)).^(1./a(9)-1).*1./exp((1-((x-a(8))./a(7)).*a(9)).^(1./a(9)))-(1-1./T) == 0;
sol = vpasolve(fun,x)

[Star Strider]

Also:

sol = vpasolve(fun,x, 'random',1)

will discover several roots in different iterations.

[hassan tahir]

@Stephan and @Star Strider : I tried both codes and the answer is

20.35396878970819364858892254613 - 322.89084343234018555725784362549i

It is in the form of 'i', which isn't right

[Walter Roberson]

syms x real
a(1) = 0.89;
a(2) = 6.7;
a(3) = 26.12;
a(4) = 0.07;
a(5) = 78.05;
a(6) = 34.01;
a(7) = 0.33;
a(8) = 22.49;
a(9) = -0.06;
T = 50;
fun = a(1)./a(2).*exp(-(x-a(3))./a(2)).*exp(-exp(-(x-a(3))./a(2)))+a(4).*a(5).*a(6).^a(5).*x.^(-1-a(5)).*1./exp((a(6)./x).^a(5))+(1-a(4)-a(1))./a(7).*(1-((x-a(8))./a(7)).*a(9)).^(1./a(9)-1).*1./exp((1-((x-a(8))./a(7)).*a(9)).^(1./a(9)))-(1-1./T) == 0;
fplot(lhs(fun), [30 40])

It peaks at well below 0. There is no real-valued solution.

[hassan tahir]

Thank you Walter