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

DIOADE

! r9 {8 X8 }5 @( [4 W0 v" e, X计算一个方程里的参数，公式如图，为什么一直出错，运行不出来QAQ，麻烦大神指点一下，非常感谢！

f1=@(x,a)(x.^(-a)-1)./a;
. K* u# V3 n* T4 _* d* Yf2=@(x,a)diff(f1,x);
: A5 X% f0 i! j9 q4 R) P7 lf3=@(x,a)f1(x,a)./f2(x,a);
f4=@(x,a)integral(@(x)f3(x,a),x,0,1);
* O  G+ m1 D% b: C5 o- e% Wf5=@(a)1+4*f4(x,a)-6/2475;
& H2 S1 \" K% I' O( pa=fsolve(f5);
, U$ i  b2 L# X# |0 V$ T- u/ s

A-Lin

clear; clc; close all;
* Z/ l3 Z. ^( o# vsyms t alpha real
1 O, R% R6 T2 E+ m. P+ s% {* K7 Xsyms theta(t)
9 g. v. [7 t. T5 H1 k2 y% gtheta = (t^(-alpha) - 1)/alpha;
Dtheta = diff( theta, t );
7 I8 L. O( x0 j1 o( f* _. b8 mtau = 4*int( theta/Dtheta, t, 0, 1 ) + 1;
+ W' N3 ]+ E. p0 ~. o3 ]  Q+ kAlpha = 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
1 A1 ~/ J# ?+ X, K. h! H
! g$ A3 E2 H* p9 ?9 y9 n: y% g# k; U0 zAlpha =
4/817
7 ^/ M, a; A* O' w& g; r
1 p6 I) A! a. g6 pRHS_LHS =
0