solving an equation for a variable that has changing varibles in it

Question

joshua payne am 9 Okt. 2022

0
Verknüpfen

Direkter Link zu dieser Frage

https://de.mathworks.com/matlabcentral/answers/1820615-solving-an-equation-for-a-variable-that-has-changing-varibles-in-it

Kommentiert: Ghazwan am 9 Okt. 2022

im trying to solve for M1 based off a changing Pr

index1=1:304;
index2=697:1000;
index3=305:696;
figure
plot(t(index1),M(index1), 'b')
hold on
plot(t(index2),M(index2), 'b')
%plot(t(index3),M(index3), 'r')
M2=(((M1).^(2))+(2/(gamma-1)))./((((2*gamma)/(gamma-1)).*(M1).^2)-1)
pr=Po(index3)./P(index3)
eqn=((1+((gamma-1)/2).*(M2).^2)).^(gamma/(gamma-1)).*(((gamma+1).*(M1).^2)./(2+(gamma-1).*(M1).^2))==pr
superM=vpasolve(eqn,M, 'real')
plot(t(index3),superM, 'r')

im trying to get only real vaules for M1

8 Kommentare
6 ältere Kommentare anzeigen6 ältere Kommentare ausblenden

joshua payne am 9 Okt. 2022

In MATLAB Online öffnen

