|
clear; clc; close all;- E q' S: F# e
syms t alpha real
^7 }2 [; t5 ~$ E2 P- W; }3 usyms theta(t)% H5 D8 j1 r N5 g1 ^) ^' ?/ a
theta = (t^(-alpha) - 1)/alpha;5 n6 F0 }8 O4 T2 d, z! u4 a
Dtheta = diff( theta, t );7 f: Z* I, }7 E( g
tau = 4*int( theta/Dtheta, t, 0, 1 ) + 1;/ t d' v4 v% e8 v6 T
Alpha = solve([tau == 6/2457 ],[ alpha ])
+ i# f7 m& p" Y& [% \3 N3 vRHS_LHS = 4*int( ( (t^(-Alpha) - 1)/Alpha )/diff( ( (t^(-Alpha) - 1)/Alpha ), t ), t, 0, 1 ) + 1 - 6/2457 y' I' }4 F8 ~8 G7 ]. s9 ^
$ v2 n" U6 h6 Q, W4 I0 sAlpha =
3 R% b: U/ K' A1 ]& P2 j4/8172 b; f( [+ J3 \- `/ R$ ]: |$ d
1 P9 f* M2 A& a, xRHS_LHS =* b. e$ [( G! I" R+ M( S
0 |
|