-
Notifications
You must be signed in to change notification settings - Fork 1
/
MN.m
63 lines (62 loc) · 10 KB
/
MN.m
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
function M = MN(N, b)
% MN(N, b) lists all components of N-th order structure tensors corresponding
% to fibre density "b". Includes tensors up to "N"=6. Results for small values of b have to taken with care as they might involve severe numerical errors.
switch N
case 1
M(1)=1/2+ 1/(8*b) -1/4 *sqrt(2/(pi*b)) *exp(2*b)/erfi(sqrt(2*b));
case 2
M(1) = (1/16).*(2.*pi).^(-1/2).*erfi(2.^(1/2).*b.^(1/2)).^(-1).*(4.*b.^(1/2).*((-3)+4.*b).*exp(1).^(2.*b)+3.*(2.*pi).^(1/2).*erfi(2.^(1/2).*b.^(1/2)));
M(2) = (1/32).*((-3)+(-4).*b+6.*b.^(1/2).*exp(1).^(2.*b).*(2.*pi.^(-1)).^(1/2).*erfi(2.^(1/2).*b.^(1/2)).^(-1));
M(3) = (3/128).*(2.*pi).^(-1/2).*erfi(2.^(1/2).*b.^(1/2)).^(-1).*((-4).*b.^(1/2).*(3+4.*b).*exp(1).^(2.*b)+(3+8.*b+16.*b.^2).*(2.*pi).^(1/2).*erfi(2.^(1/2).*b.^(1/2)));
M(4) = (1/128).*(2.*pi).^(-1/2).*erfi(2.^(1/2).*b.^(1/2)).^(-1).*((-4).*b.^(1/2).*(3+4.*b).*exp(1).^(2.*b)+(3+8.*b+16.*b.^2).*(2.*pi).^(1/2).*erfi(2.^(1/2).*b.^(1/2)));
M=M.*b.^(-2);
case 3
M(1)=(1/64).*(4.*b.^(1/2).*(15+(-20).*b+16.*b.^2).*exp(1).^(2.*b)+(-15).*(2.*pi).^(1/2).*erfi(2.^(1/2).*b.^(1/2)));
M(2)=(1/128).*(4.*b.^(1/2).*((-15)+8.*b).*exp(1).^(2.*b)+3.*(5+4.*b).*(2.*pi).^(1/2).*erfi(2.^(1/2).*b.^(1/2)));
M(3)=(-3/512).*((-4).*b.^(1/2).*(15+4.*b).*exp(1).^(2.*b)+(15+24.*b+16.*b.^2).*(2.*pi).^(1/2).*erfi(2.^(1/2).*b.^(1/2)));
M(4)=(-1/512).*((-4).*b.^(1/2).*(15+4.*b).*exp(1).^(2.*b)+(15+24.*b+16.*b.^2).*(2.*pi).^(1/2).*erfi(2.^(1/2).*b.^(1/2)));
M(5)=(5/1024).*((-4).*b.^(1/2).*(15+16.*b+16.*b.^2).*exp(1).^(2.*b)+(15+36.*b+48.*b.^2+64.*b.^3).*(2.*pi).^(1/2).*erfi(2.^(1/2).*b.^(1/2)));
M(6)=(1/1024).*((-4).*b.^(1/2).*(15+16.*b+16.*b.^2).*exp(1).^(2.*b)+(15+36.*b+48.*b.^2+64.*b.^3).*(2.*pi).^(1/2).*erfi(2.^(1/2).*b.^(1/2)));
M=M.*b.^(-3).*(2.*pi).^(-1/2).*erfi(2.^(1/2).*b.^(1/2)).^(-1);
case 4
M(1)=(1/256).*(4.*b.^(1/2).*((-105)+140.*b+(-112).*b.^2+64.*b.^3).*exp(1).^(2.*b)+105.*(2.*pi).^(1/2).*erfi(2.^(1/2).*b.^(1/2)));
M(2)=(1/512).*(4.*b.^(1/2).*(105+(-80).*b+32.*b.^2).*exp(1).^(2.*b)+(-15).*(7+4.*b).*(2.*pi).^(1/2).*erfi(2.^(1/2).*b.^(1/2)));
M(3)=(3/2048).*(20.*b.^(1/2).*((-21)+4.*b).*exp(1).^(2.*b)+3.*(35+40.*b+16.*b.^2).*(2.*pi).^(1/2).*erfi(2.^(1/2).*b.^(1/2)));
M(4)=(1/2048).*(20.*b.^(1/2).*((-21)+4.*b).*exp(1).^(2.*b)+3.*(35+40.*b+16.*b.^2).*(2.*pi).^(1/2).*erfi(2.^(1/2).*b.^(1/2)));
M(5)=(5/4096).*(4.*b.^(1/2).*(105+40.*b+16.*b.^2).*exp(1).^(2.*b)+(-1).*(105+180.*b+144.*b.^2+64.*b.^3).*(2.*pi).^(1/2).*erfi(2.^(1/2).*b.^(1/2)));
M(6)=(1/4096).*(4.*b.^(1/2).*(105+40.*b+16.*b.^2).*exp(1).^(2.*b)+(-1).*(105+180.*b+144.*b.^2+64.*b.^3).*(2.*pi).^(1/2).*erfi(2.^(1/2).*b.^(1/2)));
M(7)=(35/32768).*((-4).*b.^(1/2).*(105+100.*b+80.*b.^2+64.*b.^3).*exp(1).^(2.*b)+(105+240.*b+288.*b.^2+256.*b.^3+256.*b.^4).*(2.*pi).^(1/2).*erfi(2.^(1/2).*b.^(1/2)));
M(8)=(5/32768).*((-4).*b.^(1/2).*(105+100.*b+80.*b.^2+64.*b.^3).*exp(1).^(2.*b)+(105+240.*b+288.*b.^2+256.*b.^3+256.*b.^4).*(2.*pi).^(1/2).*erfi(2.^(1/2).*b.^(1/2)));
M(9)=(3/32768).*((-4).*b.^(1/2).*(105+100.*b+80.*b.^2+64.*b.^3).*exp(1).^(2.*b)+(105+240.*b+288.*b.^2+256.*b.^3+256.*b.^4).*(2.*pi).^(1/2).*erfi(2.^(1/2).*b.^(1/2)));
M=M.*b.^(-4).*(2.*pi).^(-1/2).*erfi(2.^(1/2).*b.^(1/2)).^(-1);
case 5
M(1)=(1/1024).*b.^(-5).*(2.*pi).^(-1/2).*erfi(2.^(1/2).*b.^(1/2)).^(-1).*(4.*b.^(1/2).*(945+(-1260).*b+1008.*b.^2+(-576).*b.^3+256.*b.^4).*exp(1).^(2.*b)+(-945).*(2.*pi).^(1/2).*erfi(2.^(1/2).*b.^(1/2)));
M(2)=(1/2048).*b.^(-5).*(2.*pi).^(-1/2).*erfi(2.^(1/2).*b.^(1/2)).^(-1).*(4.*b.^(1/2).*((-945)+840.*b+(-448).*b.^2+128.*b.^3).*exp(1).^(2.*b)+105.*(9+4.*b).*(2.*pi).^(1/2).*erfi(2.^(1/2).*b.^(1/2)));
M(3)=(3/8192).*b.^(-5).*(2.*pi).^(-1/2).*erfi(2.^(1/2).*b.^(1/2)).^(-1).*(4.*b.^(1/2).*(945+(-420).*b+128.*b.^2).*exp(1).^(2.*b)+(-15).*(63+56.*b+16.*b.^2).*(2.*pi).^(1/2).*erfi(2.^(1/2).*b.^(1/2)));
M(4)=(1/8192).*b.^(-5).*(2.*pi).^(-1/2).*erfi(2.^(1/2).*b.^(1/2)).^(-1).*(4.*b.^(1/2).*(945+(-420).*b+128.*b.^2).*exp(1).^(2.*b)+(-15).*(63+56.*b+16.*b.^2).*(2.*pi).^(1/2).*erfi(2.^(1/2).*b.^(1/2)));
M(5)=(15/16384).*b.^(-5).*(2.*pi).^(-1/2).*erfi(2.^(1/2).*b.^(1/2)).^(-1).*((-4).*b.^(1/2).*(315+16.*b.^2).*exp(1).^(2.*b)+(315+420.*b+240.*b.^2+64.*b.^3).*(2.*pi).^(1/2).*erfi(2.^(1/2).*b.^(1/2)));
M(6)=(3/16384).*b.^(-5).*(2.*pi).^(-1/2).*erfi(2.^(1/2).*b.^(1/2)).^(-1).*((-4).*b.^(1/2).*(315+16.*b.^2).*exp(1).^(2.*b)+(315+420.*b+240.*b.^2+64.*b.^3).*(2.*pi).^(1/2).*erfi(2.^(1/2).*b.^(1/2)));
M(7)=(35/131072).*b.^(-5).*(2.*pi).^(-1/2).*erfi(2.^(1/2).*b.^(1/2)).^(-1).*(4.*b.^(1/2).*(945+420.*b+208.*b.^2+64.*b.^3).*exp(1).^(2.*b)+(-1).*(945+1680.*b+1440.*b.^2+768.*b.^3+256.*b.^4).*(2.*pi).^(1/2).*erfi(2.^(1/2).*b.^(1/2)));
M(8)=(5/131072).*b.^(-5).*(2.*pi).^(-1/2).*erfi(2.^(1/2).*b.^(1/2)).^(-1).*(4.*b.^(1/2).*(945+420.*b+208.*b.^2+64.*b.^3).*exp(1).^(2.*b)+(-1).*(945+1680.*b+1440.*b.^2+768.*b.^3+256.*b.^4).*(2.*pi).^(1/2).*erfi(2.^(1/2).*b.^(1/2)));
M(9)=(3/131072).*b.^(-5).*(2.*pi).^(-1/2).*erfi(2.^(1/2).*b.^(1/2)).^(-1).*(4.*b.^(1/2).*(945+420.*b+208.*b.^2+64.*b.^3).*exp(1).^(2.*b)+(-1).*(945+1680.*b+1440.*b.^2+768.*b.^3+256.*b.^4).*(2.*pi).^(1/2).*erfi(2.^(1/2).*b.^(1/2)));
M(10)=(63/262144).*b.^(-5).*(2.*pi).^(-1/2).*erfi(2.^(1/2).*b.^(1/2)).^(-1).*((-4).*b.^(1/2).*(945+840.*b+608.*b.^2+384.*b.^3+256.*b.^4).*exp(1).^(2.*b)+(945+2100.*b+2400.*b.^2+1920.*b.^3+1280.*b.^4+1024.*b.^5).*(2.*pi).^(1/2).*erfi(2.^(1/2).*b.^(1/2)));
M(11)=(7/262144).*b.^(-5).*(2.*pi).^(-1/2).*erfi(2.^(1/2).*b.^(1/2)).^(-1).*((-4).*b.^(1/2).*(945+840.*b+608.*b.^2+384.*b.^3+256.*b.^4).*exp(1).^(2.*b)+(945+2100.*b+2400.*b.^2+1920.*b.^3+1280.*b.^4+1024.*b.^5).*(2.*pi).^(1/2).*erfi(2.^(1/2).*b.^(1/2)));
M(12)=(3/262144).*b.^(-5).*(2.*pi).^(-1/2).*erfi(2.^(1/2).*b.^(1/2)).^(-1).*((-4).*b.^(1/2).*(945+840.*b+608.*b.^2+384.*b.^3+256.*b.^4).*exp(1).^(2.*b)+(945+2100.*b+2400.*b.^2+1920.*b.^3+1280.*b.^4+1024.*b.^5).*(2.*pi).^(1/2).*erfi(2.^(1/2).*b.^(1/2)));
case 6
M(1)=(1/4096).*b.^(-6).*(2.*pi).^(-1/2).*erfi(2.^(1/2).*b.^(1/2)).^(-1).*(4.*b.^(1/2).*((-10395)+13860.*b+(-11088).*b.^2+6336.*b.^3+(-2816).*b.^4+1024.*b.^5).*exp(1).^(2.*b)+10395.*(2.*pi).^(1/2).*erfi(2.^(1/2).*b.^(1/2)));
M(2)=(1/8192).*b.^(-6).*(2.*pi).^(-1/2).*erfi(2.^(1/2).*b.^(1/2)).^(-1).*(4.*b.^(1/2).*(10395+(-10080).*b+6048.*b.^2+(-2304).*b.^3+512.*b.^4).*exp(1).^(2.*b)+(-945).*(11+4.*b).*(2.*pi).^(1/2).*erfi(2.^(1/2).*b.^(1/2)));
M(3)=(3/32768).*b.^(-6).*(2.*pi).^(-1/2).*erfi(2.^(1/2).*b.^(1/2)).^(-1).*(4.*b.^(1/2).*((-10395)+6300.*b+(-2688).*b.^2+512.*b.^3).*exp(1).^(2.*b)+105.*(99+72.*b+16.*b.^2).*(2.*pi).^(1/2).*erfi(2.^(1/2).*b.^(1/2)));
M(4)=(1/32768).*b.^(-6).*(2.*pi).^(-1/2).*erfi(2.^(1/2).*b.^(1/2)).^(-1).*(4.*b.^(1/2).*((-10395)+6300.*b+(-2688).*b.^2+512.*b.^3).*exp(1).^(2.*b)+105.*(99+72.*b+16.*b.^2).*(2.*pi).^(1/2).*erfi(2.^(1/2).*b.^(1/2)));
M(5)=(-15/65536).*b.^(-6).*(2.*pi).^(-1/2).*erfi(2.^(1/2).*b.^(1/2)).^(-1).*((-84).*b.^(1/2).*(165+(-40).*b+16.*b.^2).*exp(1).^(2.*b)+5.*(693+756.*b+336.*b.^2+64.*b.^3).*(2.*pi).^(1/2).*erfi(2.^(1/2).*b.^(1/2)));
M(6)=(-3/65536).*b.^(-6).*(2.*pi).^(-1/2).*erfi(2.^(1/2).*b.^(1/2)).^(-1).*((-84).*b.^(1/2).*(165+(-40).*b+16.*b.^2).*exp(1).^(2.*b)+5.*(693+756.*b+336.*b.^2+64.*b.^3).*(2.*pi).^(1/2).*erfi(2.^(1/2).*b.^(1/2)));
M(7)=(105/524288).*b.^(-6).*(2.*pi).^(-1/2).*erfi(2.^(1/2).*b.^(1/2)).^(-1).*((-4).*b.^(1/2).*(3465+420.*b+336.*b.^2+64.*b.^3).*exp(1).^(2.*b)+(3465+5040.*b+3360.*b.^2+1280.*b.^3+256.*b.^4).*(2.*pi).^(1/2).*erfi(2.^(1/2).*b.^(1/2)));
M(8)=(15/524288).*b.^(-6).*(2.*pi).^(-1/2).*erfi(2.^(1/2).*b.^(1/2)).^(-1).*((-4).*b.^(1/2).*(3465+420.*b+336.*b.^2+64.*b.^3).*exp(1).^(2.*b)+(3465+5040.*b+3360.*b.^2+1280.*b.^3+256.*b.^4).*(2.*pi).^(1/2).*erfi(2.^(1/2).*b.^(1/2)));
M(9)=(9/524288).*b.^(-6).*(2.*pi).^(-1/2).*erfi(2.^(1/2).*b.^(1/2)).^(-1).*((-4).*b.^(1/2).*(3465+420.*b+336.*b.^2+64.*b.^3).*exp(1).^(2.*b)+(3465+5040.*b+3360.*b.^2+1280.*b.^3+256.*b.^4).*(2.*pi).^(1/2).*erfi(2.^(1/2).*b.^(1/2)));
M(10)=(63/1048576).*b.^(-6).*(2.*pi).^(-1/2).*erfi(2.^(1/2).*b.^(1/2)).^(-1).*(4.*b.^(1/2).*(10395+5040.*b+2688.*b.^2+1024.*b.^3+256.*b.^4).*exp(1).^(2.*b)+(-1).*(10395+18900.*b+16800.*b.^2+9600.*b.^3+3840.*b.^4+1024.*b.^5).*(2.*pi).^(1/2).*erfi(2.^(1/2).*b.^(1/2)));
M(11)=(7/1048576).*b.^(-6).*(2.*pi).^(-1/2).*erfi(2.^(1/2).*b.^(1/2)).^(-1).*(4.*b.^(1/2).*(10395+5040.*b+2688.*b.^2+1024.*b.^3+256.*b.^4).*exp(1).^(2.*b)+(-1).*(10395+18900.*b+16800.*b.^2+9600.*b.^3+3840.*b.^4+1024.*b.^5).*(2.*pi).^(1/2).*erfi(2.^(1/2).*b.^(1/2)));
M(12)=(3/1048576).*b.^(-6).*(2.*pi).^(-1/2).*erfi(2.^(1/2).*b.^(1/2)).^(-1).*(4.*b.^(1/2).*(10395+5040.*b+2688.*b.^2+1024.*b.^3+256.*b.^4).*exp(1).^(2.*b)+(-1).*(10395+18900.*b+16800.*b.^2+9600.*b.^3+3840.*b.^4+1024.*b.^5).*(2.*pi).^(1/2).*erfi(2.^(1/2).*b.^(1/2)));
M(13)=(231/4194304).*b.^(-6).*(2.*pi).^(-1/2).*erfi(2.^(1/2).*b.^(1/2)).^(-1).*((-4).*b.^(1/2).*(10395+8820.*b+6048.*b.^2+3456.*b.^3+1792.*b.^4+1024.*b.^5).*exp(1).^(2.*b)+(10395+22680.*b+25200.*b.^2+19200.*b.^3+11520.*b.^4+6144.*b.^5+4096.*b.^6).*(2.*pi).^(1/2).*erfi(2.^(1/2).*b.^(1/2)));
M(14)=(21/4194304).*b.^(-6).*(2.*pi).^(-1/2).*erfi(2.^(1/2).*b.^(1/2)).^(-1).*((-4).*b.^(1/2).*(10395+8820.*b+6048.*b.^2+3456.*b.^3+1792.*b.^4+1024.*b.^5).*exp(1).^(2.*b)+(10395+22680.*b+25200.*b.^2+19200.*b.^3+11520.*b.^4+6144.*b.^5+4096.*b.^6).*(2.*pi).^(1/2).*erfi(2.^(1/2).*b.^(1/2)));
M(15)=(7/4194304).*b.^(-6).*(2.*pi).^(-1/2).*erfi(2.^(1/2).*b.^(1/2)).^(-1).*((-4).*b.^(1/2).*(10395+8820.*b+6048.*b.^2+3456.*b.^3+1792.*b.^4+1024.*b.^5).*exp(1).^(2.*b)+(10395+22680.*b+25200.*b.^2+19200.*b.^3+11520.*b.^4+6144.*b.^5+4096.*b.^6).*(2.*pi).^(1/2).*erfi(2.^(1/2).*b.^(1/2)));
M(16)=(5/4194304).*b.^(-6).*(2.*pi).^(-1/2).*erfi(2.^(1/2).*b.^(1/2)).^(-1).*((-4).*b.^(1/2).*(10395+8820.*b+6048.*b.^2+3456.*b.^3+1792.*b.^4+1024.*b.^5).*exp(1).^(2.*b)+(10395+22680.*b+25200.*b.^2+19200.*b.^3+11520.*b.^4+6144.*b.^5+4096.*b.^6).*(2.*pi).^(1/2).*erfi(2.^(1/2).*b.^(1/2)));
end