麻烦大神指点一下，非常感谢！

DIOADE

计算一个方程里的参数，公式如图，为什么一直出错，运行不出来QAQ，麻烦大神指点一下，非常感谢！
* |& z3 ^# w" L: H, O0 j
) y- }- g! Y" Y- P" hf1=@(x,a)(x.^(-a)-1)./a;
9 p1 [) ^, V3 Af2=@(x,a)diff(f1,x);
% ~0 K: f( ]' D; e0 S5 `5 [f3=@(x,a)f1(x,a)./f2(x,a);
! R3 o( g/ M5 zf4=@(x,a)integral(@(x)f3(x,a),x,0,1);
f5=@(a)1+4*f4(x,a)-6/2475;
a=fsolve(f5);

A-Lin

clear; clc; close all;
* N* `# E) O- vsyms t alpha real
syms theta(t)
theta = (t^(-alpha) - 1)/alpha;
* ~* ~, t: `, ?6 t, p3 d' @Dtheta = diff( theta, t );
; W+ d  i* w9 `3 ^tau = 4*int( theta/Dtheta, t, 0, 1 ) + 1;
4 ~& t1 ^; h) ^Alpha = solve([tau == 6/2457 ],[ alpha ])
RHS_LHS = 4*int( ( (t^(-Alpha) - 1)/Alpha )/diff( ( (t^(-Alpha) - 1)/Alpha ), t ), t, 0, 1 ) + 1 - 6/2457
" p" K, o$ O8 w: l  W
3 |5 _) g) f) L* S: }) W4 A! P1 ^Alpha =
( d, r0 R! s. K5 U: B$ @4/817
' T4 W" \  p' v
RHS_LHS =
0