|
|
clear; clc; close all;) b% K, f$ Z9 W" R+ _+ P7 ?
syms t alpha real
9 e) a' p$ B6 Q4 v* q" v1 jsyms theta(t)
& G7 l/ p6 I; Z8 g d+ G4 F. Atheta = (t^(-alpha) - 1)/alpha;
+ m3 x! r1 ]! m( y; ZDtheta = diff( theta, t );
% u1 }0 K( W, u* vtau = 4*int( theta/Dtheta, t, 0, 1 ) + 1;
/ W) b8 w3 A+ l2 \; E7 PAlpha = solve([tau == 6/2457 ],[ alpha ])/ e* W- Z, }5 G! v+ w
RHS_LHS = 4*int( ( (t^(-Alpha) - 1)/Alpha )/diff( ( (t^(-Alpha) - 1)/Alpha ), t ), t, 0, 1 ) + 1 - 6/2457
- y+ O/ |, W1 Q6 u0 b
( y/ k6 \5 _1 pAlpha =5 {% o( C, _. R: e
4/817) v: t' T0 A0 U6 p. i2 L
$ l4 M$ o4 j: \: w/ \" aRHS_LHS =
, O9 p5 j! P7 C% W0 |
|