clear; clc; close all;
; m( K. M: @, T2 Z6 }syms t alpha real
! B' L6 E0 t3 |1 J& Xsyms theta(t)
theta = (t^(-alpha) - 1)/alpha;
Dtheta = diff( theta, t );
tau = 4*int( theta/Dtheta, t, 0, 1 ) + 1;
Alpha = solve([tau == 6/2457 ],[ alpha ])
6 x$ q* w( y  B$ X! J  qRHS_LHS = 4*int( ( (t^(-Alpha) - 1)/Alpha )/diff( ( (t^(-Alpha) - 1)/Alpha ), t ), t, 0, 1 ) + 1 - 6/2457
: I  r' Y# K% ?/ ^( @; Z$ y" R7 @
Alpha =
! m% v6 e7 _% b+ X2 k! B- P4/817

RHS_LHS =
0 F8 {! P- m7 B0 z3 T- ^: d( z6 f0