How to write the codes for this?

i am trying to get the value for b.

where

t = [89;132;202;263;321;362;421;473;575;663]
n=10;
syms b
[b] = solve((1/b) + sum(log(t))/10 - 1.*sum(t.^b).*log(t)/(sum(t.^b))==0)

i tried this code to find 'b' but i am not getting the value of 'b'. How can i get the value for 'b'?

Answer 1

[Ameer Hamza]

If you want to find value of 'b' for each value of 't' then try this

tv = [89;132;202;263;321;362;421;473;575;663];
syms b t
b_sol = solve((1/b) + sum(log(t))/10 - 1/(sum(t.^b)/10).*sum(t.^b).*log(t)==0, b);
b_vals = subs(b_sol, t, tv)

[MD. Rokibujjaman sovon]

i am not looking for each t. 'b' will have a single value and it would be 2.0924

[Ameer Hamza]

Can you show your equation in latex form?

[MD. Rokibujjaman sovon]

i added the question.

[Ameer Hamza]

You were also dividing the term in the denominator by 10. Try the following code

t = [89;132;202;263;321;362;421;473;575;663]
n=10;
syms b
b_sol = solve((1/b) + sum(log(t))/10 - 1/(sum(t.^b)).*sum((t.^b).*log(t))==0)

Result

b_sol =
2.0923855034345368572888674897445

[MD Rokibujjaman Sabuj]

Thanks.

Answer 2

[Walter Roberson]

Your t is a vector.

[b] = solve((1/b) + sum(log(t))/10 - 1/(sum(t.^b)/10).*sum(t.^b).*log(t)==0)

You have log(t) there at the end, in a context that is not inside a sum() like the rest of the places you use t. Therefore you are generating a vector of equations, and you are asking solve() to find a single t that solves all of the equations simultaneously.

You need to re-check your last term. For example there is a numeric solution for

vpasolve((1/b) + sum(log(t))/10 - 1/(sum(t.^b)/10).*sum(t.^b.*log(t))==0, 0.02)