-
Notifications
You must be signed in to change notification settings - Fork 2
Expand file tree
/
Copy pathExample.nb
More file actions
1403 lines (1347 loc) · 60 KB
/
Example.nb
File metadata and controls
1403 lines (1347 loc) · 60 KB
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
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
125
126
127
128
129
130
131
132
133
134
135
136
137
138
139
140
141
142
143
144
145
146
147
148
149
150
151
152
153
154
155
156
157
158
159
160
161
162
163
164
165
166
167
168
169
170
171
172
173
174
175
176
177
178
179
180
181
182
183
184
185
186
187
188
189
190
191
192
193
194
195
196
197
198
199
200
201
202
203
204
205
206
207
208
209
210
211
212
213
214
215
216
217
218
219
220
221
222
223
224
225
226
227
228
229
230
231
232
233
234
235
236
237
238
239
240
241
242
243
244
245
246
247
248
249
250
251
252
253
254
255
256
257
258
259
260
261
262
263
264
265
266
267
268
269
270
271
272
273
274
275
276
277
278
279
280
281
282
283
284
285
286
287
288
289
290
291
292
293
294
295
296
297
298
299
300
301
302
303
304
305
306
307
308
309
310
311
312
313
314
315
316
317
318
319
320
321
322
323
324
325
326
327
328
329
330
331
332
333
334
335
336
337
338
339
340
341
342
343
344
345
346
347
348
349
350
351
352
353
354
355
356
357
358
359
360
361
362
363
364
365
366
367
368
369
370
371
372
373
374
375
376
377
378
379
380
381
382
383
384
385
386
387
388
389
390
391
392
393
394
395
396
397
398
399
400
401
402
403
404
405
406
407
408
409
410
411
412
413
414
415
416
417
418
419
420
421
422
423
424
425
426
427
428
429
430
431
432
433
434
435
436
437
438
439
440
441
442
443
444
445
446
447
448
449
450
451
452
453
454
455
456
457
458
459
460
461
462
463
464
465
466
467
468
469
470
471
472
473
474
475
476
477
478
479
480
481
482
483
484
485
486
487
488
489
490
491
492
493
494
495
496
497
498
499
500
501
502
503
504
505
506
507
508
509
510
511
512
513
514
515
516
517
518
519
520
521
522
523
524
525
526
527
528
529
530
531
532
533
534
535
536
537
538
539
540
541
542
543
544
545
546
547
548
549
550
551
552
553
554
555
556
557
558
559
560
561
562
563
564
565
566
567
568
569
570
571
572
573
574
575
576
577
578
579
580
581
582
583
584
585
586
587
588
589
590
591
592
593
594
595
596
597
598
599
600
601
602
603
604
605
606
607
608
609
610
611
612
613
614
615
616
617
618
619
620
621
622
623
624
625
626
627
628
629
630
631
632
633
634
635
636
637
638
639
640
641
642
643
644
645
646
647
648
649
650
651
652
653
654
655
656
657
658
659
660
661
662
663
664
665
666
667
668
669
670
671
672
673
674
675
676
677
678
679
680
681
682
683
684
685
686
687
688
689
690
691
692
693
694
695
696
697
698
699
700
701
702
703
704
705
706
707
708
709
710
711
712
713
714
715
716
717
718
719
720
721
722
723
724
725
726
727
728
729
730
731
732
733
734
735
736
737
738
739
740
741
742
743
744
745
746
747
748
749
750
751
752
753
754
755
756
757
758
759
760
761
762
763
764
765
766
767
768
769
770
771
772
773
774
775
776
777
778
779
780
781
782
783
784
785
786
787
788
789
790
791
792
793
794
795
796
797
798
799
800
801
802
803
804
805
806
807
808
809
810
811
812
813
814
815
816
817
818
819
820
821
822
823
824
825
826
827
828
829
830
831
832
833
834
835
836
837
838
839
840
841
842
843
844
845
846
847
848
849
850
851
852
853
854
855
856
857
858
859
860
861
862
863
864
865
866
867
868
869
870
871
872
873
874
875
876
877
878
879
880
881
882
883
884
885
886
887
888
889
890
891
892
893
894
895
896
897
898
899
900
901
902
903
904
905
906
907
908
909
910
911
912
913
914
915
916
917
918
919
920
921
922
923
924
925
926
927
928
929
930
931
932
933
934
935
936
937
938
939
940
941
942
943
944
945
946
947
948
949
950
951
952
953
954
955
956
957
958
959
960
961
962
963
964
965
966
967
968
969
970
971
972
973
974
975
976
977
978
979
980
981
982
983
984
985
986
987
988
989
990
991
992
993
994
995
996
997
998
999
1000
(* Content-type: application/vnd.wolfram.mathematica *)
(*** Wolfram Notebook File ***)
(* http://www.wolfram.com/nb *)
(* CreatedBy='Mathematica 11.3' *)
(*CacheID: 234*)
(* Internal cache information:
NotebookFileLineBreakTest
NotebookFileLineBreakTest
NotebookDataPosition[ 158, 7]
NotebookDataLength[ 61236, 1395]
NotebookOptionsPosition[ 56917, 1320]
NotebookOutlinePosition[ 57251, 1335]
CellTagsIndexPosition[ 57208, 1332]
WindowFrame->Normal*)
(* Beginning of Notebook Content *)
Notebook[{
Cell[CellGroupData[{
Cell["StateMint Example", "Title",
CellChangeTimes->{{3.749517748387988*^9, 3.7495177578584023`*^9}, {
3.752074679813195*^9, 3.7520746861317587`*^9}, {3.7531991258476343`*^9,
3.753199126839426*^9}},ExpressionUUID->"d3849ac4-c012-46c7-b86b-\
5e1ba9703926"],
Cell["\<\
This notebook will focus on how to interact with StateMint in a Mathematica \
notebook.\
\>", "Text",
CellChangeTimes->{{3.749517768393689*^9, 3.74951781569664*^9}, {
3.752074695638905*^9, 3.752074729204295*^9}, {3.753199100036727*^9,
3.753199101053008*^9}},ExpressionUUID->"dee37b68-b003-47bb-83aa-\
34f50c8d94d0"],
Cell["Before we begin we need to import the StateMint package.", "Text",
CellChangeTimes->{{3.749517831842804*^9, 3.74951784732657*^9}, {
3.752074736641693*^9, 3.752074738766007*^9}, {3.753199097196782*^9,
3.753199098364512*^9}},ExpressionUUID->"736874a7-447e-420a-8fe9-\
247fc95dd459"],
Cell[BoxData[{
RowBox[{
RowBox[{"SetDirectory", "[",
RowBox[{"NotebookDirectory", "[", "]"}], "]"}],
";"}], "\[IndentingNewLine]",
RowBox[{"<<", "StateMint`"}]}], "Input",
CellChangeTimes->{{3.749517850744534*^9, 3.749517904428472*^9}, {
3.752074859105585*^9, 3.752074860517879*^9}, 3.752076465465398*^9, {
3.753199093907756*^9,
3.753199094921145*^9}},ExpressionUUID->"4ae703ef-2121-4e4a-857f-\
262964453034"],
Cell[TextData[{
"Now that StateMint has been imported we can start working on a problem. In \
this case we will be working on the problem set forth in the ",
ButtonBox["tutorial",
BaseStyle->"Hyperlink",
ButtonData->{
URL["https://github.com/CameronDevine/StateMint/blob/master/tutorial.md"],
None},
ButtonNote->
"https://github.com/CameronDevine/StateMint/blob/master/tutorial.md"],
". This problem consists of a motor powering a pump through a flexible \
shaft. This pump pushes water through an elbow with a known resistance and \
out into the atmosphere. A diagram of the physical system can be seen below."
}], "Text",
CellChangeTimes->{{3.7495179168731947`*^9, 3.7495179367746696`*^9}, {
3.749517973149335*^9, 3.749518010190508*^9}, {3.752074918533415*^9,
3.7520749200987577`*^9}, {3.752074980552095*^9, 3.7520749808721447`*^9}, {
3.7531991046505327`*^9,
3.753199106022131*^9}},ExpressionUUID->"9a3a49fd-e894-426a-9ec8-\
6ca841f06f7b"],
Cell[BoxData[
GraphicsBox[
TagBox[RasterBox[CompressedData["
1:eJztnQd4FEUbx/936RUCofcivfcOAtIEETAKIlKEAEpRQEBEqUqQIkVEVNCg
ovDRpEgJSBNQQZEmKFV6DS2EkDbftrvba3t3e427vL/nIfvuztzOJNx/d8o7
75TqN7xrvBbAyFDuR9e+Y54eMaLvO91ycydxw0YOHjhsQP/2w0YNGDhgRIN+
fLaz3L9eAQAjCIIgCIIgCIIgCIIgCIIgCIIgCIIgCN/is9oSjZ6f+q+3K0MQ
/kTarW3AR1eTbx79rmXA0HRvV4cg/Ik9CEoRjKxWGOLluhCEXzEFDSRrPoLu
eLUqBOFftMIYyZoJnPRqVQjCr0gLw0+S2RUxaV6tC0H4FbsQcE+0roRgpnfr
QhB+xSTUEY30Dhie7d26EIRf8TRG8oesn+tX3ubtuhCEP/EoFC3j43s3z//s
rixv18UC/V6Ie6lPXFzX8d6uCEE4zA4E3ucON4YETPF2VSxx6c94aD76/fht
b1eEIBxmAuqLRmOs9m5NrDAKVbxdBYJQRXPdLNfreNG7NbFCHQzzdhUIQg1c
l0ua5WqOrt6timXuBjyhb1SCsMF2scvFWGYU+nu5LhZZD+1Nb9eBINTwHuqJ
xhHgiRyRG4ka3q4CQaiiKd4WjVXAAsbuf+zd6phTC296uwoEoYZLQVgpWnuA
7xn7cJ1362MG1+X60dt1IAiHud+yVgSQv44wXJj9DMawNa96u06m/Ahtsrfr
QBCOcyM5OfNBcvJd4ST9m769F2R6uUZmvIVa3q4CQfgnNUTPR4IgXEyyFut1
9pxj3qwJQfgZaxCgCzyQ1eJJdCwmCF9luG6xGWOr3/VmRQjC36ium4ZjqRX+
8WpNCMK/uKXFRtF6/GJD71aFIPyLr6X4Hlk76mGhl+tCEP7Dzq7NQqHtEhfX
tUl+IJzmkgnCVTw4KOect6tDEARBEARBEARBEARBEARBEARBEARBEARhhXdK
d1NIfab0B9YTb5cvvdHl9SEIf+F2BHDQauo2IE+K1dQPgGbuqBNB+AVTAPSw
mtqOS51vLTGtEJf6q1tqRRC+DyeQhgj8z0rqMQ0aoIy1eE9LEFYd3d1VM4Lw
cRYj/FpZjLCS2hdVLwRhlZXUahi8EgFn3FU1guBJXvhK53Ydpl3xdj0cJpsT
CPsEUXcspl4PxVLWQxcq3pSfoDmRWRZvubF6RI7n1oCikzldXX4jYp63q+Io
P0F7gj2MxQyLqe+iyGN2ENhnMbUVnmdsvjVhEoQLWJW3zXXRmgiFcesnEkEg
vIjSLSRyopvOHVrA4lj9EQ32cHny4iO31tCM7OUD3j6kO1kZP0o/1Llm4Ijf
dfaPg96ikRY/YI72FV20vMxyAX94tS6OclgQCN/8+9ZCqtRcXI+AUxZSe4kh
zcbxbzcPkv4cgEAxWkZGN84OEJsLmS9xtnaWYGe9wtvTPVktwh3MRyvDl2sy
4rxYFcd5BXWFYz9UyzZLzCorbs2TXQlDzD96KRjL+eM1y8J0GyOgfbketLt5
eyw03RtCs52334fmxSbQbObtKcALTaHZ4Ml6Ea7n58BI2TDGJkR69DHuJJxA
VgjGCQ22m6Wu0o0Efo7wW2apo1EyQzD6WBKm2zgdyPUP05sJey+cD+aa4Rmt
UIUr/1IIJrDMdqjA2dfCMI5lPouyFH/X3SwpGsNTLE+xGv23uPrmqUUxRXb6
O3DS1UW4kbd1AmHt0MEstZFue/G0guZ9yfu5MUe0jloSptsYhQqZQpn7+BZp
Ga76JzXYyXd3S3CPtdNabOMdR4qkMXYuULeVM+FGLgOfcIcH+waijYujb01F
TKrsdIu1kbUnEoNAWBI0x01SuQfFfsmciAKPTFI/RvQ9yWxrQZhuoxyEblUj
jGWsCj7k7Rb8bEAtTOTtZzCUsXri7rIdMNBz9cqxHNDvjtYPz7r0zmm5MUB+
/iXwt0sLcCuzkfuBzq5ptpf4C2iqM29HYIlxYmZpjNHZW82F6TauAsJGJpPQ
mCVrIAwfJaAOe6AVnwQzUYOlBgijNGweKnmqWjmYGfoNqZcBZ1155/XAz/Lz
sdDes5b3iYMTyFj9yVKEXDVKPRuItfqTQUK/RsYPCLqgP6lhJky3sQMhgjPW
WuRh+xCYxttcRzf7ILQPeTsJYVl/QSP8J+xCcIan6pVz6YDhkvUZcMCVdx6A
kDT5eSPdJofrFwmsPuzK0lwMJ5CL+pP0YnjPKHUonjKMB/yrxWaj1PqQ7d76
takw3cf3KCUc/wQer0ZhwT4O3N+AWMH+B0jejNyCfQ645qFq5VwyorFaMnsg
90NX3rolGshPrwfifdE68DkqbU7auLBj9d2uLM8RFsy1QW2jyeHpiJX/aZIj
8ans9Dk8I7/3buCQ4exxYRNhfj/LVtmO8F7fcfo7L0YV4fgvJ6Jv8ZRg/wdc
XYESgs11ry+uRhHBvgno/R7T+vad7kCR8n40ocRv0NwUrcuhmOvSWz+FXvLT
+QjRjcv/InawWYdwb212GNqmiw0aySeC7kQZqWka8sq1tgv4S3baGa2ZUWYj
YWYUDG9qq2wHqM01AXUsQxnh+BeQtgpFBfsEcG898gv2KeD2Jik/Lzrd5y4C
be0vMeS643/uHMp06WnHUlsHTHTtrcsbua1mlOHHqkSmYIdwHIrZri3SbkL/
dSz/myhtWFuSXtTkZWTUDvzHpJVo8pr7AkGudJVfhwJ6+yfkFY47Eca2I1Kw
9yMgew+CBfsANJn8OW8fBvRjnCeBu/aXmJ/EZS/t8AZ/SPmhYrM/XXzrVkaD
hV+gqH44o2WI2LSoj10uLtNeHBXX+UCs0Z+YdaOMRjAGmo5vDJN30B6XwuuO
Fa3MLoTq7b0IFEpei4LsIDTCgMUmxLAjgPD3TkIU+5vrg5l+7ldoHZhRJnHZ
S0YUWsbHv1onuJXxvoMHZjkf/2GK2BXZ+zX/xbucO1Q3McTSwsRx7JVGc8we
xVFxsZfQRG+bDQDKx97NR+bPyYcW5yH0InMhZ7h+lM7mmnrn+ONM1GHXpBn7
+ajKbgNHeHshKrK7GggP0cVSE5LnW6kJaR8kLnvZD+1t7pC5IKCH/OH1edVf
ijjtFHchJPgW1/RbsHvAf+x+vZDV+oSdGJGcfHzjoKd3OFuEahwW1wHDrPEW
86kr2ayxhTll2aRYahGMdLBkZbJzYZ3+JFYcnOqJeMYK4Xve7ofejJVAIm/H
oydjZbCYt4fIHD1H4jkHSiRx2cs0VBeNOHxsuJoeO3VnkPPuMVMwnG1Yyj1d
373aNI/MC2gC3kxIeKtkRy+udXJYXKyZ/tvYxnyy3eDOYckbSubOMR2RLv5u
NsMkvd1S7AtWwgLG2mM0b9fku7WdRSfjuvwCmG6io3FjTDX9mJ2QuOylrW6W
a4x8iGsnNjAX+NhmD8OcD7jmZmqZ/F3k65Cbh/M3v1ugpfNFqMVxcf2IgNOC
cVRjPDcuoPfUXWTJj7exzhHxQT7R+ciFjIZhs+Nxgh//ecE5YxKqcfZlLe9y
Ng0VOPt6AHbzTgNlOft2EO9nKHIvFI70AUhcdpIeoeunt5PP1cyAi9bkb2qa
p/1HQ+tG8vM+2culi6khrYRj7UDXFKIGx8WVXRHDBMOio7tujQmXy8IKlNU6
YU5EblfvnvorNOd19kGAK+cjlOIqeBQ4wdgcFMkSZo+PMrYABTMZOyt0ur5A
rN5BIxG50yzf2yIkLjvZC634nM2ONUQqOnewY/TBPyytvlXBnYO7T1/Q8N+o
H3SjF9swmT+kBhWw/jF347i42EJE8H8rK0u0pNWRltdOZpUVGwh3crt+NXZ2
CSToT0rzjb2aeJu3K/Bz9vXFgqvxvlyNxXHK2nyv72nZmEw7+USCbUhcdvKB
0HZgwvzim7qLn42JKj/mPfPGjWpO8LPTOwvo3objIThm7Edt7pnqJQc3FeJ6
VADTuMM7lhcXH9EIv5WVVf+fIIIfNxqD2PsOl2uL8Sigj0w6kStgAzRHeftD
xNzdCo3gvjsT0bd3AL/x9lxE3tgL0X2X54CDq2JIXHbSWmrqsK3gI0I8EH1V
UwNGu7SUkwhsU1+/xCGjTpDQCtnL9fIeP+elJXsqxMXeR+HHfFgMy/FqWqEz
sxqvRgy4cSPSHZPmt6MNry7u1TipjqTv+3nxXkN0EuyUAhjbDO0EO7UwRrSG
ob/bXjbJYA8kLvs4HYovRItrlS/kHnY/Cie/uXhtenokgEqSE0C/Fq1bt+JH
JrMG55nf63fFD7oPNeK6HoZENs9aQCc+lpr1SGtCqKhhKOQOv7z3EHZUZ0+G
Btojoj2NszVS2JJZnK2L/zuft/UvrsVwcDknicsO0quXzBsTk6+Y2Ncaghfv
fRkvpnyOI64t6quYyB4WfLCTj3nNB1SNuNgAVM0sYzVGaDUMvmg1Rigf5PBS
GBaqKNUmD6ug4m3JflQJ+om0tGrQe5yl1wIGSXZGfaCf7sN/hKOvY8WRuBxn
34Sh/5NGwd4I9qVYF6pQJa4TGgxC4DkrqYsR3tt6dOt+qDEApdzzd/0nChWl
EcM1YUAZ6XdbGw6UlhxvtkYBpU6IdlI0UFjyNN6eC5UdXApB4nKKxjVs5/Fx
VImLdeSaty9ZS0wrqLQvwzGuiYav1RRqB1uiEJNwmWUd6KVBtWjk+vASy/qj
txZVciN6ygWWfaifFpXyIGrSfyz7r/4BqBCLiPfPMXZkUCDKn3OwMBKXM2RH
9/ZYWRk/eWdEQ524dihuXDIViHlgNbUdUM5tQ6OHinI1y8W9tdDxweFi3CGa
t9vfP1aKO0RxbzA8c/fv0vx13n46+d+ndNebOjwqTOJyhtP43GNlPUKR2Q4s
d3AZ6sTFaittuXUrAuOspyYBy62nOsv99/kdi1CPL+LBhMK8Xed7rpWfMrkI
b9f6hrMfTuUliBqJ3PMsNaEEb1db4vizjcTlDKtcG0pDkUdAWHBfzy+ZVCmu
1GSlIZj7yUoz7+fOuPUlnX1yy66renvrTv2imH+27ryss/+V2aeSdqjyzydx
qeZsInuvrtnVG6sm94+zwo/OFPeIf34GB9Zc4eHJZJXiIkhcTtA1PKXqGuNL
mUubBlTvM3nWIsv8ZflG9iGIC9CExbx/w5n7OAqJSzUkLtUsa9fLJEb/+srl
F9x0V3GSuDjCAzt7MGQoiUs1JC71ZBjPxaR0K7TQjU02g7iAwOBKi0zXGboL
EpdqSFyu4lLtZ+xdIzGpuQqawghN3sku8sa3AYlLNSQuF3G12FC7X1t7vlHB
EiNtBYeMPG+7IFdA4lINics1PGrQ290lGJQVEFxukUuDkipB4lINics1xDdx
t5OhXlyhgc8lubksOSQu1ZC4XMKhCP0s47mZXRrWtkyiM0VI4grL9fYF25ld
CIlLNSQul/DMO5Jxf2B4pzmrkizjVC/pkYdHCXWQuFRD4nIF+2OkZYGXqrQ9
7aYyHiE0qMdB2/lcDYlLNSQuVzBYirv8sFZ/a6uUnOZRkQ/dNkGtBIlLNSQu
V1BCimw3oZH7ZpGz3SZbZUhcqiFxuYCLgaIH+MOo/TZy+iAkLtX4k7hOb/v2
40veKDipnHhcW8EbpbsZEpdq/ElcAxtpPLi+SsY3LcTjpH7K+XwSEpdq/Elc
LDOsuFfK/ayjeHz9HeV8PgmJSzV+Ja6DxrufeowF0r4y8e97pXj3QuJSjV+J
aya+9Eq5JC7CEn4lro6wENrfA5C4CEv4k7gycxfyTsEkLsIS/iSuA8LePhlr
Ezy4CF6AxEVYwp/ENYMPMH7mlRXf5POwCx6Ji7CEoriuXE92CzfdM9X7LP5m
Jwc8uBfh4q1HbELiIixhVVz73m0WC7cR0+jtHa4O/JiZK3/2id5p7HHdttZD
JbsFEhdhCcviyl5RMU/frw8W3+qeQo+GHvo2vmCpJa51Rf0dL/zRRxbm9bLH
Vj6RuAhLWBTXvnrFv+JXrZfa5Z5C/wnjfmQsL1fFpcvVP0KFfjKn9FuRn7ny
7kr4rbiuHdy+LujTpAPnzfcNJ2xjSVyLIqeLD323ioux9E+jE2zkdIQOmPB8
4RH6rSiWw2N7MfqjuC5/279ebn0zPrzGy4voDeYo5uLKHl7gF8l0s7gYO1zy
FZctgMqILsa9cyOr6wLFDI70WFB1vxPXucmVJFVF5C9dKFqyi41wKh53zsNc
XFNK/6cz3S4udqP6UKWMjvAbenM/X+D3YV/HziQkFC+eMNMzYTP9TFzZ61to
OCnle3HWhlPi8yn7v6QFfcWNdBZ76E/qF5iJa3W0YZdf94uLXS7sqp7RdGE3
wvY4yx4/z5K3f4vBSb/Y/JBr8CdxZa+ozomozMQjZt2sM3NqcynFPyF52Yup
uG7kWmc48YC42N7wM665aXvwIcdaBmeyNV8xvsvlOT8NPxLXseaAttsvVgYw
jg0OBSr+7Nkq+S6m4hrWRXbiCXGx+B4uuWdG1FP8YUJ01qOO/Hj865Gee8L6
jbjS3wmCpuffCjmujAgBet7xWI18GhNxnQ0/ITvziLguR7rEWelu3g/4w8Oe
cd0O8Ealdq64q334i7guNAaq22pLn27PtQ097b3pm5iIa3RP+ZlHxMXe7OOG
Im5opjHmqQesn4hrSwwCp9kxxLosF4I9Nofoy5iIq/gW+ZlnxPVntBsacOu4
Llf6866/r0X8Q1yrQlB0j105zzcAxri5Nv6Asbj+CUuTn3pGXNkF3TCqtxi3
2bjdrr+vRfxCXJ9oUd/eqKNp3YHB5LVhC2Nxfd/QKNEz4mKd57q+iJQOfQZ7
bFTLH8T1rRat7ltK2GPJwTp7JDDEzTXyfYzFNbW3UaKHxDXKZRPJclztdK+A
H4hrfSA6Wm6dR1t2e3obmOHOCvkBFzoaNZ1GvmWU6iFxfeidmE0qST2UtHLF
up/PysTr++I6HoHGqZaTwo5bvJzdB5r1bqyR77MjuknMFNn5G2ONkj0krjnd
3FOM67md2Lu0JqZi7Wa1SweH1BudJK2Z8XlxPayC8tY2dLYiLpbRBnk9u1OY
j9FsHjsR+NFcPd2No1qSuIzY1jW0wfgkqR2dcWrZgGKFRgl+mD4vrgEIO2wt
zZq42I0iaOylzSF8gpor2R1thy56OpG4rJJUP9/7vKfWNcNuxFnbu4e9dtn3
xbVDg0VWE62Ki+0OxHy31Mc/WFBwXD35ZNAbJC4rXO0eMytFsKqtlV8/9Ur0
7CwfF1dGVbS2nmpdXOwtRF9xQ338hbWvz5XPbJG4rLC1QM9rktlzinHSrgqt
pvm2uGYh5B/rqQriul8E7vCu8VNIXJaZFb1Ubye8aJKY8mquTqLlm+JKLQCl
DSQUxMWWIdBdu9T6H74krnOrEzfq1vT/tyZxo8674MKaxA03JPvi2sQNuom8
y2sT119VUxBjYwseMpxsrGiUtud1xt7dIdq+Ka55iLDgmbEqTiKgrWTMNM+U
VRn93V9BP8F3xLWGX7gHbSf+W7++nmB34L3rNzYQ7HZ8/IwtjXhb04bf7nFb
U8FupcbTalIx+eP5YoB8QuhBWb5PL8WZeqLEdfFicvKlM2ekqp09c/X6mRuW
8mUWxygLl48tkgiaKBk7LeRaiuDLrquyf+Mr4krvzSklf7loIGBexgB+TXq5
XJykZmUO4uxYwZ6eNYyz85bLzUlqatYIwY7h7AkO13ZJXuOGUZ4/ZSe9n5G5
2D1R4urdWAuUbS9W9l6TQKD6h5bybUGgoj6UmoUsoyg+Ul/FnIWPiCulDfDS
X9xDd0tdoDzQ7Q/OTuLeWeWA5w9w9vbGgt3pN67lsqOpYLffn82ydrcAXnPQ
L+pIpMkju2miwV4dKx8ue6LExdjzqGpQ/tj8lt48HD3RQfEuiuJiY1DZ8Yrl
THxEXK8i4BPRyow3vIwyB8vsNwx29mjutTVG/JplTwBMRvts8LjyZJMrrxta
UZfzLpenPGHiGoOSeju7uRU35pRwfK94F2VxHQMoJJR9+Ia4vgaW6OzvgFDd
/+4KIFi3qnkNEPSrZG8AAvV9rbEIsPIIt8y0qhnXXr8rv7Kwrc7K7tDXKO8T
Jq6FCNT74ia+bSXTJoQ+tJIkoiwuVgXTHa5YzsRBcU1v3SUuLq5Ta2Fd/f/4
k3ZWvWjkOCeuOzFCEDWB+7EIQiPRTinA2XXEN1RqYc6uIdppxTi7iq4xmNkM
FRxw2rkeuY9d71Bih+zS3sI6a26pe0aZnzBxbYZ+//VbdayFzx+Np5XvYkNc
w+HBeAo+jYPiunsmERh++ILwZX0wFGXnHbHre+ucuCYgVr/waCpitmmwWbCn
I9fPWojBq2YjcqcWqwR7PsJ3BeAH3UdOB+Eb+wsbK0xhJUYNM8y139dIo25/
m/bGnjBxnQS2SWZvqw7sdWHa7DXBhrjWIuKxYgZCwuFm4QZob0vm3SYJ9oa4
dUpc93IZBqhSYjCVdUZj3n4UiwksDnV5O60AxrGXUYO30wvjbdZb1u/uj/J2
L6JNid4vHM82rmyY6pLeY+l1x5vkNhbXet1MUVz3oV9Z8zl3J4+0+EK0dlqN
sJURDBtrSm2I6xZw1OGa5UgcFtco8QvMcayzXS1CAafE9S1ypejs5Yi6z34F
eIfaNQi7w/4E+PBV6xF6mx0V/9u3IOgG/xD/Q/eh0xo+Kq99LNX9ehkJUQm6
lmXHecJhTC3TR7axuG4e+B/wxkGOTROL5P7C7l/QdRTmHjE8afWtDrafAmzM
rBe8qJyeT2ogEDZwWFy18aZofOdILDunxPUcXtXb3cA/kssKfeoe4G9TSWjl
9EJn7mc18F/219BeqKhhqVoDi7OmFmk7W2/+Vq6ltHxp3AD+556oE6a5TZuF
/wJSaPC7LeGGWAa2aAzxjTX+c6tZNiDaxk3SbKQ3gcXpM8IUR8V1NwBr+GPa
yARHwpWYiGtg0Zfsg3eNSA3Fj7oPpoULT82xaMAHBcUyxnfIavI7jwvdqimo
ylh2LL5kfIesvL68WSjDHw51jbNF12CZb0bqsFziwozvG/C/eslPzX4vU3F9
gVI682Z0sIJzrJvoifr84WRz61N7n6CWk4X0xQAn75BDcFRc66Dl3fsuvLDX
oWJMxNUl1wx7mAp+bO5XaPUDX39Aw1/6GSHpfCOQ95Dbi8BUdgLgPQp/Q0AK
OyO+PQ5Bo3+1cln5WjeJS7DFgIJGFd1cKI7vYR6L4L6sr7Q1f5yYiutlGIbq
+3ghoMt45Od+ZrdV6DV9iFZOFvIWXnLyDjkER8U1AtW5nz/Vsei0Zh1TcRWx
61P3BHEtRAX9lS8ghLC+o8FhthTFeDtFi9/ZdxDumBqIfWwl8vL242Doh/Yy
QoVRtJZf2yxyhokH/I3OBX/iPh+yL3NVrHFXJev06utm4ioMQxFzxbp6lC8B
7kG06D2FLO/A2eCOE2x4eBASjoqrJoaz7ARtiNKu5RZwRlyD0F1/ZQjihGNJ
JLKRENd9lMfnbDSeFezK+JS9Jy0ErIk5+s/VBe/hbYe4XjMbW0+Mjr/SV4O8
ESv0l25sm9O/XkRorV9NxfUPcE5/sgRaj6+J384P6Vyrp7R97TD0VEi1h5lo
6uQdcggOiutOAFbf7Ls21NEurTPi6iyNgPF0w2jh+DQSuCbYcMFuh0msl9QE
64j3WH+pS9BVFiP2RYxkdomrlXmWc03DhZ3gPubsx8cSx3QsjUIdxyQe5Pv9
JuJaJHM/YjOEt4hn4ZrEP7Lu25SyjHC6UfchWjp5hxyCg+JaC82mXle4Lm0R
xwJSOyOu5jI37FaSqrtiLL9pq2D3wFusE8SWUC8MZS9IAnwNA/WfG4jXmF3i
qrPG/FqiuM1iyKRelQJiGsfP2WPQjIm4usvX6fZHlM3iXE1GIOZsMtQh87u3
3lho8j/lfKPO+YZlDsFBcb2J4HcyGD+5tEI5owlqxdWySZMmEbrBf8aPr4sj
dv0wiDUS3iWMDUZf1hTiwr6h6MVaQ9j9hGs2Gh7RzyKWu1Mh2+KquNX8Wn9p
E9NnZyddM0kyEVchgwskp1NUs1mcyymBAXX1/eHsjp9msPmjjXPMcLpRNwQ+
FXbSezgorhpS86upg/9DasW1dNmyZRPHGfYRryqMsvOC6sN9eRcI9pvowepD
nOYdxfXJmklCGydMfYlsHvchd6fqtsVV9Sf52cPfv3yzVT6NJC4LfgnG4jph
cO1jLDnYG1sXtECAYX3MZn7os8Ebxjm+Qlkny+gme9oRCjgmrttabBSM/wEH
HCnGmWahnAbSIMUrnMibS97Zr3F9rJZSc3Eg9xZrLzUXh5r13O1oFjbSvZEv
r0vQtQJ/CxK0ldvC0nhjcS1EccPJEsAl25A5Rh+0MEwYLMOzX566YzLltRcB
tmaJbVAZC527QU7BMXGtRqD4dc8sKZvQOS04JZxSWpLoKnF1wQjhyDsJdIf4
TG7Fiakn10zkaYvxXJNRrFonmDSI7BHXc/PZjaRZfeuEhdfpN1sXDnQor618
LQtuNMtuLK6X5O2lNmhrmtsDTA07ZTh52J6rdx+TIcvbgNJGkrbJCrXlnEiI
OCauYbrFHmwadKPx2e/MqPwZP9SxxOrHXCeud9GCP2Tkwlo2WaxMViyWsw9Q
T6hLAXzHZkrej8Xwlcmn7RDXqKFs4FPdJq78V/aouFtiTP+S9W/yg/KmC6GM
xMWV/qX+5JAm0BtrCq/JIhJkPGaXV3c0U0I+wbFFPccBiqJhF46Jq5p+VPxW
iDRswBJ/YK3aMHY+v6nHuBxXiWsrQni3i128n8hOBPJuF79Cc4X9ggBe6weB
i+x3aC5x9hHRu1eOHeJabGGt08vtstln3K/IzjWr+IdxkpG4uG+d3ncqS9fz
8yINm3M/0gIPmVx+HvFO3Xa+6EtG2MQhcZ3VQL/IvTYKiktPemVeD5zIHd93
5M0VWNYeSpuJKy2v8J3tzr/A0gsIHa3eaMg9pAsL7ruv8etPsoqD1/lgwZfE
iEFFbBZZPMhsAnZlvqtcT6UQb2YmhE8wamUZiWsBiurtCejt9a3iiq3kfizv
anp5rpMjGl3JtdBO7BfXf51axgAlWwsvqMnVudZ8wQ6iQ9AnGn5Dp/EnFYox
EdfEZj/bh5mLwzvIfYpt1AjuuxMQfZIlaYVWzhREHmc7AsCH8pyO8KNsT6Cs
jSaRYkeJhUxbUZfy8KsxdeslD1dtKI+6ZiSuOP0ISsYYvOXB3cKssHXw2MmD
J5o9LI4Bljfgso+03E42K3MO9osr8/SZM8mXz5wRpnqunTlz/eqZ0+LDuRnv
M84Uwxy7LOLu/WLI+0Iw+EYaSymFPHGh4vBYalnExIWhCf+VflQeuV4MR31V
zkdDTNpMWS3FCyUl0aWOySXbwkAurtS8EBd6ZCyvVWSlmrJdTqbFJcPlpeFU
daxEqKe2dfd1XBCgJjuCd0PaO1spj+vCWR8twr0z64rLfP8uxtk1xRHyf0pw
djVxkOVMac6urC7Y7u+5jbeEm1lGjDDQSb86a2uRbrq4vwZxpQ+KKwc8HR8f
3+uZPPXmWdwO9UlhMko60WTtDNPo3oQVXCCurHCu6/P4ZSVnUVfGir+3YOQy
XXSBBwtHfqPz7kn5bMRSnf3w8xFfq43zUMtok5xjkeKqf/auIYrznR4FNkim
4c21P0li23ErGzY+MZzXYovqD18IBm0vaSeuCK32QeW/T8YfU8ziQxsxrCwu
k2VatYmS9UN9WZ4VMfFi4IEnLECNnXREc0uXe5U2pp6lCClDUdrewCk5HpfE
LTw4J9FGOBYfEldWzQ8MJyPq6F6GxyPkIxRnG9cXWla+Ka5fgd0WLp/dlmSE
pTA018IUds0jjPGNoKAeZX+kfn5sV5R+pX5GiNHWOZniDLFviou1Ql2Vg5n9
UdRJ56kcBInLnGH1pIbhneKyOC81LCxG8VVx/a41n6ewi/1aM7cXwiokLnPS
akvBL3o8J7vaa5KFrD4qLjYIsQ4GahBIr4kmXp8b9x1IXBY4X3gaf/g2n3z5
1kcvyLNIXzFfFVdyfrRT0TAciaAjrq+M30LissSRPNMZu5zbaMx5U3nZyV1p
AMdXxcV+0mCawx/aoKE9GByBxGWRw4WGZ2QlGV26FGBwiT9fSmoj+qy42GgE
rnPwI0dy41lqFCpiPEtB4rLM2RqNL5hcyqtf+7g2T3NpkbPviiu9OcL2OPSJ
80VQ+pbtbDmY620CysvjeZK4rPBocPRs4+dQ86/E4+WXcv9vvhSixXfFxe5U
Q4zdEfQ5LpRDfmccfnMAvV68NKfQ1Rt6BpG4rLG7Srmv5JGTZnzH/7wyMurV
a8wPxMWulELkZrtzHy+GKIciO+RAKm1lmVoYaEbiskrGN5UKjTggH1R7sLpL
aFd+re+8LuIFXxYXu1ARwfZuxZKUB7GOvOdyJAPaH5lQKjVNz+skLgWyfu4T
k+f58d8lHfxz36aFwxsHVZwkem/4hbhYcmPghbu287HMCQEoobRej+C50yOm
oTy2w1DjIC4eEtfsOPcU4wYyD8zq3/SpQjElKncYkajfuGqutL7Xt8XFHvYF
ntpuM9vRJkBLdSt4chj5jcK8v2Mc1c5D4prUzz3FeAzdq9fHxcVYYgTQ45Ji
lnsjgxAw0eNh730SY3HNNY5M7CFxxb/rnmI8ht+Ii51qB4S+ft5q+u0JMUDN
fR6skS9jLK4txtvceEhcjR3YAPyJ5GPJF8r3xcXYqhJAUNx6S6u0sn+JjwJi
5tFry06MxXUv0CgGnWfElRJs/UnpG+g6jf4gLpa++CkA+V/77orR5Ts/DufD
JeSfZhqCi7CKsbhYC6NV7Z4R1w+V3VOK5/ArcTGWuapzMD8xU7LNkNlfrNi8
YsncUZ3KBXAXNM2+NI2DSihgIq4lVeQTOZ4RV4sE95TiOWZJQVr8RFwctz6L
i4Uxke1n+HoLw9OYiCv9KXkYWo+Ia2M+n29ozJR2J/IfcXFkH/5u/IsNq5bS
FK5U7/nRi/dTrAyHMREXW15cFq3IE+LKrDbfekYfwS/FpSOUXAjVYiqu7Oav
Gk48Ia6RVR3blfJJZIa0TTOJi5BjKi52q7RhW1QPiCsx7ymFjD7CRz3EI4mL
kGMmLvZXtD5arvvF9XXkDvcU4VFIXIQlzMXF/iqhi5frbnFljslr25XNB5j+
sngkcRFyLIiLXW1YeZtguFlcm6pXNd0zyzchcRGWsCQuljEnb+Ol990rrofL
n8433/fHMgQSpM2CSFyEHIviYuzBnJqB9fvHDFnkFt4NGtgkuFLCE735hyNM
e0U8krgIOVbExXHp+/EV2sbpaNQtTkbXxnF20LSz0WmjrjqrU6lxS8968rd0
M3ekP+JNGzHyfRISl2qsi8uIxzDKdybIng8V22t0qvnP3koRTxAkLtWQuAhl
SFyqIXERypC4VEPiIpQhcamGxEUoQ+JSDYmLUIbEpRpjcf1vmBXewAD5aW+t
tYxyouKMTtHXSr5V3vrlCTsgcanGWFxrRllhZN3h8tMh9a1llNM03ui03lAr
+dZ665cn7IDEpRo7m4VEjoXEpRoSF6EMiUs1JC5CGRKXakhchDIkLtWQuAhl
SFyqIXERypC4VEPiIpQhcamGxEUoQ+JSDYmLUIbEpRoSF6EMiUs1JC5CGRKX
akhchDIkLtWQuAhlSFyqIXERypC4VEPiIpQhcamGxEUoQ+JSDYmLUIbEpRoS
F6EMiUs1JC5CGRKXakhchDIkLtWQuAhlSFyqIXERypC4VEPiIpQhcamGxEUo
Q+JSDYmLUIbEpRoSF6EMiUs1JC5CmdA9lwh1xJK4CEVCQaiFxOUH/B9lMR2H
"], {{0, 222}, {861, 0}}, {0, 255},
ColorFunction->GrayLevel],
BoxForm`ImageTag["Byte", ColorSpace -> "Grayscale", Interleaving -> None],
Selectable->False],
DefaultBaseStyle->"ImageGraphics",
ImageSizeRaw->{861, 222},
PlotRange->{{0, 861}, {0, 222}}]], "Input",ExpressionUUID->"049bf644-cb53-\
4e85-b08b-365e0b6afffe"],
Cell["\<\
The following equations for this system were found in the tutorial .\
\>", "Text",
CellChangeTimes->{{3.7495182918191347`*^9, 3.749518296746046*^9}, {
3.752075750946992*^9,
3.7520757511794786`*^9}},ExpressionUUID->"558df971-96b2-4ea3-88fa-\
d91934866a49"],
Cell[BoxData[{
RowBox[{
RowBox[{"inVars", "=",
RowBox[{"{",
RowBox[{"Vs", "[", "t", "]"}], "}"}]}], ";"}], "\n",
RowBox[{
RowBox[{"primaryVars", " ", "=", " ",
RowBox[{"{", "\[IndentingNewLine]",
RowBox[{
RowBox[{"tk", "[", "t", "]"}], ",", "\[IndentingNewLine]",
RowBox[{"vR", "[", "t", "]"}], ",", "\[IndentingNewLine]",
RowBox[{"vL", "[", "t", "]"}], ",", "\[IndentingNewLine]",
RowBox[{"i1", "[", "t", "]"}], ",", "\[IndentingNewLine]",
RowBox[{"w2", "[", "t", "]"}], ",", "\[IndentingNewLine]",
RowBox[{"w3", "[", "t", "]"}], ",", "\[IndentingNewLine]",
RowBox[{"P4", "[", "t", "]"}], ",", "\[IndentingNewLine]",
RowBox[{"QR", "[", "t", "]"}]}], "\[IndentingNewLine]", "}"}]}],
";"}], "\[IndentingNewLine]",
RowBox[{
RowBox[{"elementalEqs", "=",
RowBox[{"{", "\[IndentingNewLine]",
RowBox[{
RowBox[{
RowBox[{
RowBox[{"tk", "'"}], "[", "t", "]"}], "\[Equal]",
RowBox[{"kt", " ",
RowBox[{"wk", "[", "t", "]"}]}]}], ",", "\[IndentingNewLine]",
RowBox[{
RowBox[{"vR", "[", "t", "]"}], "\[Equal]",
RowBox[{"R", " ",
RowBox[{"iR", "[", "t", "]"}]}]}], ",", "\[IndentingNewLine]",
RowBox[{
RowBox[{"vL", "[", "t", "]"}], "\[Equal]",
RowBox[{"L", " ",
RowBox[{
RowBox[{"iL", "'"}], "[", "t", "]"}]}]}], ",", "\[IndentingNewLine]",
RowBox[{
RowBox[{"i1", "[", "t", "]"}], "\[Equal]",
RowBox[{
RowBox[{"-", "Kv"}], " ",
RowBox[{"t2", "[", "t", "]"}]}]}], ",", "\[IndentingNewLine]",
RowBox[{
RowBox[{"w2", "[", "t", "]"}], "\[Equal]",
RowBox[{"Kv", " ",
RowBox[{"v1", "[", "t", "]"}]}]}], ",", "\[IndentingNewLine]",
RowBox[{
RowBox[{"w3", "[", "t", "]"}], "\[Equal]",
RowBox[{
RowBox[{"Q4", "[", "t", "]"}], "/",
RowBox[{"-", "DD"}]}]}], ",", "\[IndentingNewLine]",
RowBox[{
RowBox[{"P4", "[", "t", "]"}], "\[Equal]",
RowBox[{
RowBox[{"t3", "[", "t", "]"}], "/", "DD"}]}], ",",
"\[IndentingNewLine]",
RowBox[{
RowBox[{"QR", "[", "t", "]"}], "\[Equal]",
RowBox[{
RowBox[{"PR", "[", "t", "]"}], "/", "Rf"}]}]}], "\[IndentingNewLine]",
"}"}]}], ";"}], "\n",
RowBox[{
RowBox[{"constraintEqs", "=",
RowBox[{"{", "\[IndentingNewLine]",
RowBox[{
RowBox[{
RowBox[{"iL", "[", "t", "]"}], "\[Equal]",
RowBox[{"i1", "[", "t", "]"}]}], ",", "\[IndentingNewLine]",
RowBox[{
RowBox[{"iR", "[", "t", "]"}], "\[Equal]",
RowBox[{"i1", "[", "t", "]"}]}], ",", "\[IndentingNewLine]",
RowBox[{
RowBox[{"t2", "[", "t", "]"}], "\[Equal]",
RowBox[{"-",
RowBox[{"tk", "[", "t", "]"}]}]}], ",", "\[IndentingNewLine]",
RowBox[{
RowBox[{"t3", "[", "t", "]"}], "\[Equal]",
RowBox[{"tk", "[", "t", "]"}]}], ",", "\[IndentingNewLine]",
RowBox[{
RowBox[{"Q4", "[", "t", "]"}], "\[Equal]",
RowBox[{"QR", "[", "t", "]"}]}], ",", "\[IndentingNewLine]",
RowBox[{
RowBox[{"v1", "[", "t", "]"}], "\[Equal]",
RowBox[{
RowBox[{"Vs", "[", "t", "]"}], "-",
RowBox[{"vR", "[", "t", "]"}], "-",
RowBox[{"vL", "[", "t", "]"}]}]}], ",", "\[IndentingNewLine]",
RowBox[{
RowBox[{"wk", "[", "t", "]"}], "\[Equal]",
RowBox[{
RowBox[{"w2", "[", "t", "]"}], "-",
RowBox[{"w3", "[", "t", "]"}]}]}], ",", "\[IndentingNewLine]",
RowBox[{
RowBox[{"PR", "[", "t", "]"}], "\[Equal]",
RowBox[{"P4", "[", "t", "]"}]}]}], "\[IndentingNewLine]", "}"}]}],
";"}], "\n",
RowBox[{
RowBox[{"outExps", "=",
RowBox[{"{",
RowBox[{"QR", "[", "t", "]"}], "}"}]}], ";"}]}], "Input",
CellChangeTimes->{{3.7495183069606*^9, 3.749518538126053*^9}, {
3.7495186172483473`*^9, 3.749518640291297*^9}, {3.752075776698476*^9,
3.752075926490069*^9}, {3.752076585839171*^9, 3.752076586211344*^9}},
CellLabel->"In[16]:=",ExpressionUUID->"0c27e491-35b1-4cda-acb1-b9a5c08763cd"],
Cell["\<\
Using the following command, the state equations can be found.\
\>", "Text",
CellChangeTimes->{{3.749518799424985*^9, 3.749518810605731*^9}, {
3.7520759769970713`*^9, 3.752075987481391*^9}, {3.752076366531571*^9,
3.752076369261776*^9}},ExpressionUUID->"243d6832-afeb-4a11-b0b2-\
10d3ab515d3d"],
Cell[CellGroupData[{
Cell[BoxData[{
RowBox[{
RowBox[{"stateSolution", "=",
RowBox[{
RowBox[{"stateEquations", "[",
RowBox[{
"inVars", ",", "primaryVars", ",", "elementalEqs", ",",
"constraintEqs"}], "]"}], "//", "Flatten"}]}],
";"}], "\[IndentingNewLine]",
RowBox[{
RowBox[{"stateVars", " ", "=", " ",
RowBox[{"\"\<state variables\>\"", "/.", "stateSolution"}]}],
";"}], "\[IndentingNewLine]",
RowBox[{
RowBox[{"stateEqList", " ", "=", " ",
RowBox[{"\"\<RHS\>\"", "/.", "stateSolution"}]}],
";"}], "\[IndentingNewLine]",
RowBox[{
RowBox[{"\"\<state equations\>\"", "/.", "stateSolution"}], "//",
"TableForm"}]}], "Input",
CellChangeTimes->{{3.7018797195836763`*^9, 3.7018797204767036`*^9}, {
3.701879810803969*^9, 3.701879844544981*^9}, {3.7018942050184517`*^9,
3.701894566874161*^9}, 3.701898020551093*^9, {3.701903282327601*^9,
3.701903283500523*^9}, {3.701911200698804*^9, 3.701911202601316*^9}, {
3.701915682054678*^9, 3.701915720675168*^9}, {3.701915759772274*^9,
3.7019157927706213`*^9}, {3.701915923602686*^9, 3.701915923900078*^9},
3.7019209193048553`*^9, 3.701920950145496*^9, {3.701921306764739*^9,
3.701921308673313*^9}, {3.701921408184596*^9, 3.7019214156809187`*^9},
3.7019220747800493`*^9, {3.701922189786292*^9, 3.701922218088579*^9}, {
3.701922248592429*^9, 3.701922251876824*^9}, 3.701922322852202*^9, {
3.701922383801173*^9, 3.7019224315654707`*^9}, {3.701922463501308*^9,
3.7019224895050097`*^9}, 3.7019225472984467`*^9, {3.701922594846767*^9,
3.701922598092758*^9}, {3.701922662563005*^9, 3.701922735877761*^9}, {
3.7019231046383343`*^9, 3.701923105684675*^9}, {3.7019639715923243`*^9,
3.7019639718841667`*^9}, {3.7019642654040737`*^9, 3.701964352167139*^9}, {
3.701964387110903*^9, 3.7019644147887087`*^9}, {3.702757419429731*^9,
3.7027575039883213`*^9}, {3.7511302149210863`*^9, 3.751130222225973*^9},
3.751131245052647*^9, {3.751758020590618*^9, 3.7517580314940643`*^9}, {
3.751758096034569*^9, 3.751758150178234*^9}, {3.751858834585458*^9,
3.751858867892303*^9}, {3.7518589221835537`*^9, 3.7518590012649117`*^9}, {
3.751859049763453*^9, 3.75185909971*^9}, {3.7518592622618513`*^9,
3.75185926422972*^9}, {3.751859306929682*^9, 3.751859367871738*^9}, {
3.7518594346930857`*^9, 3.751859464138123*^9}, {3.751859792301899*^9,
3.7518598121041403`*^9}, {3.7518598489196587`*^9, 3.751859923477129*^9},
3.75186099185213*^9, 3.751937475686306*^9},
CellLabel->"In[8]:=",
CellID->1449591737,ExpressionUUID->"68f14f0a-68cc-47c8-b1bd-da3eb46ddfb3"],
Cell[BoxData[
TagBox[
TagBox[GridBox[{
{
RowBox[{
RowBox[{
SuperscriptBox["tk", "\[Prime]",
MultilineFunction->None], "[", "t", "]"}], "\[Equal]",
RowBox[{
FractionBox[
RowBox[{
RowBox[{"(",
RowBox[{"kt", "-",
RowBox[{
SuperscriptBox["DD", "2"], " ", "kt", " ",
SuperscriptBox["Kv", "2"], " ", "R", " ", "Rf"}]}], ")"}], " ",
RowBox[{"tk", "[", "t", "]"}]}],
RowBox[{
SuperscriptBox["DD", "2"], " ",
RowBox[{"(",
RowBox[{"1", "+",
RowBox[{"kt", " ",
SuperscriptBox["Kv", "2"], " ", "L"}]}], ")"}], " ", "Rf"}]],
"+",
FractionBox[
RowBox[{"kt", " ", "Kv", " ",
RowBox[{"Vs", "[", "t", "]"}]}],
RowBox[{"1", "+",
RowBox[{"kt", " ",
SuperscriptBox["Kv", "2"], " ", "L"}]}]]}]}]}
},
GridBoxAlignment->{
"Columns" -> {{Left}}, "ColumnsIndexed" -> {}, "Rows" -> {{Baseline}},
"RowsIndexed" -> {}},
GridBoxSpacings->{"Columns" -> {
Offset[0.27999999999999997`], {
Offset[0.5599999999999999]},
Offset[0.27999999999999997`]}, "ColumnsIndexed" -> {}, "Rows" -> {
Offset[0.2], {
Offset[0.4]},
Offset[0.2]}, "RowsIndexed" -> {}}],
Column],
Function[BoxForm`e$,
TableForm[BoxForm`e$]]]], "Output",
CellChangeTimes->{{3.7518593365167*^9, 3.751859368149544*^9},
3.751859469567102*^9, 3.751859641381487*^9, 3.751859929862179*^9,
3.751859981942204*^9, {3.751861010775159*^9, 3.75186103657938*^9},
3.752076373756929*^9, 3.752076484471745*^9},
CellLabel->
"Out[11]//TableForm=",ExpressionUUID->"9255e3e8-baa2-48e7-93e4-\
d02a205005c6"]
}, Open ]],
Cell["In a similar way, the output equations can be found.", "Text",
CellChangeTimes->{{3.7520766271595907`*^9,
3.7520766467332363`*^9}},ExpressionUUID->"fee207ac-a92c-40f9-8955-\
3a0ef126ebcd"],
Cell[CellGroupData[{
Cell[BoxData[{
RowBox[{
RowBox[{"outEqList", "=",
RowBox[{"outEquations", "[",
RowBox[{"inVars", ",", "stateVars", ",", "outExps", ",",
RowBox[{"elementalEqs", "~", "Join", "~", "constraintEqs"}]}], "]"}]}],
";"}], "\[IndentingNewLine]",
RowBox[{
RowBox[{"outEqList", "//",
RowBox[{
RowBox[{"Thread", "[",
RowBox[{
RowBox[{"Table", "[",
RowBox[{
RowBox[{"y", "[", "i", "]"}], ",",
RowBox[{"{",
RowBox[{"i", ",", "1", ",",
RowBox[{"Length", "[", "outExps", "]"}]}], "}"}]}], "]"}],
"\[Equal]", "#"}], "]"}], "&"}]}], "//", "TableForm"}]}], "Input",
CellChangeTimes->{{3.75115407229806*^9, 3.7511540907068577`*^9}, {
3.7511541283727303`*^9, 3.751154145140655*^9}, {3.751154185947949*^9,
3.751154186454541*^9}, {3.75115440553708*^9, 3.751154408241644*^9}, {
3.751157061478198*^9, 3.75115706228808*^9}, {3.751157258647379*^9,
3.751157263470278*^9}, {3.751157318329653*^9, 3.751157382599204*^9},
3.751243546881419*^9, 3.7512446518094807`*^9, {3.751758233817813*^9,
3.75175824703872*^9}},
CellLabel->"In[21]:=",
CellID->168493274,ExpressionUUID->"1e3b22c3-818c-4665-ba95-13871edd5826"],
Cell[BoxData[
TagBox[
TagBox[GridBox[{
{
RowBox[{
RowBox[{"y", "[", "1", "]"}], "\[Equal]",
FractionBox[
RowBox[{"tk", "[", "t", "]"}],
RowBox[{"DD", " ", "Rf"}]]}]}
},
GridBoxAlignment->{
"Columns" -> {{Left}}, "ColumnsIndexed" -> {}, "Rows" -> {{Baseline}},
"RowsIndexed" -> {}},
GridBoxSpacings->{"Columns" -> {
Offset[0.27999999999999997`], {
Offset[0.5599999999999999]},
Offset[0.27999999999999997`]}, "ColumnsIndexed" -> {}, "Rows" -> {
Offset[0.2], {
Offset[0.4]},
Offset[0.2]}, "RowsIndexed" -> {}}],
Column],
Function[BoxForm`e$,
TableForm[BoxForm`e$]]]], "Output",
CellChangeTimes->{
3.751154146341402*^9, {3.7511541874520807`*^9, 3.75115421035389*^9},
3.751154409289606*^9, 3.751155063159431*^9, 3.7511551578703537`*^9, {
3.751155631171876*^9, 3.7511556669733973`*^9}, 3.751156831153225*^9, {
3.7511570444961233`*^9, 3.751157063011318*^9}, 3.7511573845757217`*^9,
3.7511576929574337`*^9, 3.751237248545403*^9, 3.75123731846103*^9,
3.751237558009218*^9, 3.75124126512782*^9, 3.751241341691312*^9,
3.7512414602930183`*^9, {3.7512415890389347`*^9, 3.751241614550478*^9},
3.751243391757332*^9, {3.751243525864511*^9, 3.7512435473077707`*^9},
3.7512436090180597`*^9, 3.751243754645124*^9, 3.751243860169832*^9,
3.751244182055914*^9, 3.751244399162704*^9, 3.751244486514542*^9,
3.751244523619996*^9, 3.7512445593242817`*^9, {3.751244640162719*^9,
3.7512446526642714`*^9}, 3.7512450693463917`*^9, 3.751245129916429*^9,
3.751245199944006*^9, 3.7517581589332533`*^9, 3.751758247484844*^9,
3.7518117541350803`*^9, 3.7518120056910477`*^9, {3.7518120796408854`*^9,
3.751812087500362*^9}, {3.751812187588044*^9, 3.751812209829607*^9},
3.751824464419462*^9, 3.751827255574244*^9, 3.751859641437621*^9,
3.7518599517115088`*^9, 3.751859982075001*^9, {3.751861010824903*^9,
3.75186103667798*^9}, 3.752076506380468*^9, 3.7520765899367*^9},
CellLabel->
"Out[22]//TableForm=",ExpressionUUID->"6a1f7b24-08fa-4272-bc18-\
ecbf92b7613a"]
}, Open ]],
Cell["\<\
Linearize the state and output equations about the origin. Note that they \
equations are already linear, but this step is useful for finding the linear \
operators.\
\>", "Text",
CellChangeTimes->{{3.7027577023090067`*^9, 3.7027577064193897`*^9}, {
3.702757795179463*^9, 3.702757797848589*^9}, {3.702758517676444*^9,
3.702758517969734*^9}, {3.751157404122614*^9, 3.751157407793623*^9}, {
3.752076666333744*^9, 3.7520767184530783`*^9}},
CellID->1243308122,ExpressionUUID->"58bca072-25c9-412d-92e7-5062006193ff"],
Cell[CellGroupData[{
Cell[BoxData[{
RowBox[{
RowBox[{"abe", "=",
RowBox[{"linearizeState", "[",
RowBox[{"inVars", ",", "stateVars", ",", "stateEqList"}], "]"}]}],
";"}], "\[IndentingNewLine]",
RowBox[{
RowBox[{"cdf", "=",
RowBox[{"linearizeOutput", "[",
RowBox[{"inVars", ",", "stateVars", ",", "outEqList"}], "]"}]}],
";"}], "\[IndentingNewLine]",
RowBox[{
RowBox[{
RowBox[{"{",
RowBox[{"a", ",", "b", ",", "e"}], "}"}], " ", "=", " ", "abe"}],
";"}], "\[IndentingNewLine]",
RowBox[{
RowBox[{
RowBox[{"{",
RowBox[{"c", ",", "d", ",", "f"}], "}"}], " ", "=", " ", "cdf"}],
";"}], "\[IndentingNewLine]",
RowBox[{
RowBox[{"Table", "[", "\[IndentingNewLine]",
RowBox[{
RowBox[{"Print", "[",
RowBox[{
RowBox[{
RowBox[{
RowBox[{"{",
RowBox[{"\"\<A\>\"", ",", "\"\<B\>\"", ",", "\"\<E\>\""}], "}"}],
"[",
RowBox[{"[", "i", "]"}], "]"}], "<>", "\"\< = \>\""}], ",",
"\[IndentingNewLine]",
RowBox[{
RowBox[{"abe", "[",
RowBox[{"[", "i", "]"}], "]"}], "//", "MatrixForm"}]}], "]"}], ",",
"\[IndentingNewLine]",
RowBox[{"{",
RowBox[{"i", ",", "1", ",", "3"}], "}"}]}], "\[IndentingNewLine]", "]"}],
";"}], "\[IndentingNewLine]",
RowBox[{
RowBox[{"Table", "[", "\[IndentingNewLine]",
RowBox[{
RowBox[{"Print", "[",
RowBox[{
RowBox[{
RowBox[{
RowBox[{"{",
RowBox[{"\"\<C\>\"", ",", "\"\<D\>\"", ",", "\"\<F\>\""}], "}"}],
"[",
RowBox[{"[", "i", "]"}], "]"}], "<>", "\"\< = \>\""}], ",",
"\[IndentingNewLine]",
RowBox[{
RowBox[{"cdf", "[",
RowBox[{"[", "i", "]"}], "]"}], "//", "MatrixForm"}]}], "]"}], ",",
"\[IndentingNewLine]",
RowBox[{"{",
RowBox[{"i", ",", "1", ",", "3"}], "}"}]}], "\[IndentingNewLine]", "]"}],
";"}]}], "Input",
CellChangeTimes->{{3.7027439116175957`*^9, 3.702744033267743*^9}, {
3.702744268313121*^9, 3.702744305447536*^9}, {3.702756777004896*^9,
3.70275678344866*^9}, {3.702756831990849*^9, 3.702756860110464*^9}, {
3.702757530029274*^9, 3.702757693348208*^9}, {3.70275831097122*^9,
3.702758346655436*^9}, {3.702758404323063*^9, 3.7027584657658863`*^9}, {
3.702759733356806*^9, 3.702760026609954*^9}, {3.70276007934348*^9,
3.702760085749899*^9}, {3.702760162502097*^9, 3.7027602713332367`*^9}, {
3.7031743380777597`*^9, 3.703174365043808*^9}, {3.7032126025366993`*^9,
3.703212608677622*^9}, {3.7275822132423773`*^9, 3.727582355849462*^9},
3.751154035567876*^9, {3.751157216290072*^9, 3.751157229742868*^9}, {
3.751157416529111*^9, 3.751157421062249*^9}, {3.7511576066364393`*^9,
3.751157656881137*^9}, {3.7511577043529177`*^9, 3.7511577145152082`*^9}, {
3.7512372971265383`*^9, 3.7512373336944847`*^9}, {3.751237381927462*^9,
3.751237410777441*^9}, {3.75207677021564*^9, 3.75207683862414*^9}},
CellLabel->"In[24]:=",
CellID->1263721144,ExpressionUUID->"776bde41-1ddb-4711-9f56-dcb55f78bb22"],
Cell[CellGroupData[{
Cell[BoxData[
InterpretationBox[
RowBox[{"\<\"A = \"\>", "\[InvisibleSpace]",
TagBox[
RowBox[{"(", "\[NoBreak]", GridBox[{
{
FractionBox[
RowBox[{"kt", "-",
RowBox[{
SuperscriptBox["DD", "2"], " ", "kt", " ",
SuperscriptBox["Kv", "2"], " ", "R", " ", "Rf"}]}],
RowBox[{
SuperscriptBox["DD", "2"], " ",
RowBox[{"(",
RowBox[{"1", "+",
RowBox[{"kt", " ",
SuperscriptBox["Kv", "2"], " ", "L"}]}], ")"}], " ", "Rf"}]]}
},
GridBoxAlignment->{
"Columns" -> {{Center}}, "ColumnsIndexed" -> {},
"Rows" -> {{Baseline}}, "RowsIndexed" -> {}},
GridBoxSpacings->{"Columns" -> {
Offset[0.27999999999999997`], {
Offset[0.7]},
Offset[0.27999999999999997`]}, "ColumnsIndexed" -> {}, "Rows" -> {
Offset[0.2], {
Offset[0.4]},
Offset[0.2]}, "RowsIndexed" -> {}}], "\[NoBreak]", ")"}],
Function[BoxForm`e$,
MatrixForm[BoxForm`e$]]]}],
SequenceForm["A = ",
MatrixForm[{{$CellContext`DD^(-2) (
1 + $CellContext`kt $CellContext`Kv^2 $CellContext`L)^(-1) \
$CellContext`Rf^(-1) ($CellContext`kt - $CellContext`DD^2 $CellContext`kt \
$CellContext`Kv^2 $CellContext`R $CellContext`Rf)}}]],
Editable->False]], "Print",
CellChangeTimes->{
3.75124360906213*^9, 3.7512437547528667`*^9, 3.751243860273984*^9,
3.751244182134783*^9, 3.751244399232733*^9, 3.751244486581511*^9,
3.75124455936265*^9, 3.7512446402133417`*^9, 3.751244670741263*^9,
3.751245069381319*^9, 3.751245129946912*^9, 3.7512451999957867`*^9,
3.751758253230949*^9, 3.751811754203699*^9, 3.751812005763603*^9, {
3.751812079703279*^9, 3.751812087609988*^9}, {3.751812187660625*^9,
3.751812209901985*^9}, 3.751824464483783*^9, 3.751827255619883*^9,
3.7518596414861507`*^9, 3.751859982209362*^9, {3.751861010877256*^9,
3.7518610367788982`*^9}, 3.752076842423059*^9},
CellLabel->
"During evaluation of \
In[24]:=",ExpressionUUID->"5a584e1e-1820-4ebb-9b5b-e42761396b5c"],
Cell[BoxData[
InterpretationBox[
RowBox[{"\<\"B = \"\>", "\[InvisibleSpace]",
TagBox[
RowBox[{"(", "\[NoBreak]", GridBox[{
{
FractionBox[
RowBox[{"kt", " ", "Kv"}],
RowBox[{"1", "+",
RowBox[{"kt", " ",
SuperscriptBox["Kv", "2"], " ", "L"}]}]]}
},
GridBoxAlignment->{
"Columns" -> {{Center}}, "ColumnsIndexed" -> {},
"Rows" -> {{Baseline}}, "RowsIndexed" -> {}},
GridBoxSpacings->{"Columns" -> {
Offset[0.27999999999999997`], {
Offset[0.7]},
Offset[0.27999999999999997`]}, "ColumnsIndexed" -> {}, "Rows" -> {
Offset[0.2], {
Offset[0.4]},
Offset[0.2]}, "RowsIndexed" -> {}}], "\[NoBreak]", ")"}],
Function[BoxForm`e$,
MatrixForm[BoxForm`e$]]]}],
SequenceForm["B = ",
MatrixForm[{{$CellContext`kt $CellContext`Kv/(
1 + $CellContext`kt $CellContext`Kv^2 $CellContext`L)}}]],
Editable->False]], "Print",
CellChangeTimes->{
3.75124360906213*^9, 3.7512437547528667`*^9, 3.751243860273984*^9,
3.751244182134783*^9, 3.751244399232733*^9, 3.751244486581511*^9,
3.75124455936265*^9, 3.7512446402133417`*^9, 3.751244670741263*^9,
3.751245069381319*^9, 3.751245129946912*^9, 3.7512451999957867`*^9,
3.751758253230949*^9, 3.751811754203699*^9, 3.751812005763603*^9, {
3.751812079703279*^9, 3.751812087609988*^9}, {3.751812187660625*^9,
3.751812209901985*^9}, 3.751824464483783*^9, 3.751827255619883*^9,
3.7518596414861507`*^9, 3.751859982209362*^9, {3.751861010877256*^9,
3.7518610367788982`*^9}, 3.752076842435828*^9},
CellLabel->
"During evaluation of \
In[24]:=",ExpressionUUID->"cf9bae7a-6ac1-4757-a1f1-03f12fc238a0"],
Cell[BoxData[
InterpretationBox[
RowBox[{"\<\"E = \"\>", "\[InvisibleSpace]",
TagBox[
RowBox[{"(", "\[NoBreak]", GridBox[{
{"0"}
},
GridBoxAlignment->{
"Columns" -> {{Center}}, "ColumnsIndexed" -> {},
"Rows" -> {{Baseline}}, "RowsIndexed" -> {}},
GridBoxSpacings->{"Columns" -> {
Offset[0.27999999999999997`], {
Offset[0.7]},
Offset[0.27999999999999997`]}, "ColumnsIndexed" -> {}, "Rows" -> {
Offset[0.2], {
Offset[0.4]},
Offset[0.2]}, "RowsIndexed" -> {}}], "\[NoBreak]", ")"}],
Function[BoxForm`e$,
MatrixForm[BoxForm`e$]]]}],
SequenceForm["E = ",
MatrixForm[{{0}}]],
Editable->False]], "Print",
CellChangeTimes->{
3.75124360906213*^9, 3.7512437547528667`*^9, 3.751243860273984*^9,
3.751244182134783*^9, 3.751244399232733*^9, 3.751244486581511*^9,
3.75124455936265*^9, 3.7512446402133417`*^9, 3.751244670741263*^9,
3.751245069381319*^9, 3.751245129946912*^9, 3.7512451999957867`*^9,
3.751758253230949*^9, 3.751811754203699*^9, 3.751812005763603*^9, {
3.751812079703279*^9, 3.751812087609988*^9}, {3.751812187660625*^9,
3.751812209901985*^9}, 3.751824464483783*^9, 3.751827255619883*^9,
3.7518596414861507`*^9, 3.751859982209362*^9, {3.751861010877256*^9,
3.7518610367788982`*^9}, 3.752076842445363*^9},
CellLabel->
"During evaluation of \
In[24]:=",ExpressionUUID->"21af23fc-80c2-4ab0-862c-09a11983b55d"],
Cell[BoxData[
InterpretationBox[
RowBox[{"\<\"C = \"\>", "\[InvisibleSpace]",
TagBox[
RowBox[{"(", "\[NoBreak]", GridBox[{
{
FractionBox["1",
RowBox[{"DD", " ", "Rf"}]]}
},
GridBoxAlignment->{
"Columns" -> {{Center}}, "ColumnsIndexed" -> {},
"Rows" -> {{Baseline}}, "RowsIndexed" -> {}},
GridBoxSpacings->{"Columns" -> {
Offset[0.27999999999999997`], {
Offset[0.7]},
Offset[0.27999999999999997`]}, "ColumnsIndexed" -> {}, "Rows" -> {
Offset[0.2], {
Offset[0.4]},
Offset[0.2]}, "RowsIndexed" -> {}}], "\[NoBreak]", ")"}],
Function[BoxForm`e$,
MatrixForm[BoxForm`e$]]]}],
SequenceForm["C = ",
MatrixForm[{{$CellContext`DD^(-1)/$CellContext`Rf}}]],
Editable->False]], "Print",
CellChangeTimes->{
3.75124360906213*^9, 3.7512437547528667`*^9, 3.751243860273984*^9,
3.751244182134783*^9, 3.751244399232733*^9, 3.751244486581511*^9,
3.75124455936265*^9, 3.7512446402133417`*^9, 3.751244670741263*^9,
3.751245069381319*^9, 3.751245129946912*^9, 3.7512451999957867`*^9,
3.751758253230949*^9, 3.751811754203699*^9, 3.751812005763603*^9, {
3.751812079703279*^9, 3.751812087609988*^9}, {3.751812187660625*^9,
3.751812209901985*^9}, 3.751824464483783*^9, 3.751827255619883*^9,
3.7518596414861507`*^9, 3.751859982209362*^9, {3.751861010877256*^9,
3.7518610367788982`*^9}, 3.752076842450294*^9},
CellLabel->
"During evaluation of \
In[24]:=",ExpressionUUID->"8cbf0e3d-ede6-49e2-af98-02eac40d7f4a"],
Cell[BoxData[
InterpretationBox[
RowBox[{"\<\"D = \"\>", "\[InvisibleSpace]",
TagBox[
RowBox[{"(", "\[NoBreak]", GridBox[{
{"0"}
},
GridBoxAlignment->{
"Columns" -> {{Center}}, "ColumnsIndexed" -> {},
"Rows" -> {{Baseline}}, "RowsIndexed" -> {}},
GridBoxSpacings->{"Columns" -> {
Offset[0.27999999999999997`], {
Offset[0.7]},
Offset[0.27999999999999997`]}, "ColumnsIndexed" -> {}, "Rows" -> {
Offset[0.2], {
Offset[0.4]},
Offset[0.2]}, "RowsIndexed" -> {}}], "\[NoBreak]", ")"}],
Function[BoxForm`e$,
MatrixForm[BoxForm`e$]]]}],
SequenceForm["D = ",
MatrixForm[{{0}}]],
Editable->False]], "Print",
CellChangeTimes->{
3.75124360906213*^9, 3.7512437547528667`*^9, 3.751243860273984*^9,
3.751244182134783*^9, 3.751244399232733*^9, 3.751244486581511*^9,
3.75124455936265*^9, 3.7512446402133417`*^9, 3.751244670741263*^9,
3.751245069381319*^9, 3.751245129946912*^9, 3.7512451999957867`*^9,
3.751758253230949*^9, 3.751811754203699*^9, 3.751812005763603*^9, {
3.751812079703279*^9, 3.751812087609988*^9}, {3.751812187660625*^9,
3.751812209901985*^9}, 3.751824464483783*^9, 3.751827255619883*^9,
3.7518596414861507`*^9, 3.751859982209362*^9, {3.751861010877256*^9,
3.7518610367788982`*^9}, 3.7520768424550543`*^9},
CellLabel->
"During evaluation of \
In[24]:=",ExpressionUUID->"fe6560ac-d336-46b8-9697-b73087b050c9"],
Cell[BoxData[
InterpretationBox[
RowBox[{"\<\"F = \"\>", "\[InvisibleSpace]",
TagBox[
RowBox[{"(", "\[NoBreak]", GridBox[{
{"0"}
},
GridBoxAlignment->{
"Columns" -> {{Center}}, "ColumnsIndexed" -> {},
"Rows" -> {{Baseline}}, "RowsIndexed" -> {}},
GridBoxSpacings->{"Columns" -> {
Offset[0.27999999999999997`], {
Offset[0.7]},
Offset[0.27999999999999997`]}, "ColumnsIndexed" -> {}, "Rows" -> {
Offset[0.2], {
Offset[0.4]},
Offset[0.2]}, "RowsIndexed" -> {}}], "\[NoBreak]", ")"}],
Function[BoxForm`e$,
MatrixForm[BoxForm`e$]]]}],
SequenceForm["F = ",
MatrixForm[{{0}}]],
Editable->False]], "Print",
CellChangeTimes->{
3.75124360906213*^9, 3.7512437547528667`*^9, 3.751243860273984*^9,
3.751244182134783*^9, 3.751244399232733*^9, 3.751244486581511*^9,
3.75124455936265*^9, 3.7512446402133417`*^9, 3.751244670741263*^9,
3.751245069381319*^9, 3.751245129946912*^9, 3.7512451999957867`*^9,
3.751758253230949*^9, 3.751811754203699*^9, 3.751812005763603*^9, {
3.751812079703279*^9, 3.751812087609988*^9}, {3.751812187660625*^9,
3.751812209901985*^9}, 3.751824464483783*^9, 3.751827255619883*^9,
3.7518596414861507`*^9, 3.751859982209362*^9, {3.751861010877256*^9,
3.7518610367788982`*^9}, 3.752076842459983*^9},
CellLabel->
"During evaluation of \
In[24]:=",ExpressionUUID->"495b5093-ca5c-49d9-96fc-a915ecf2516d"]
}, Open ]]
}, Open ]],
Cell["\<\
The results can also be viewed in other languages such as LaTeX,\
\>", "Text",
CellChangeTimes->{{3.749519004775931*^9,
3.7495190192240267`*^9}},ExpressionUUID->"a6fbf9b4-2792-4e36-ba8a-\
ccb4c87d8b63"],
Cell[CellGroupData[{
Cell[BoxData[
RowBox[{"a", "//", "TeXForm"}]], "Input",
CellChangeTimes->{{3.749519022571344*^9, 3.749519061905691*^9}},
CellLabel->"In[30]:=",ExpressionUUID->"17be4c9f-eade-4683-8f67-61704fc7aa5a"],
Cell["\<\
\\left(
\\begin{array}{c}
\\frac{\\text{kt}-\\text{DD}^2 \\text{kt} \\text{Kv}^2 R
\\text{Rf}}{\\text{DD}^2 \\left(\\text{kt} L \\text{Kv}^2+1\\right)
\\text{Rf}} \\\\
\\end{array}
\\right)\
\>", "Output",
CellChangeTimes->{3.7495190296879463`*^9, 3.749519062396727*^9,
3.749568910749227*^9, 3.749569418643159*^9, 3.749569500845234*^9,
3.749582797308194*^9, 3.74958311222946*^9, 3.7520768651458063`*^9},
CellLabel->
"Out[30]//TeXForm=",ExpressionUUID->"cb47caff-2dec-4068-acc0-2c02afc2be95"]
}, Open ]],
Cell["\<\
Other Mathematica functions can be used to manipulate the equations. At this \
point let\[CloseCurlyQuote]s solve the differential equation,\
\>", "Text",
CellChangeTimes->{{3.749519149185121*^9, 3.7495191771852207`*^9}, {
3.749519858089747*^9,
3.7495198881295958`*^9}},ExpressionUUID->"7c2257b8-5345-427b-8414-\
02f236e5c37e"],
Cell[CellGroupData[{
Cell[BoxData[{
RowBox[{
RowBox[{"input", "=",
RowBox[{