clc
clear all
t=[0.00
20
40
60
80
01
21
41
61
81
01
21
41
62
82
02
22
42
62
82
02
23
43
63
83
03
23
43
63
83
04
24
44
64
84
04
24
44
65
85
05
25
45
65
85
05
26
46
66
86
06
26
46
66
86
07
27
47
67
87
07
27
47
68
88
08
28
48
68
88
08
29
49
69
89
09
29
49
69
89
10
30
50
70
90
10
30
50
71
91
11
31
51
71
91
11
32
52
72
92
12
32
52
72
92
13
33
53
73
93
13
33
53
74
94
14
34
54
74
94
14
35
55
75
95
15
35
55
75
95
16
36
56
76
96
16
36
56
77
97
17
37
57
77
97
17
38
58
78
98
18
38
58
78
98
19
39
59
79
99
19
39
59
80
00
20
40
60
80
00
20
41
61
81
01
21
41
61
81
02
22
42
62
82
02
22
42
62
83
03
23
43
63
83
03
23
44
64
84
04
24
44
64
84
05
25
45
65
85
05
25
45
65
86
06
26
46
66
86
06
26
47
67
87
07
27
47
67
87
08
28
48
68
88
08
28
48
68
89
09
29
49
69
89
09
29
50
70
90
10
30
50
70
90
11
31
51
71
91
11
31
51
71
92
12
32
52
72
92
12
32
53
73
93
13
33
53
73
93
14
34
54
74
94
14
34
54
74
95
15
35
55
75
95
15
35
56
76
96
16
36
56
76
96
17
37
57
77
97
17
37
57
77
98
18
38
58
78
98
18
38
59
79
99
19
39
59
79
99
20
40
60
80
00
20
40
60
80
01
21
41
61
81
01
21
41
62
82
02
22
42
62
82
02
23
43
63
83
03
23
43
63
83
04
24
44
64
84
04
24
44
65
85
05
25
45
65
85
05
26
46
66
86
06
26
46
66
86
07
27
47
67
87
07
27
47
68
88
08
28
48
68
88
08
29
49
69
89
09
29
49
69
89
10
30
50
70
90
10
30
50
71
91
11
31
51
71
91
11
32
52
72
92
12
32
52
72
92
13
33
53
73
93
13
33
53
74
94
14
34
54
74
94
14
35
55
75
95
15
35
55
75
95
16
36
56
76
96
16
36
56
77
97
17
37
57
77
97
17
38
58
78
98
18
38
58
78
98
19
39
59
79
99
19
39
59
80
00
20
40
60
80
00
20
41
61
81
01
21
41
61
81
02
22
42
62
82
02
22
42
62
83
03
23
43
63
83
03
23
44
64
84
04
24
44
64
84
05
25
45
65
85
05
25
45
65
86
06
26
46
66
86
06
26
47
67
87
07
27
47
67
87
08
28
48
68
88
08
28
48
68
89
09
29
49
69
89
09
29
50
70
90
10
30
50
70
90
11
31
51
71
91
11
31
51
71
92
12
32
52
72
92
12
32
53
73
93
13
33
53
73
93
14
34
54
74
94
14
34
54
74
95
15
35
55
75
95
15
35
56
76
96
16
36
56
76
96
17
37
57
77
97
17
37
57
77
98
18
38
58
78
98
18
38
59
79
99
19
39
59
79
99
20
40
60
80
00
20
40
60
80
01
21
41
61
81
01
21
41
62
82
02
22
42
62
82
02
23
43
63
83
03
23
43
63
83
04
24
44
64
84
04
24
44
65
85
05
25
45
65
85
05
26
46
66
86
06
26
46
66
86
07
27
47
67
87
07
27
47
68
88
08
28
48
68
88
08
29
49
69
89
09
29
49
69
89
10
30
50
70
90
10
30
50
71
91
11
31
51
71
91
11
32
52
72
92
12
32
52
72
92
13
33
53
73
93
13
33
53
74
94
14
34
54
74
94
14
35
55
75
95
15
35
55
75
95
16
36
56
76
96
16
36
56
77
97
17
37
57
77
97
17
38
58
78
98
18
38
58
78
98
19
39
59
79
99
19
39
59
80
00
20
40
60
80
00
20
41
61
81
01
21
41
61
81
02
22
42
62
82
02
22
42
62
83
03
23
43
63
83
03
23
44
64
84
04
24
44
64
84
05
25
45
65
85
05
25
45
65
86
06
26
46
66
86
06
26
47
67
87
07
27
47
67
87
08
28
48
68
88
08
28
48
68
89
09
29
49
69
89
09
29
50
70
90
10
30
50
70
90
11
31
51
71
91
11
31
51
71
92
12
32
52
72
92
12
32
53
73
93
13
33
53
73
93
14
34
54
74
94
14
34
54
74
95
15
35
55
75
95
15
35
56
76
96
16
36
56
76
96
17
37
57
77
97
17
37
57
77
98
18
38
58
78
98
18
38
59
79
99
19
39
59
79
99
20
40
60
80
00];
P=[101000.00
100821.64
100643.29
100464.93
100286.57
100108.22
86
50
15
79
43
08
72
36
01
65
29
94
58
22
87
51
15
80
44
08
73
37
01
66
30
94
59
23
87
52
16
80
44
09
73
37
02
66
30
95
59
23
88
52
16
81
45
09
74
38
02
67
31
95
60
24
88
53
17
81
46
10
74
39
03
67
32
96
60
25
89
53
18
82
46
11
75
39
04
68
32
97
61
25
90
54
18
83
47
11
76
40
04
69
33
97
62
26
90
55
19
83
47
12
76
40
05
69
33
98
62
26
91
55
19
84
48
12
77
41
05
70
34
98
63
27
91
56
20
84
49
13
77
42
06
70
35
99
63
28
92
56
21
85
49
14
78
42
07
71
35
00
64
28
93
57
21
86
50
14
79
43
07
72
36
00
65
29
93
58
22
86
51
15
79
43
08
72
36
01
65
29
94
58
22
87
51
15
80
44
08
73
37
01
66
30
94
59
23
87
52
16
80
45
09
73
38
02
66
31
95
59
24
88
52
17
81
45
10
74
38
03
67
31
96
60
24
89
53
17
82
46
10
75
39
03
68
32
96
61
25
89
54
18
82
46
11
75
39
04
68
32
97
61
25
90
54
18
83
47
11
76
40
04
69
33
97
62
26
90
55
19
83
48
12
76
41
05
69
34
98
62
27
91
55
20
84
48
13
77
41
06
70
34
99
63
27
92
56
20
85
49
13
78
42
06
71
35
99
64
28
92
57
21
85
49
14
78
42
07
71
35
00
64
28
93
57
21
86
50
14
79
43
07
72
36
00
65
29
93
58
22
86
51
15
79
44
08
72
37
01
65
30
94
58
23
87
51
16
80
44
09
73
37
02
66
30
95
59
23
88
52
16
81
45
09
74
38
02
67
31
95
60
24
88
53
17
81
45
10
74
38
03
67
31
96
60
24
89
53
17
82
46
10
75
39
03
68
32
96
61
25
89
54
18
82
47
11
75
40
04
68
33
97
61
26
90
54
19
83
47
12
76
40
05
69
33
98
62
26
91
55
19
84
48
12
77
41
05
70
34
98
63
27
91
56
20
84
48
13
77
41
06
70
34
99
63
27
92
56
20
85
49
13
78
42
06
71
35
99
64
28
92
57
21
85
50
14
78
43
07
71
36
00
00
36
71
07
43
78
14
50
85
21
57
92
28
64
99
35
71
06
42
78
13
49
85
20
56
92
27
63
99
34
70
06
41
77
13
48
84
20
56
91
27
63
98
34
70
05
41
77
12
48
84
19
55
91
26
62
98
33
69
05
40
76
12
47
83
19
54
90
26
61
97
33
68
04
40
75
11
47
82
18
54
89
25
61
96
32
68
03
39
75
10
46
82
17
53
89
24
60
96
31
67
03
38
74
10
45
81
17
53
88
24
60
95
31
67
02
38
74
09
45
81
16
52
88
23
59
95
30
66
02
37
73
09
44
80
16
51
87
23
58
94
30
65
01
37
72
08
44
79
15
51
86
22
58
93
29
65
00
36
72
07
43
79
14
50
86
21
57
93
28
64
00
35
71
07
42
78
14
49
85
21
57
92
28
64
99
35
71
06
42
78
13
49
85
20
56
92
27
63
99
34
70
06
41
77
13
48
84
20
55
91
27
62
98
34
69
05
41
76
12
48
83
19
55
90
26
62
97
33
69
04
40
76
11
47
83
18
54
90
25
61
97
32
68
04
39
75
11
46
82
18
54
89
25
61
96
32
68
03
39
75
10
46
82
17
53
89
24
60
96
31
67
03
38
74
10
45
81
17
52
88
24
59
95
31
66
02
38
73
09
45
80
16
52
87
23
59
94
30
66
01
37
73
08
44
80
15
51
87
22
58
94
29
65
01
36
72
08
43
79
15
51
86
22
58
93
29
65
00
36
72
07
43
79
14
50
86
21
57
93
28
64
00
35
71
07
42
78
14
49
85
21
56
92
28
63
99
35
70
06
42
77
13
49
84
20
56
91
27
63
98
34
70
05
41
77
12
48
84
19
55
91
26
62
98
33
69
05
40
76
12
47
83
19
55
90
26
62
97
33
69
04
40
76
11
47
83
18
54
90
25
61
97
32
68
04
39
75
11
46
82
18
53
89
25
60
96
32
67
03
39
74
10
46
81
17
53
88
24
60
95
31
67
02
38
74
09
45
81
16
52
88
23
59
95
30
66
02
37
73
09
44
80
16
52
87
23
59
94
30
66
01
37
73
08
44
80
15
51
87
22
58
94
29
65
01
36
72
08
43
79
15
50
86
100108.22
100286.57
100464.93
100643.29
100821.64
101000.00];
Po=[101034.72
100857.57
100680.46
100503.40
100326.37
100149.39
45
56
71
91
16
47
82
23
70
22
80
44
15
91
75
65
62
67
78
98
25
60
03
55
15
85
63
51
48
56
74
02
41
91
52
25
09
07
16
39
74
24
87
64
56
63
86
24
79
50
38
43
67
08
69
48
48
67
07
69
52
57
85
36
11
10
35
85
61
63
94
52
39
55
01
77
86
26
99
05
46
21
33
81
66
90
53
55
99
83
11
81
96
57
63
17
19
69
71
23
27
85
97
64
88
70
11
11
73
97
85
38
56
42
97
22
18
86
29
47
42
14
67
01
17
17
03
76
38
90
34
71
04
33
61
89
19
52
91
37
92
58
36
29
39
66
15
85
79
00
49
28
39
85
67
88
49
54
03
00
47
45
97
06
74
02
94
52
78
74
44
89
13
17
04
77
38
89
35
76
17
58
04
57
20
95
85
93
23
76
55
64
06
83
98
55
56
04
02
54
63
31
62
58
24
62
75
67
40
98
45
83
16
47
80
17
62
19
91
81
92
29
94
91
24
95
08
68
77
38
56
34
75
82
61
13
42
53
48
31
06
76
44
15
92
77
76
91
25
83
67
82
30
15
40
09
25
91
12
89
26
26
93
30
39
24
87
33
63
80
87
87
83
77
72
69
72
83
04
36
83
46
27
27
49
93
62
56
77
25
02
09
45
12
10
39
99
89
10
55
02
26
04
13
30
32
96
02
27
50
50
05
95
99
98
69
95
53
26
93
34
30
62
10
100003.56
100458.80
100913.63
101367.87
101821.32
102273.80
102725.13
103175.11
103623.56
104070.30
104515.14
104957.90
105398.40
105836.45
106271.88
106704.50
107134.13
107560.60
107983.72
108403.32
108819.22
109231.24
109639.20
110042.94
110442.28
110837.04
111227.06
111612.15
111992.15
112366.89
112736.20
113099.91
113457.85
113809.86
114155.78
114495.43
114828.65
115155.30
115475.19
115788.18
116094.11
116392.81
116684.15
116967.95
117244.07
117512.36
117772.68
118024.86
118268.77
118504.27
118731.20
118949.44
119158.83
119359.26
119550.57
119732.64
119905.33
120068.53
120222.10
120365.91
120499.85
120623.80
120737.63
120841.24
120934.51
121017.33
121089.59
121151.19
121202.02
121241.99
121271.00
121288.95
121295.75
121291.32
121275.56
121248.40
121209.75
121159.53
121097.68
121024.12
120938.78
120841.60
120732.51
120611.46
120478.39
120333.24
120175.98
120006.55
119824.92
119631.03
119424.86
119206.37
118975.54
118732.34
118476.75
118208.74
117928.31
117635.44
117330.13
117012.36
116682.14
116339.48
115984.37
115616.83
115236.87
114844.51
114439.77
114022.67
113593.24
113151.52
112697.53
112231.32
111752.92
111262.39
110759.77
110245.12
109718.49
109179.95
108629.55
108067.37
107493.47
106907.93
106310.82
105702.24
105082.25
104450.96
103808.45
103154.82
102490.15
101814.57
101128.17
100431.05
33
13
56
73
79
84
03
47
32
70
76
64
48
45
68
33
56
52
39
31
47
02
14
99
76
62
74
31
51
52
53
72
29
41
29
11
07
36
36
07
11
29
41
29
72
53
52
51
31
74
62
76
99
14
02
47
31
39
52
56
33
68
45
48
64
76
70
32
47
03
84
79
73
56
13
33
100431.05
101128.17
101814.57
102490.15
103154.82
103808.45
104450.96
105082.25
105702.24
106310.82
106907.93
107493.47
108067.37
108629.55
109179.95
109718.49
110245.12
110759.77
111262.39
111752.92
112231.32
112697.53
113151.52
113593.24
114022.67
114439.77
114844.51
115236.87
115616.83
115984.37
116339.48
116682.14
117012.36
117330.13
117635.44
117928.31
118208.74
118476.75
118732.34
118975.54
119206.37
119424.86
119631.03
119824.92
120006.55
120175.98
120333.24
120478.39
120611.46
120732.51
120841.60
120938.78
121024.12
121097.68
121159.53
121209.75
121248.40
121275.56
121291.32
121295.75
121288.95
121271.00
121241.99
121202.02
121151.19
121089.59
121017.33
120934.51
120841.24
120737.63
120623.80
120499.85
120365.91
120222.10
120068.53
119905.33
119732.64
119550.57
119359.26
119158.83
118949.44
118731.20
118504.27
118268.77
118024.86
117772.68
117512.36
117244.07
116967.95
116684.15
116392.81
116094.11
115788.18
115475.19
115155.30
114828.65
114495.43
114155.78
113809.86
113457.85
113099.91
112736.20
112366.89
111992.15
111612.15
111227.06
110837.04
110442.28
110042.94
109639.20
109231.24
108819.22
108403.32
107983.72
107560.60
107134.13
106704.50
106271.88
105836.45
105398.40
104957.90
104515.14
104070.30
103623.56
103175.11
102725.13
102273.80
101821.32
101367.87
100913.63
100458.80
100003.56
10
62
30
34
93
26
53
95
69
98
99
95
05
50
50
27
02
96
32
30
13
04
26
02
55
10
89
99
39
10
12
45
09
02
25
77
56
62
93
49
27
27
46
83
36
04
83
72
69
72
77
83
87
87
80
63
33
87
24
39
30
93
26
26
89
12
91
25
09
40
15
30
82
67
83
25
91
76
77
92
15
44
76
06
31
48
53
42
13
61
82
75
34
56
38
77
68
08
95
24
91
94
29
92
81
91
19
62
17
80
47
16
83
45
98
40
67
75
62
24
58
62
31
63
54
02
04
56
55
98
83
06
64
55
76
23
93
85
95
20
57
04
58
17
76
35
89
38
77
04
17
13
89
44
74
78
52
94
02
74
06
97
45
47
00
03
54
49
88
67
85
39
28
49
00
79
85
15
66
39
29
36
58
92
37
91
52
19
89
61
33
04
71
34
90
38
76
03
17
17
01
67
14
42
47
29
86
18
22
97
42
56
38
85
97
73
11
11
70
88
64
97
85
27
23
71
69
19
17
63
57
96
81
11
83
99
55
53
90
66
81
33
21
46
05
99
26
86
77
01
55
39
52
94
63
61
85
35
10
11
36
85
57
52
69
07
67
48
48
69
08
67
43
38
50
79
24
86
63
56
64
87
24
74
39
16
07
09
25
52
91
41
02
74
56
48
51
63
85
15
55
03
60
25
98
78
67
62
65
75
91
15
44
80
22
70
23
82
47
16
91
71
56
45
100149.39
100326.37
100503.40
100680.46
100857.57
101034.72];
gamma=1.4;
M=((((Po./P).^((gamma-1)./gamma)-1)*2)/(gamma-1)).^(1/2);
Msonic=find(abs(M-1)<.005)
values=M(305)
values=M(696)
   
Po1=P.*(1+((gamma-1)/2).*M.^2);
M2=(((M.^2)+(2/(gamma-1)))/(((2*gamma)/(gamma-1))-1));
%Po2=
index1=1:304;
index2=697:1000;
index3=305:696;
figure
plot(t(index1),M(index1), 'b')
hold on
plot(t(index2),M(index2), 'b')
%plot(t(index3),M(index3), 'r')
%M2=((((M1).^(2))+(2/(gamma-1)))./((((2*gamma)/(gamma-1)).*(M1).^2)-1))
pr=Po(index3)./P(index3)
syms M1
eqn=((1+((gamma-1)/2).*(((((M1).^(2))+(2/(gamma-1)))./((((2*gamma)/(gamma-1)).*(M1).^2)-1))).^2)).^(gamma/(gamma-1)).*(((gamma+1).*(M1).^2)./(2+(gamma-1).*(M1).^2))==pr
superM=vpasolve(eqn,M1, 'real')
plot(t(index3),superM, 'r')
hold off