Merge branch 'fixf'
[wrf-fire-matlab.git] / quicwind / symbolic_log_feb11.txt
bloba7914e1e63f7296f0315f676f7878583bad4b664
1 X = regular_mesh([5,5,5],[1,1,1]/5,1);
3 X =
4   1×3 <a href="matlab:helpPopup cell" style="font-weight:bold">cell</a> array
5     {6×6×6 double}    {6×6×6 double}    {6×6×6 double}
6 edit mass_cons_flux_test
7 edit mass_cons_flux
8 edit mass_cons_flux
9 edit mat_gen_wind_flux_div
10 DM=mat_gen_wind_flux_div(X)
11 type =
12     'DM'
13 DM =
14    (1,1)      -0.0400
15    (1,2)       0.0400
16    (2,2)      -0.0400
17    (2,3)       0.0400
18    (3,3)      -0.0400
19    (3,4)       0.0400
20    (4,4)      -0.0400
21    (4,5)       0.0400
22    (5,5)      -0.0400
23    (5,6)       0.0400
24    (6,7)      -0.0400
25    (6,8)       0.0400
26    (7,8)      -0.0400
27    (7,9)       0.0400
28    (8,9)      -0.0400
29    (8,10)      0.0400
30    (9,10)     -0.0400
31    (9,11)      0.0400
32   (10,11)     -0.0400
33   (10,12)      0.0400
34   (11,13)     -0.0400
35   (11,14)      0.0400
36   (12,14)     -0.0400
37   (12,15)      0.0400
38   (13,15)     -0.0400
39   (13,16)      0.0400
40   (14,16)     -0.0400
41   (14,17)      0.0400
42   (15,17)     -0.0400
43   (15,18)      0.0400
44   (16,19)     -0.0400
45   (16,20)      0.0400
46   (17,20)     -0.0400
47   (17,21)      0.0400
48   (18,21)     -0.0400
49   (18,22)      0.0400
50   (19,22)     -0.0400
51   (19,23)      0.0400
52   (20,23)     -0.0400
53   (20,24)      0.0400
54   (21,25)     -0.0400
55   (21,26)      0.0400
56   (22,26)     -0.0400
57   (22,27)      0.0400
58   (23,27)     -0.0400
59   (23,28)      0.0400
60   (24,28)     -0.0400
61   (24,29)      0.0400
62   (25,29)     -0.0400
63   (25,30)      0.0400
64   (26,31)     -0.0400
65   (26,32)      0.0400
66   (27,32)     -0.0400
67   (27,33)      0.0400
68   (28,33)     -0.0400
69   (28,34)      0.0400
70   (29,34)     -0.0400
71   (29,35)      0.0400
72   (30,35)     -0.0400
73   (30,36)      0.0400
74   (31,37)     -0.0400
75   (31,38)      0.0400
76   (32,38)     -0.0400
77   (32,39)      0.0400
78   (33,39)     -0.0400
79   (33,40)      0.0400
80   (34,40)     -0.0400
81   (34,41)      0.0400
82   (35,41)     -0.0400
83   (35,42)      0.0400
84   (36,43)     -0.0400
85   (36,44)      0.0400
86   (37,44)     -0.0400
87   (37,45)      0.0400
88   (38,45)     -0.0400
89   (38,46)      0.0400
90   (39,46)     -0.0400
91   (39,47)      0.0400
92   (40,47)     -0.0400
93   (40,48)      0.0400
94   (41,49)     -0.0400
95   (41,50)      0.0400
96   (42,50)     -0.0400
97   (42,51)      0.0400
98   (43,51)     -0.0400
99   (43,52)      0.0400
100   (44,52)     -0.0400
101   (44,53)      0.0400
102   (45,53)     -0.0400
103   (45,54)      0.0400
104   (46,55)     -0.0400
105   (46,56)      0.0400
106   (47,56)     -0.0400
107   (47,57)      0.0400
108   (48,57)     -0.0400
109   (48,58)      0.0400
110   (49,58)     -0.0400
111   (49,59)      0.0400
112   (50,59)     -0.0400
113   (50,60)      0.0400
114   (51,61)     -0.0400
115   (51,62)      0.0400
116   (52,62)     -0.0400
117   (52,63)      0.0400
118   (53,63)     -0.0400
119   (53,64)      0.0400
120   (54,64)     -0.0400
121   (54,65)      0.0400
122   (55,65)     -0.0400
123   (55,66)      0.0400
124   (56,67)     -0.0400
125   (56,68)      0.0400
126   (57,68)     -0.0400
127   (57,69)      0.0400
128   (58,69)     -0.0400
129   (58,70)      0.0400
130   (59,70)     -0.0400
131   (59,71)      0.0400
132   (60,71)     -0.0400
133   (60,72)      0.0400
134   (61,73)     -0.0400
135   (61,74)      0.0400
136   (62,74)     -0.0400
137   (62,75)      0.0400
138   (63,75)     -0.0400
139   (63,76)      0.0400
140   (64,76)     -0.0400
141   (64,77)      0.0400
142   (65,77)     -0.0400
143   (65,78)      0.0400
144   (66,79)     -0.0400
145   (66,80)      0.0400
146   (67,80)     -0.0400
147   (67,81)      0.0400
148   (68,81)     -0.0400
149   (68,82)      0.0400
150   (69,82)     -0.0400
151   (69,83)      0.0400
152   (70,83)     -0.0400
153   (70,84)      0.0400
154   (71,85)     -0.0400
155   (71,86)      0.0400
156   (72,86)     -0.0400
157   (72,87)      0.0400
158   (73,87)     -0.0400
159   (73,88)      0.0400
160   (74,88)     -0.0400
161   (74,89)      0.0400
162   (75,89)     -0.0400
163   (75,90)      0.0400
164   (76,91)     -0.0400
165   (76,92)      0.0400
166   (77,92)     -0.0400
167   (77,93)      0.0400
168   (78,93)     -0.0400
169   (78,94)      0.0400
170   (79,94)     -0.0400
171   (79,95)      0.0400
172   (80,95)     -0.0400
173   (80,96)      0.0400
174   (81,97)     -0.0400
175   (81,98)      0.0400
176   (82,98)     -0.0400
177   (82,99)      0.0400
178   (83,99)     -0.0400
179   (83,100)     0.0400
180   (84,100)    -0.0400
181   (84,101)     0.0400
182   (85,101)    -0.0400
183   (85,102)     0.0400
184   (86,103)    -0.0400
185   (86,104)     0.0400
186   (87,104)    -0.0400
187   (87,105)     0.0400
188   (88,105)    -0.0400
189   (88,106)     0.0400
190   (89,106)    -0.0400
191   (89,107)     0.0400
192   (90,107)    -0.0400
193   (90,108)     0.0400
194   (91,109)    -0.0400
195   (91,110)     0.0400
196   (92,110)    -0.0400
197   (92,111)     0.0400
198   (93,111)    -0.0400
199   (93,112)     0.0400
200   (94,112)    -0.0400
201   (94,113)     0.0400
202   (95,113)    -0.0400
203   (95,114)     0.0400
204   (96,115)    -0.0400
205   (96,116)     0.0400
206   (97,116)    -0.0400
207   (97,117)     0.0400
208   (98,117)    -0.0400
209   (98,118)     0.0400
210   (99,118)    -0.0400
211   (99,119)     0.0400
212  (100,119)    -0.0400
213  (100,120)     0.0400
214  (101,121)    -0.0400
215  (101,122)     0.0400
216  (102,122)    -0.0400
217  (102,123)     0.0400
218  (103,123)    -0.0400
219  (103,124)     0.0400
220  (104,124)    -0.0400
221  (104,125)     0.0400
222  (105,125)    -0.0400
223  (105,126)     0.0400
224  (106,127)    -0.0400
225  (106,128)     0.0400
226  (107,128)    -0.0400
227  (107,129)     0.0400
228  (108,129)    -0.0400
229  (108,130)     0.0400
230  (109,130)    -0.0400
231  (109,131)     0.0400
232  (110,131)    -0.0400
233  (110,132)     0.0400
234  (111,133)    -0.0400
235  (111,134)     0.0400
236  (112,134)    -0.0400
237  (112,135)     0.0400
238  (113,135)    -0.0400
239  (113,136)     0.0400
240  (114,136)    -0.0400
241  (114,137)     0.0400
242  (115,137)    -0.0400
243  (115,138)     0.0400
244  (116,139)    -0.0400
245  (116,140)     0.0400
246  (117,140)    -0.0400
247  (117,141)     0.0400
248  (118,141)    -0.0400
249  (118,142)     0.0400
250  (119,142)    -0.0400
251  (119,143)     0.0400
252  (120,143)    -0.0400
253  (120,144)     0.0400
254  (121,145)    -0.0400
255  (121,146)     0.0400
256  (122,146)    -0.0400
257  (122,147)     0.0400
258  (123,147)    -0.0400
259  (123,148)     0.0400
260  (124,148)    -0.0400
261  (124,149)     0.0400
262  (125,149)    -0.0400
263  (125,150)     0.0400
264    (1,151)    -0.0400
265    (2,152)    -0.0400
266    (3,153)    -0.0400
267    (4,154)    -0.0400
268    (5,155)    -0.0400
269    (1,156)     0.0400
270    (6,156)    -0.0400
271    (2,157)     0.0400
272    (7,157)    -0.0400
273    (3,158)     0.0400
274    (8,158)    -0.0400
275    (4,159)     0.0400
276    (9,159)    -0.0400
277    (5,160)     0.0400
278   (10,160)    -0.0400
279    (6,161)     0.0400
280   (11,161)    -0.0400
281    (7,162)     0.0400
282   (12,162)    -0.0400
283    (8,163)     0.0400
284   (13,163)    -0.0400
285    (9,164)     0.0400
286   (14,164)    -0.0400
287   (10,165)     0.0400
288   (15,165)    -0.0400
289   (11,166)     0.0400
290   (16,166)    -0.0400
291   (12,167)     0.0400
292   (17,167)    -0.0400
293   (13,168)     0.0400
294   (18,168)    -0.0400
295   (14,169)     0.0400
296   (19,169)    -0.0400
297   (15,170)     0.0400
298   (20,170)    -0.0400
299   (16,171)     0.0400
300   (21,171)    -0.0400
301   (17,172)     0.0400
302   (22,172)    -0.0400
303   (18,173)     0.0400
304   (23,173)    -0.0400
305   (19,174)     0.0400
306   (24,174)    -0.0400
307   (20,175)     0.0400
308   (25,175)    -0.0400
309   (21,176)     0.0400
310   (22,177)     0.0400
311   (23,178)     0.0400
312   (24,179)     0.0400
313   (25,180)     0.0400
314   (26,181)    -0.0400
315   (27,182)    -0.0400
316   (28,183)    -0.0400
317   (29,184)    -0.0400
318   (30,185)    -0.0400
319   (26,186)     0.0400
320   (31,186)    -0.0400
321   (27,187)     0.0400
322   (32,187)    -0.0400
323   (28,188)     0.0400
324   (33,188)    -0.0400
325   (29,189)     0.0400
326   (34,189)    -0.0400
327   (30,190)     0.0400
328   (35,190)    -0.0400
329   (31,191)     0.0400
330   (36,191)    -0.0400
331   (32,192)     0.0400
332   (37,192)    -0.0400
333   (33,193)     0.0400
334   (38,193)    -0.0400
335   (34,194)     0.0400
336   (39,194)    -0.0400
337   (35,195)     0.0400
338   (40,195)    -0.0400
339   (36,196)     0.0400
340   (41,196)    -0.0400
341   (37,197)     0.0400
342   (42,197)    -0.0400
343   (38,198)     0.0400
344   (43,198)    -0.0400
345   (39,199)     0.0400
346   (44,199)    -0.0400
347   (40,200)     0.0400
348   (45,200)    -0.0400
349   (41,201)     0.0400
350   (46,201)    -0.0400
351   (42,202)     0.0400
352   (47,202)    -0.0400
353   (43,203)     0.0400
354   (48,203)    -0.0400
355   (44,204)     0.0400
356   (49,204)    -0.0400
357   (45,205)     0.0400
358   (50,205)    -0.0400
359   (46,206)     0.0400
360   (47,207)     0.0400
361   (48,208)     0.0400
362   (49,209)     0.0400
363   (50,210)     0.0400
364   (51,211)    -0.0400
365   (52,212)    -0.0400
366   (53,213)    -0.0400
367   (54,214)    -0.0400
368   (55,215)    -0.0400
369   (51,216)     0.0400
370   (56,216)    -0.0400
371   (52,217)     0.0400
372   (57,217)    -0.0400
373   (53,218)     0.0400
374   (58,218)    -0.0400
375   (54,219)     0.0400
376   (59,219)    -0.0400
377   (55,220)     0.0400
378   (60,220)    -0.0400
379   (56,221)     0.0400
380   (61,221)    -0.0400
381   (57,222)     0.0400
382   (62,222)    -0.0400
383   (58,223)     0.0400
384   (63,223)    -0.0400
385   (59,224)     0.0400
386   (64,224)    -0.0400
387   (60,225)     0.0400
388   (65,225)    -0.0400
389   (61,226)     0.0400
390   (66,226)    -0.0400
391   (62,227)     0.0400
392   (67,227)    -0.0400
393   (63,228)     0.0400
394   (68,228)    -0.0400
395   (64,229)     0.0400
396   (69,229)    -0.0400
397   (65,230)     0.0400
398   (70,230)    -0.0400
399   (66,231)     0.0400
400   (71,231)    -0.0400
401   (67,232)     0.0400
402   (72,232)    -0.0400
403   (68,233)     0.0400
404   (73,233)    -0.0400
405   (69,234)     0.0400
406   (74,234)    -0.0400
407   (70,235)     0.0400
408   (75,235)    -0.0400
409   (71,236)     0.0400
410   (72,237)     0.0400
411   (73,238)     0.0400
412   (74,239)     0.0400
413   (75,240)     0.0400
414   (76,241)    -0.0400
415   (77,242)    -0.0400
416   (78,243)    -0.0400
417   (79,244)    -0.0400
418   (80,245)    -0.0400
419   (76,246)     0.0400
420   (81,246)    -0.0400
421   (77,247)     0.0400
422   (82,247)    -0.0400
423   (78,248)     0.0400
424   (83,248)    -0.0400
425   (79,249)     0.0400
426   (84,249)    -0.0400
427   (80,250)     0.0400
428   (85,250)    -0.0400
429   (81,251)     0.0400
430   (86,251)    -0.0400
431   (82,252)     0.0400
432   (87,252)    -0.0400
433   (83,253)     0.0400
434   (88,253)    -0.0400
435   (84,254)     0.0400
436   (89,254)    -0.0400
437   (85,255)     0.0400
438   (90,255)    -0.0400
439   (86,256)     0.0400
440   (91,256)    -0.0400
441   (87,257)     0.0400
442   (92,257)    -0.0400
443   (88,258)     0.0400
444   (93,258)    -0.0400
445   (89,259)     0.0400
446   (94,259)    -0.0400
447   (90,260)     0.0400
448   (95,260)    -0.0400
449   (91,261)     0.0400
450   (96,261)    -0.0400
451   (92,262)     0.0400
452   (97,262)    -0.0400
453   (93,263)     0.0400
454   (98,263)    -0.0400
455   (94,264)     0.0400
456   (99,264)    -0.0400
457   (95,265)     0.0400
458  (100,265)    -0.0400
459   (96,266)     0.0400
460   (97,267)     0.0400
461   (98,268)     0.0400
462   (99,269)     0.0400
463  (100,270)     0.0400
464  (101,271)    -0.0400
465  (102,272)    -0.0400
466  (103,273)    -0.0400
467  (104,274)    -0.0400
468  (105,275)    -0.0400
469  (101,276)     0.0400
470  (106,276)    -0.0400
471  (102,277)     0.0400
472  (107,277)    -0.0400
473  (103,278)     0.0400
474  (108,278)    -0.0400
475  (104,279)     0.0400
476  (109,279)    -0.0400
477  (105,280)     0.0400
478  (110,280)    -0.0400
479  (106,281)     0.0400
480  (111,281)    -0.0400
481  (107,282)     0.0400
482  (112,282)    -0.0400
483  (108,283)     0.0400
484  (113,283)    -0.0400
485  (109,284)     0.0400
486  (114,284)    -0.0400
487  (110,285)     0.0400
488  (115,285)    -0.0400
489  (111,286)     0.0400
490  (116,286)    -0.0400
491  (112,287)     0.0400
492  (117,287)    -0.0400
493  (113,288)     0.0400
494  (118,288)    -0.0400
495  (114,289)     0.0400
496  (119,289)    -0.0400
497  (115,290)     0.0400
498  (120,290)    -0.0400
499  (116,291)     0.0400
500  (121,291)    -0.0400
501  (117,292)     0.0400
502  (122,292)    -0.0400
503  (118,293)     0.0400
504  (123,293)    -0.0400
505  (119,294)     0.0400
506  (124,294)    -0.0400
507  (120,295)     0.0400
508  (125,295)    -0.0400
509  (121,296)     0.0400
510  (122,297)     0.0400
511  (123,298)     0.0400
512  (124,299)     0.0400
513  (125,300)     0.0400
514    (1,326)     0.0400
515   (26,326)    -0.0400
516    (2,327)     0.0400
517   (27,327)    -0.0400
518    (3,328)     0.0400
519   (28,328)    -0.0400
520    (4,329)     0.0400
521   (29,329)    -0.0400
522    (5,330)     0.0400
523   (30,330)    -0.0400
524    (6,331)     0.0400
525   (31,331)    -0.0400
526    (7,332)     0.0400
527   (32,332)    -0.0400
528    (8,333)     0.0400
529   (33,333)    -0.0400
530    (9,334)     0.0400
531   (34,334)    -0.0400
532   (10,335)     0.0400
533   (35,335)    -0.0400
534   (11,336)     0.0400
535   (36,336)    -0.0400
536   (12,337)     0.0400
537   (37,337)    -0.0400
538   (13,338)     0.0400
539   (38,338)    -0.0400
540   (14,339)     0.0400
541   (39,339)    -0.0400
542   (15,340)     0.0400
543   (40,340)    -0.0400
544   (16,341)     0.0400
545   (41,341)    -0.0400
546   (17,342)     0.0400
547   (42,342)    -0.0400
548   (18,343)     0.0400
549   (43,343)    -0.0400
550   (19,344)     0.0400
551   (44,344)    -0.0400
552   (20,345)     0.0400
553   (45,345)    -0.0400
554   (21,346)     0.0400
555   (46,346)    -0.0400
556   (22,347)     0.0400
557   (47,347)    -0.0400
558   (23,348)     0.0400
559   (48,348)    -0.0400
560   (24,349)     0.0400
561   (49,349)    -0.0400
562   (25,350)     0.0400
563   (50,350)    -0.0400
564   (26,351)     0.0400
565   (51,351)    -0.0400
566   (27,352)     0.0400
567   (52,352)    -0.0400
568   (28,353)     0.0400
569   (53,353)    -0.0400
570   (29,354)     0.0400
571   (54,354)    -0.0400
572   (30,355)     0.0400
573   (55,355)    -0.0400
574   (31,356)     0.0400
575   (56,356)    -0.0400
576   (32,357)     0.0400
577   (57,357)    -0.0400
578   (33,358)     0.0400
579   (58,358)    -0.0400
580   (34,359)     0.0400
581   (59,359)    -0.0400
582   (35,360)     0.0400
583   (60,360)    -0.0400
584   (36,361)     0.0400
585   (61,361)    -0.0400
586   (37,362)     0.0400
587   (62,362)    -0.0400
588   (38,363)     0.0400
589   (63,363)    -0.0400
590   (39,364)     0.0400
591   (64,364)    -0.0400
592   (40,365)     0.0400
593   (65,365)    -0.0400
594   (41,366)     0.0400
595   (66,366)    -0.0400
596   (42,367)     0.0400
597   (67,367)    -0.0400
598   (43,368)     0.0400
599   (68,368)    -0.0400
600   (44,369)     0.0400
601   (69,369)    -0.0400
602   (45,370)     0.0400
603   (70,370)    -0.0400
604   (46,371)     0.0400
605   (71,371)    -0.0400
606   (47,372)     0.0400
607   (72,372)    -0.0400
608   (48,373)     0.0400
609   (73,373)    -0.0400
610   (49,374)     0.0400
611   (74,374)    -0.0400
612   (50,375)     0.0400
613   (75,375)    -0.0400
614   (51,376)     0.0400
615   (76,376)    -0.0400
616   (52,377)     0.0400
617   (77,377)    -0.0400
618   (53,378)     0.0400
619   (78,378)    -0.0400
620   (54,379)     0.0400
621   (79,379)    -0.0400
622   (55,380)     0.0400
623   (80,380)    -0.0400
624   (56,381)     0.0400
625   (81,381)    -0.0400
626   (57,382)     0.0400
627   (82,382)    -0.0400
628   (58,383)     0.0400
629   (83,383)    -0.0400
630   (59,384)     0.0400
631   (84,384)    -0.0400
632   (60,385)     0.0400
633   (85,385)    -0.0400
634   (61,386)     0.0400
635   (86,386)    -0.0400
636   (62,387)     0.0400
637   (87,387)    -0.0400
638   (63,388)     0.0400
639   (88,388)    -0.0400
640   (64,389)     0.0400
641   (89,389)    -0.0400
642   (65,390)     0.0400
643   (90,390)    -0.0400
644   (66,391)     0.0400
645   (91,391)    -0.0400
646   (67,392)     0.0400
647   (92,392)    -0.0400
648   (68,393)     0.0400
649   (93,393)    -0.0400
650   (69,394)     0.0400
651   (94,394)    -0.0400
652   (70,395)     0.0400
653   (95,395)    -0.0400
654   (71,396)     0.0400
655   (96,396)    -0.0400
656   (72,397)     0.0400
657   (97,397)    -0.0400
658   (73,398)     0.0400
659   (98,398)    -0.0400
660   (74,399)     0.0400
661   (99,399)    -0.0400
662   (75,400)     0.0400
663  (100,400)    -0.0400
664   (76,401)     0.0400
665  (101,401)    -0.0400
666   (77,402)     0.0400
667  (102,402)    -0.0400
668   (78,403)     0.0400
669  (103,403)    -0.0400
670   (79,404)     0.0400
671  (104,404)    -0.0400
672   (80,405)     0.0400
673  (105,405)    -0.0400
674   (81,406)     0.0400
675  (106,406)    -0.0400
676   (82,407)     0.0400
677  (107,407)    -0.0400
678   (83,408)     0.0400
679  (108,408)    -0.0400
680   (84,409)     0.0400
681  (109,409)    -0.0400
682   (85,410)     0.0400
683  (110,410)    -0.0400
684   (86,411)     0.0400
685  (111,411)    -0.0400
686   (87,412)     0.0400
687  (112,412)    -0.0400
688   (88,413)     0.0400
689  (113,413)    -0.0400
690   (89,414)     0.0400
691  (114,414)    -0.0400
692   (90,415)     0.0400
693  (115,415)    -0.0400
694   (91,416)     0.0400
695  (116,416)    -0.0400
696   (92,417)     0.0400
697  (117,417)    -0.0400
698   (93,418)     0.0400
699  (118,418)    -0.0400
700   (94,419)     0.0400
701  (119,419)    -0.0400
702   (95,420)     0.0400
703  (120,420)    -0.0400
704   (96,421)     0.0400
705  (121,421)    -0.0400
706   (97,422)     0.0400
707  (122,422)    -0.0400
708   (98,423)     0.0400
709  (123,423)    -0.0400
710   (99,424)     0.0400
711  (124,424)    -0.0400
712  (100,425)     0.0400
713  (125,425)    -0.0400
714  (101,426)     0.0400
715  (102,427)     0.0400
716  (103,428)     0.0400
717  (104,429)     0.0400
718  (105,430)     0.0400
719  (106,431)     0.0400
720  (107,432)     0.0400
721  (108,433)     0.0400
722  (109,434)     0.0400
723  (110,435)     0.0400
724  (111,436)     0.0400
725  (112,437)     0.0400
726  (113,438)     0.0400
727  (114,439)     0.0400
728  (115,440)     0.0400
729  (116,441)     0.0400
730  (117,442)     0.0400
731  (118,443)     0.0400
732  (119,444)     0.0400
733  (120,445)     0.0400
734  (121,446)     0.0400
735  (122,447)     0.0400
736  (123,448)     0.0400
737  (124,449)     0.0400
738  (125,450)     0.0400
739 DM=full(DM);
740 whos
741   Name              Size              Bytes  Class     Attributes
743   DM              125x450            450000  double              
744   X                 1x3                5520  cell                
745   ans               6x6x6              1728  double              
746   bits              1x1                   8  double              
747   c                 1x64                128  char                
748   i                 1x1                   8  double              
749   known_bits        1x1                   8  double              
750   p                 2x2x2                 8  sym                 
752 L=sym('L_',[5 5 5])
753 L(:,:,1) =
754 [ L_1_1_1, L_1_2_1, L_1_3_1, L_1_4_1, L_1_5_1]
755 [ L_2_1_1, L_2_2_1, L_2_3_1, L_2_4_1, L_2_5_1]
756 [ L_3_1_1, L_3_2_1, L_3_3_1, L_3_4_1, L_3_5_1]
757 [ L_4_1_1, L_4_2_1, L_4_3_1, L_4_4_1, L_4_5_1]
758 [ L_5_1_1, L_5_2_1, L_5_3_1, L_5_4_1, L_5_5_1]
759 L(:,:,2) =
760 [ L_1_1_2, L_1_2_2, L_1_3_2, L_1_4_2, L_1_5_2]
761 [ L_2_1_2, L_2_2_2, L_2_3_2, L_2_4_2, L_2_5_2]
762 [ L_3_1_2, L_3_2_2, L_3_3_2, L_3_4_2, L_3_5_2]
763 [ L_4_1_2, L_4_2_2, L_4_3_2, L_4_4_2, L_4_5_2]
764 [ L_5_1_2, L_5_2_2, L_5_3_2, L_5_4_2, L_5_5_2]
765 L(:,:,3) =
766 [ L_1_1_3, L_1_2_3, L_1_3_3, L_1_4_3, L_1_5_3]
767 [ L_2_1_3, L_2_2_3, L_2_3_3, L_2_4_3, L_2_5_3]
768 [ L_3_1_3, L_3_2_3, L_3_3_3, L_3_4_3, L_3_5_3]
769 [ L_4_1_3, L_4_2_3, L_4_3_3, L_4_4_3, L_4_5_3]
770 [ L_5_1_3, L_5_2_3, L_5_3_3, L_5_4_3, L_5_5_3]
771 L(:,:,4) =
772 [ L_1_1_4, L_1_2_4, L_1_3_4, L_1_4_4, L_1_5_4]
773 [ L_2_1_4, L_2_2_4, L_2_3_4, L_2_4_4, L_2_5_4]
774 [ L_3_1_4, L_3_2_4, L_3_3_4, L_3_4_4, L_3_5_4]
775 [ L_4_1_4, L_4_2_4, L_4_3_4, L_4_4_4, L_4_5_4]
776 [ L_5_1_4, L_5_2_4, L_5_3_4, L_5_4_4, L_5_5_4]
777 L(:,:,5) =
778 [ L_1_1_5, L_1_2_5, L_1_3_5, L_1_4_5, L_1_5_5]
779 [ L_2_1_5, L_2_2_5, L_2_3_5, L_2_4_5, L_2_5_5]
780 [ L_3_1_5, L_3_2_5, L_3_3_5, L_3_4_5, L_3_5_5]
781 [ L_4_1_5, L_4_2_5, L_4_3_5, L_4_4_5, L_4_5_5]
782 [ L_5_1_5, L_5_2_5, L_5_3_5, L_5_4_5, L_5_5_5]
783 help grad3z
784   u=<strong>grad3z</strong>(f,h)
785   compute gradient in 3d assuming zero boundary conditions
786   u=<strong>grad3z</strong>(f,h,b) 
787   in 3rd coordinate use reflection at bottom instead of zero
788   input:
789      f       3d array
790      h(1:3)  stepsize
791      b       use zero boundary condition on output at the bottom
793 edit grad3z
794 u=grad3z(L,[1,1,1])
795 {\bThe following error occurred converting from sym to double:
796 Unable to convert expression into double array.
797 Error in <a href="matlab:matlab.internal.language.introspective.errorDocCallback('grad3z', '/Users/jmandel/Projects/wrf-fire-matlab/quicwind/grad3z.m', 18)" style="font-weight:bold">grad3z</a> (<a href="matlab: opentoline('/Users/jmandel/Projects/wrf-fire-matlab/quicwind/grad3z.m',18,0)">line 18</a>)
798 fz(2:n(1)+1,2:n(2)+1,2:n(3)+1)=f;}\b 
799 matlab.internal.language.introspective.errorDocCallback('grad3z', '/Users/jmandel/Projects/wrf-fire-matlab/quicwind/grad3z.m', 18)
800 opentoline('/Users/jmandel/Projects/wrf-fire-matlab/quicwind/grad3z.m',18,0)
801 u=grad3z(L,[1,1,1])
802 u =
803   1×3 <a href="matlab:helpPopup cell" style="font-weight:bold">cell</a> array
804     {5×4×4 sym}    {4×5×4 sym}    {4×4×5 sym}
805 u{1}
806 ans(:,:,1) =
807 [           L_1_1_1,           L_1_2_1,           L_1_3_1,           L_1_4_1]
808 [ L_2_1_1 - L_1_1_1, L_2_2_1 - L_1_2_1, L_2_3_1 - L_1_3_1, L_2_4_1 - L_1_4_1]
809 [ L_3_1_1 - L_2_1_1, L_3_2_1 - L_2_2_1, L_3_3_1 - L_2_3_1, L_3_4_1 - L_2_4_1]
810 [ L_4_1_1 - L_3_1_1, L_4_2_1 - L_3_2_1, L_4_3_1 - L_3_3_1, L_4_4_1 - L_3_4_1]
811 [ L_5_1_1 - L_4_1_1, L_5_2_1 - L_4_2_1, L_5_3_1 - L_4_3_1, L_5_4_1 - L_4_4_1]
812 ans(:,:,2) =
813 [           L_1_1_2,           L_1_2_2,           L_1_3_2,           L_1_4_2]
814 [ L_2_1_2 - L_1_1_2, L_2_2_2 - L_1_2_2, L_2_3_2 - L_1_3_2, L_2_4_2 - L_1_4_2]
815 [ L_3_1_2 - L_2_1_2, L_3_2_2 - L_2_2_2, L_3_3_2 - L_2_3_2, L_3_4_2 - L_2_4_2]
816 [ L_4_1_2 - L_3_1_2, L_4_2_2 - L_3_2_2, L_4_3_2 - L_3_3_2, L_4_4_2 - L_3_4_2]
817 [ L_5_1_2 - L_4_1_2, L_5_2_2 - L_4_2_2, L_5_3_2 - L_4_3_2, L_5_4_2 - L_4_4_2]
818 ans(:,:,3) =
819 [           L_1_1_3,           L_1_2_3,           L_1_3_3,           L_1_4_3]
820 [ L_2_1_3 - L_1_1_3, L_2_2_3 - L_1_2_3, L_2_3_3 - L_1_3_3, L_2_4_3 - L_1_4_3]
821 [ L_3_1_3 - L_2_1_3, L_3_2_3 - L_2_2_3, L_3_3_3 - L_2_3_3, L_3_4_3 - L_2_4_3]
822 [ L_4_1_3 - L_3_1_3, L_4_2_3 - L_3_2_3, L_4_3_3 - L_3_3_3, L_4_4_3 - L_3_4_3]
823 [ L_5_1_3 - L_4_1_3, L_5_2_3 - L_4_2_3, L_5_3_3 - L_4_3_3, L_5_4_3 - L_4_4_3]
824 ans(:,:,4) =
825 [           L_1_1_4,           L_1_2_4,           L_1_3_4,           L_1_4_4]
826 [ L_2_1_4 - L_1_1_4, L_2_2_4 - L_1_2_4, L_2_3_4 - L_1_3_4, L_2_4_4 - L_1_4_4]
827 [ L_3_1_4 - L_2_1_4, L_3_2_4 - L_2_2_4, L_3_3_4 - L_2_3_4, L_3_4_4 - L_2_4_4]
828 [ L_4_1_4 - L_3_1_4, L_4_2_4 - L_3_2_4, L_4_3_4 - L_3_3_4, L_4_4_4 - L_3_4_4]
829 [ L_5_1_4 - L_4_1_4, L_5_2_4 - L_4_2_4, L_5_3_4 - L_4_3_4, L_5_4_4 - L_4_4_4]
830 edit div3
831 d=div3(u,[1 1 1])
832 d(:,:,1) =
833 [           L_1_1_2 - 6*L_1_1_1 + L_1_2_1 + L_2_1_1,           L_1_1_1 - 6*L_1_2_1 + L_1_2_2 + L_1_3_1 + L_2_2_1,           L_1_2_1 - 6*L_1_3_1 + L_1_3_2 + L_1_4_1 + L_2_3_1,           L_1_3_1 - 6*L_1_4_1 + L_1_4_2 + L_1_5_1 + L_2_4_1]
834 [ L_1_1_1 - 6*L_2_1_1 + L_2_1_2 + L_2_2_1 + L_3_1_1, L_1_2_1 + L_2_1_1 - 6*L_2_2_1 + L_2_2_2 + L_2_3_1 + L_3_2_1, L_1_3_1 + L_2_2_1 - 6*L_2_3_1 + L_2_3_2 + L_2_4_1 + L_3_3_1, L_1_4_1 + L_2_3_1 - 6*L_2_4_1 + L_2_4_2 + L_2_5_1 + L_3_4_1]
835 [ L_2_1_1 - 6*L_3_1_1 + L_3_1_2 + L_3_2_1 + L_4_1_1, L_2_2_1 + L_3_1_1 - 6*L_3_2_1 + L_3_2_2 + L_3_3_1 + L_4_2_1, L_2_3_1 + L_3_2_1 - 6*L_3_3_1 + L_3_3_2 + L_3_4_1 + L_4_3_1, L_2_4_1 + L_3_3_1 - 6*L_3_4_1 + L_3_4_2 + L_3_5_1 + L_4_4_1]
836 [ L_3_1_1 - 6*L_4_1_1 + L_4_1_2 + L_4_2_1 + L_5_1_1, L_3_2_1 + L_4_1_1 - 6*L_4_2_1 + L_4_2_2 + L_4_3_1 + L_5_2_1, L_3_3_1 + L_4_2_1 - 6*L_4_3_1 + L_4_3_2 + L_4_4_1 + L_5_3_1, L_3_4_1 + L_4_3_1 - 6*L_4_4_1 + L_4_4_2 + L_4_5_1 + L_5_4_1]
837 d(:,:,2) =
838 [           L_1_1_1 - 6*L_1_1_2 + L_1_1_3 + L_1_2_2 + L_2_1_2,           L_1_1_2 + L_1_2_1 - 6*L_1_2_2 + L_1_2_3 + L_1_3_2 + L_2_2_2,           L_1_2_2 + L_1_3_1 - 6*L_1_3_2 + L_1_3_3 + L_1_4_2 + L_2_3_2,           L_1_3_2 + L_1_4_1 - 6*L_1_4_2 + L_1_4_3 + L_1_5_2 + L_2_4_2]
839 [ L_1_1_2 + L_2_1_1 - 6*L_2_1_2 + L_2_1_3 + L_2_2_2 + L_3_1_2, L_1_2_2 + L_2_1_2 + L_2_2_1 - 6*L_2_2_2 + L_2_2_3 + L_2_3_2 + L_3_2_2, L_1_3_2 + L_2_2_2 + L_2_3_1 - 6*L_2_3_2 + L_2_3_3 + L_2_4_2 + L_3_3_2, L_1_4_2 + L_2_3_2 + L_2_4_1 - 6*L_2_4_2 + L_2_4_3 + L_2_5_2 + L_3_4_2]
840 [ L_2_1_2 + L_3_1_1 - 6*L_3_1_2 + L_3_1_3 + L_3_2_2 + L_4_1_2, L_2_2_2 + L_3_1_2 + L_3_2_1 - 6*L_3_2_2 + L_3_2_3 + L_3_3_2 + L_4_2_2, L_2_3_2 + L_3_2_2 + L_3_3_1 - 6*L_3_3_2 + L_3_3_3 + L_3_4_2 + L_4_3_2, L_2_4_2 + L_3_3_2 + L_3_4_1 - 6*L_3_4_2 + L_3_4_3 + L_3_5_2 + L_4_4_2]
841 [ L_3_1_2 + L_4_1_1 - 6*L_4_1_2 + L_4_1_3 + L_4_2_2 + L_5_1_2, L_3_2_2 + L_4_1_2 + L_4_2_1 - 6*L_4_2_2 + L_4_2_3 + L_4_3_2 + L_5_2_2, L_3_3_2 + L_4_2_2 + L_4_3_1 - 6*L_4_3_2 + L_4_3_3 + L_4_4_2 + L_5_3_2, L_3_4_2 + L_4_3_2 + L_4_4_1 - 6*L_4_4_2 + L_4_4_3 + L_4_5_2 + L_5_4_2]
842 d(:,:,3) =
843 [           L_1_1_2 - 6*L_1_1_3 + L_1_1_4 + L_1_2_3 + L_2_1_3,           L_1_1_3 + L_1_2_2 - 6*L_1_2_3 + L_1_2_4 + L_1_3_3 + L_2_2_3,           L_1_2_3 + L_1_3_2 - 6*L_1_3_3 + L_1_3_4 + L_1_4_3 + L_2_3_3,           L_1_3_3 + L_1_4_2 - 6*L_1_4_3 + L_1_4_4 + L_1_5_3 + L_2_4_3]
844 [ L_1_1_3 + L_2_1_2 - 6*L_2_1_3 + L_2_1_4 + L_2_2_3 + L_3_1_3, L_1_2_3 + L_2_1_3 + L_2_2_2 - 6*L_2_2_3 + L_2_2_4 + L_2_3_3 + L_3_2_3, L_1_3_3 + L_2_2_3 + L_2_3_2 - 6*L_2_3_3 + L_2_3_4 + L_2_4_3 + L_3_3_3, L_1_4_3 + L_2_3_3 + L_2_4_2 - 6*L_2_4_3 + L_2_4_4 + L_2_5_3 + L_3_4_3]
845 [ L_2_1_3 + L_3_1_2 - 6*L_3_1_3 + L_3_1_4 + L_3_2_3 + L_4_1_3, L_2_2_3 + L_3_1_3 + L_3_2_2 - 6*L_3_2_3 + L_3_2_4 + L_3_3_3 + L_4_2_3, L_2_3_3 + L_3_2_3 + L_3_3_2 - 6*L_3_3_3 + L_3_3_4 + L_3_4_3 + L_4_3_3, L_2_4_3 + L_3_3_3 + L_3_4_2 - 6*L_3_4_3 + L_3_4_4 + L_3_5_3 + L_4_4_3]
846 [ L_3_1_3 + L_4_1_2 - 6*L_4_1_3 + L_4_1_4 + L_4_2_3 + L_5_1_3, L_3_2_3 + L_4_1_3 + L_4_2_2 - 6*L_4_2_3 + L_4_2_4 + L_4_3_3 + L_5_2_3, L_3_3_3 + L_4_2_3 + L_4_3_2 - 6*L_4_3_3 + L_4_3_4 + L_4_4_3 + L_5_3_3, L_3_4_3 + L_4_3_3 + L_4_4_2 - 6*L_4_4_3 + L_4_4_4 + L_4_5_3 + L_5_4_3]
847 d(:,:,4) =
848 [           L_1_1_3 - 6*L_1_1_4 + L_1_1_5 + L_1_2_4 + L_2_1_4,           L_1_1_4 + L_1_2_3 - 6*L_1_2_4 + L_1_2_5 + L_1_3_4 + L_2_2_4,           L_1_2_4 + L_1_3_3 - 6*L_1_3_4 + L_1_3_5 + L_1_4_4 + L_2_3_4,           L_1_3_4 + L_1_4_3 - 6*L_1_4_4 + L_1_4_5 + L_1_5_4 + L_2_4_4]
849 [ L_1_1_4 + L_2_1_3 - 6*L_2_1_4 + L_2_1_5 + L_2_2_4 + L_3_1_4, L_1_2_4 + L_2_1_4 + L_2_2_3 - 6*L_2_2_4 + L_2_2_5 + L_2_3_4 + L_3_2_4, L_1_3_4 + L_2_2_4 + L_2_3_3 - 6*L_2_3_4 + L_2_3_5 + L_2_4_4 + L_3_3_4, L_1_4_4 + L_2_3_4 + L_2_4_3 - 6*L_2_4_4 + L_2_4_5 + L_2_5_4 + L_3_4_4]
850 [ L_2_1_4 + L_3_1_3 - 6*L_3_1_4 + L_3_1_5 + L_3_2_4 + L_4_1_4, L_2_2_4 + L_3_1_4 + L_3_2_3 - 6*L_3_2_4 + L_3_2_5 + L_3_3_4 + L_4_2_4, L_2_3_4 + L_3_2_4 + L_3_3_3 - 6*L_3_3_4 + L_3_3_5 + L_3_4_4 + L_4_3_4, L_2_4_4 + L_3_3_4 + L_3_4_3 - 6*L_3_4_4 + L_3_4_5 + L_3_5_4 + L_4_4_4]
851 [ L_3_1_4 + L_4_1_3 - 6*L_4_1_4 + L_4_1_5 + L_4_2_4 + L_5_1_4, L_3_2_4 + L_4_1_4 + L_4_2_3 - 6*L_4_2_4 + L_4_2_5 + L_4_3_4 + L_5_2_4, L_3_3_4 + L_4_2_4 + L_4_3_3 - 6*L_4_3_4 + L_4_3_5 + L_4_4_4 + L_5_3_4, L_3_4_4 + L_4_3_4 + L_4_4_3 - 6*L_4_4_4 + L_4_4_5 + L_4_5_4 + L_5_4_4]
853 d(:,:,1) =
854 [           L_1_1_2 - 6*L_1_1_1 + L_1_2_1 + L_2_1_1,           L_1_1_1 - 6*L_1_2_1 + L_1_2_2 + L_1_3_1 + L_2_2_1,           L_1_2_1 - 6*L_1_3_1 + L_1_3_2 + L_1_4_1 + L_2_3_1,           L_1_3_1 - 6*L_1_4_1 + L_1_4_2 + L_1_5_1 + L_2_4_1]
855 [ L_1_1_1 - 6*L_2_1_1 + L_2_1_2 + L_2_2_1 + L_3_1_1, L_1_2_1 + L_2_1_1 - 6*L_2_2_1 + L_2_2_2 + L_2_3_1 + L_3_2_1, L_1_3_1 + L_2_2_1 - 6*L_2_3_1 + L_2_3_2 + L_2_4_1 + L_3_3_1, L_1_4_1 + L_2_3_1 - 6*L_2_4_1 + L_2_4_2 + L_2_5_1 + L_3_4_1]
856 [ L_2_1_1 - 6*L_3_1_1 + L_3_1_2 + L_3_2_1 + L_4_1_1, L_2_2_1 + L_3_1_1 - 6*L_3_2_1 + L_3_2_2 + L_3_3_1 + L_4_2_1, L_2_3_1 + L_3_2_1 - 6*L_3_3_1 + L_3_3_2 + L_3_4_1 + L_4_3_1, L_2_4_1 + L_3_3_1 - 6*L_3_4_1 + L_3_4_2 + L_3_5_1 + L_4_4_1]
857 [ L_3_1_1 - 6*L_4_1_1 + L_4_1_2 + L_4_2_1 + L_5_1_1, L_3_2_1 + L_4_1_1 - 6*L_4_2_1 + L_4_2_2 + L_4_3_1 + L_5_2_1, L_3_3_1 + L_4_2_1 - 6*L_4_3_1 + L_4_3_2 + L_4_4_1 + L_5_3_1, L_3_4_1 + L_4_3_1 - 6*L_4_4_1 + L_4_4_2 + L_4_5_1 + L_5_4_1]
858 d(:,:,2) =
859 [           L_1_1_1 - 6*L_1_1_2 + L_1_1_3 + L_1_2_2 + L_2_1_2,           L_1_1_2 + L_1_2_1 - 6*L_1_2_2 + L_1_2_3 + L_1_3_2 + L_2_2_2,           L_1_2_2 + L_1_3_1 - 6*L_1_3_2 + L_1_3_3 + L_1_4_2 + L_2_3_2,           L_1_3_2 + L_1_4_1 - 6*L_1_4_2 + L_1_4_3 + L_1_5_2 + L_2_4_2]
860 [ L_1_1_2 + L_2_1_1 - 6*L_2_1_2 + L_2_1_3 + L_2_2_2 + L_3_1_2, L_1_2_2 + L_2_1_2 + L_2_2_1 - 6*L_2_2_2 + L_2_2_3 + L_2_3_2 + L_3_2_2, L_1_3_2 + L_2_2_2 + L_2_3_1 - 6*L_2_3_2 + L_2_3_3 + L_2_4_2 + L_3_3_2, L_1_4_2 + L_2_3_2 + L_2_4_1 - 6*L_2_4_2 + L_2_4_3 + L_2_5_2 + L_3_4_2]
861 [ L_2_1_2 + L_3_1_1 - 6*L_3_1_2 + L_3_1_3 + L_3_2_2 + L_4_1_2, L_2_2_2 + L_3_1_2 + L_3_2_1 - 6*L_3_2_2 + L_3_2_3 + L_3_3_2 + L_4_2_2, L_2_3_2 + L_3_2_2 + L_3_3_1 - 6*L_3_3_2 + L_3_3_3 + L_3_4_2 + L_4_3_2, L_2_4_2 + L_3_3_2 + L_3_4_1 - 6*L_3_4_2 + L_3_4_3 + L_3_5_2 + L_4_4_2]
862 [ L_3_1_2 + L_4_1_1 - 6*L_4_1_2 + L_4_1_3 + L_4_2_2 + L_5_1_2, L_3_2_2 + L_4_1_2 + L_4_2_1 - 6*L_4_2_2 + L_4_2_3 + L_4_3_2 + L_5_2_2, L_3_3_2 + L_4_2_2 + L_4_3_1 - 6*L_4_3_2 + L_4_3_3 + L_4_4_2 + L_5_3_2, L_3_4_2 + L_4_3_2 + L_4_4_1 - 6*L_4_4_2 + L_4_4_3 + L_4_5_2 + L_5_4_2]
863 d(:,:,3) =
864 [           L_1_1_2 - 6*L_1_1_3 + L_1_1_4 + L_1_2_3 + L_2_1_3,           L_1_1_3 + L_1_2_2 - 6*L_1_2_3 + L_1_2_4 + L_1_3_3 + L_2_2_3,           L_1_2_3 + L_1_3_2 - 6*L_1_3_3 + L_1_3_4 + L_1_4_3 + L_2_3_3,           L_1_3_3 + L_1_4_2 - 6*L_1_4_3 + L_1_4_4 + L_1_5_3 + L_2_4_3]
865 [ L_1_1_3 + L_2_1_2 - 6*L_2_1_3 + L_2_1_4 + L_2_2_3 + L_3_1_3, L_1_2_3 + L_2_1_3 + L_2_2_2 - 6*L_2_2_3 + L_2_2_4 + L_2_3_3 + L_3_2_3, L_1_3_3 + L_2_2_3 + L_2_3_2 - 6*L_2_3_3 + L_2_3_4 + L_2_4_3 + L_3_3_3, L_1_4_3 + L_2_3_3 + L_2_4_2 - 6*L_2_4_3 + L_2_4_4 + L_2_5_3 + L_3_4_3]
866 [ L_2_1_3 + L_3_1_2 - 6*L_3_1_3 + L_3_1_4 + L_3_2_3 + L_4_1_3, L_2_2_3 + L_3_1_3 + L_3_2_2 - 6*L_3_2_3 + L_3_2_4 + L_3_3_3 + L_4_2_3, L_2_3_3 + L_3_2_3 + L_3_3_2 - 6*L_3_3_3 + L_3_3_4 + L_3_4_3 + L_4_3_3, L_2_4_3 + L_3_3_3 + L_3_4_2 - 6*L_3_4_3 + L_3_4_4 + L_3_5_3 + L_4_4_3]
867 [ L_3_1_3 + L_4_1_2 - 6*L_4_1_3 + L_4_1_4 + L_4_2_3 + L_5_1_3, L_3_2_3 + L_4_1_3 + L_4_2_2 - 6*L_4_2_3 + L_4_2_4 + L_4_3_3 + L_5_2_3, L_3_3_3 + L_4_2_3 + L_4_3_2 - 6*L_4_3_3 + L_4_3_4 + L_4_4_3 + L_5_3_3, L_3_4_3 + L_4_3_3 + L_4_4_2 - 6*L_4_4_3 + L_4_4_4 + L_4_5_3 + L_5_4_3]
868 d(:,:,4) =
869 [           L_1_1_3 - 6*L_1_1_4 + L_1_1_5 + L_1_2_4 + L_2_1_4,           L_1_1_4 + L_1_2_3 - 6*L_1_2_4 + L_1_2_5 + L_1_3_4 + L_2_2_4,           L_1_2_4 + L_1_3_3 - 6*L_1_3_4 + L_1_3_5 + L_1_4_4 + L_2_3_4,           L_1_3_4 + L_1_4_3 - 6*L_1_4_4 + L_1_4_5 + L_1_5_4 + L_2_4_4]
870 [ L_1_1_4 + L_2_1_3 - 6*L_2_1_4 + L_2_1_5 + L_2_2_4 + L_3_1_4, L_1_2_4 + L_2_1_4 + L_2_2_3 - 6*L_2_2_4 + L_2_2_5 + L_2_3_4 + L_3_2_4, L_1_3_4 + L_2_2_4 + L_2_3_3 - 6*L_2_3_4 + L_2_3_5 + L_2_4_4 + L_3_3_4, L_1_4_4 + L_2_3_4 + L_2_4_3 - 6*L_2_4_4 + L_2_4_5 + L_2_5_4 + L_3_4_4]
871 [ L_2_1_4 + L_3_1_3 - 6*L_3_1_4 + L_3_1_5 + L_3_2_4 + L_4_1_4, L_2_2_4 + L_3_1_4 + L_3_2_3 - 6*L_3_2_4 + L_3_2_5 + L_3_3_4 + L_4_2_4, L_2_3_4 + L_3_2_4 + L_3_3_3 - 6*L_3_3_4 + L_3_3_5 + L_3_4_4 + L_4_3_4, L_2_4_4 + L_3_3_4 + L_3_4_3 - 6*L_3_4_4 + L_3_4_5 + L_3_5_4 + L_4_4_4]
872 [ L_3_1_4 + L_4_1_3 - 6*L_4_1_4 + L_4_1_5 + L_4_2_4 + L_5_1_4, L_3_2_4 + L_4_1_4 + L_4_2_3 - 6*L_4_2_4 + L_4_2_5 + L_4_3_4 + L_5_2_4, L_3_3_4 + L_4_2_4 + L_4_3_3 - 6*L_4_3_4 + L_4_3_5 + L_4_4_4 + L_5_3_4, L_3_4_4 + L_4_3_4 + L_4_4_3 - 6*L_4_4_4 + L_4_4_5 + L_4_5_4 + L_5_4_4]
873 d(3,3,3)
874 ans =
875 L_2_3_3 + L_3_2_3 + L_3_3_2 - 6*L_3_3_3 + L_3_3_4 + L_3_4_3 + L_4_3_3
876 d(3,3,1)
877 ans =
878 L_2_3_1 + L_3_2_1 - 6*L_3_3_1 + L_3_3_2 + L_3_4_1 + L_4_3_1
879 u=grad3z(L,[1,1,1],0)
880 u =
881   1×3 <a href="matlab:helpPopup cell" style="font-weight:bold">cell</a> array
882     {5×4×4 sym}    {4×5×4 sym}    {4×4×5 sym}
883 d=div3(u,[1 1 1])
884 d(:,:,1) =
885 [           L_1_1_2 - 6*L_1_1_1 + L_1_2_1 + L_2_1_1,           L_1_1_1 - 6*L_1_2_1 + L_1_2_2 + L_1_3_1 + L_2_2_1,           L_1_2_1 - 6*L_1_3_1 + L_1_3_2 + L_1_4_1 + L_2_3_1,           L_1_3_1 - 6*L_1_4_1 + L_1_4_2 + L_1_5_1 + L_2_4_1]
886 [ L_1_1_1 - 6*L_2_1_1 + L_2_1_2 + L_2_2_1 + L_3_1_1, L_1_2_1 + L_2_1_1 - 6*L_2_2_1 + L_2_2_2 + L_2_3_1 + L_3_2_1, L_1_3_1 + L_2_2_1 - 6*L_2_3_1 + L_2_3_2 + L_2_4_1 + L_3_3_1, L_1_4_1 + L_2_3_1 - 6*L_2_4_1 + L_2_4_2 + L_2_5_1 + L_3_4_1]
887 [ L_2_1_1 - 6*L_3_1_1 + L_3_1_2 + L_3_2_1 + L_4_1_1, L_2_2_1 + L_3_1_1 - 6*L_3_2_1 + L_3_2_2 + L_3_3_1 + L_4_2_1, L_2_3_1 + L_3_2_1 - 6*L_3_3_1 + L_3_3_2 + L_3_4_1 + L_4_3_1, L_2_4_1 + L_3_3_1 - 6*L_3_4_1 + L_3_4_2 + L_3_5_1 + L_4_4_1]
888 [ L_3_1_1 - 6*L_4_1_1 + L_4_1_2 + L_4_2_1 + L_5_1_1, L_3_2_1 + L_4_1_1 - 6*L_4_2_1 + L_4_2_2 + L_4_3_1 + L_5_2_1, L_3_3_1 + L_4_2_1 - 6*L_4_3_1 + L_4_3_2 + L_4_4_1 + L_5_3_1, L_3_4_1 + L_4_3_1 - 6*L_4_4_1 + L_4_4_2 + L_4_5_1 + L_5_4_1]
889 d(:,:,2) =
890 [           L_1_1_1 - 6*L_1_1_2 + L_1_1_3 + L_1_2_2 + L_2_1_2,           L_1_1_2 + L_1_2_1 - 6*L_1_2_2 + L_1_2_3 + L_1_3_2 + L_2_2_2,           L_1_2_2 + L_1_3_1 - 6*L_1_3_2 + L_1_3_3 + L_1_4_2 + L_2_3_2,           L_1_3_2 + L_1_4_1 - 6*L_1_4_2 + L_1_4_3 + L_1_5_2 + L_2_4_2]
891 [ L_1_1_2 + L_2_1_1 - 6*L_2_1_2 + L_2_1_3 + L_2_2_2 + L_3_1_2, L_1_2_2 + L_2_1_2 + L_2_2_1 - 6*L_2_2_2 + L_2_2_3 + L_2_3_2 + L_3_2_2, L_1_3_2 + L_2_2_2 + L_2_3_1 - 6*L_2_3_2 + L_2_3_3 + L_2_4_2 + L_3_3_2, L_1_4_2 + L_2_3_2 + L_2_4_1 - 6*L_2_4_2 + L_2_4_3 + L_2_5_2 + L_3_4_2]
892 [ L_2_1_2 + L_3_1_1 - 6*L_3_1_2 + L_3_1_3 + L_3_2_2 + L_4_1_2, L_2_2_2 + L_3_1_2 + L_3_2_1 - 6*L_3_2_2 + L_3_2_3 + L_3_3_2 + L_4_2_2, L_2_3_2 + L_3_2_2 + L_3_3_1 - 6*L_3_3_2 + L_3_3_3 + L_3_4_2 + L_4_3_2, L_2_4_2 + L_3_3_2 + L_3_4_1 - 6*L_3_4_2 + L_3_4_3 + L_3_5_2 + L_4_4_2]
893 [ L_3_1_2 + L_4_1_1 - 6*L_4_1_2 + L_4_1_3 + L_4_2_2 + L_5_1_2, L_3_2_2 + L_4_1_2 + L_4_2_1 - 6*L_4_2_2 + L_4_2_3 + L_4_3_2 + L_5_2_2, L_3_3_2 + L_4_2_2 + L_4_3_1 - 6*L_4_3_2 + L_4_3_3 + L_4_4_2 + L_5_3_2, L_3_4_2 + L_4_3_2 + L_4_4_1 - 6*L_4_4_2 + L_4_4_3 + L_4_5_2 + L_5_4_2]
894 d(:,:,3) =
895 [           L_1_1_2 - 6*L_1_1_3 + L_1_1_4 + L_1_2_3 + L_2_1_3,           L_1_1_3 + L_1_2_2 - 6*L_1_2_3 + L_1_2_4 + L_1_3_3 + L_2_2_3,           L_1_2_3 + L_1_3_2 - 6*L_1_3_3 + L_1_3_4 + L_1_4_3 + L_2_3_3,           L_1_3_3 + L_1_4_2 - 6*L_1_4_3 + L_1_4_4 + L_1_5_3 + L_2_4_3]
896 [ L_1_1_3 + L_2_1_2 - 6*L_2_1_3 + L_2_1_4 + L_2_2_3 + L_3_1_3, L_1_2_3 + L_2_1_3 + L_2_2_2 - 6*L_2_2_3 + L_2_2_4 + L_2_3_3 + L_3_2_3, L_1_3_3 + L_2_2_3 + L_2_3_2 - 6*L_2_3_3 + L_2_3_4 + L_2_4_3 + L_3_3_3, L_1_4_3 + L_2_3_3 + L_2_4_2 - 6*L_2_4_3 + L_2_4_4 + L_2_5_3 + L_3_4_3]
897 [ L_2_1_3 + L_3_1_2 - 6*L_3_1_3 + L_3_1_4 + L_3_2_3 + L_4_1_3, L_2_2_3 + L_3_1_3 + L_3_2_2 - 6*L_3_2_3 + L_3_2_4 + L_3_3_3 + L_4_2_3, L_2_3_3 + L_3_2_3 + L_3_3_2 - 6*L_3_3_3 + L_3_3_4 + L_3_4_3 + L_4_3_3, L_2_4_3 + L_3_3_3 + L_3_4_2 - 6*L_3_4_3 + L_3_4_4 + L_3_5_3 + L_4_4_3]
898 [ L_3_1_3 + L_4_1_2 - 6*L_4_1_3 + L_4_1_4 + L_4_2_3 + L_5_1_3, L_3_2_3 + L_4_1_3 + L_4_2_2 - 6*L_4_2_3 + L_4_2_4 + L_4_3_3 + L_5_2_3, L_3_3_3 + L_4_2_3 + L_4_3_2 - 6*L_4_3_3 + L_4_3_4 + L_4_4_3 + L_5_3_3, L_3_4_3 + L_4_3_3 + L_4_4_2 - 6*L_4_4_3 + L_4_4_4 + L_4_5_3 + L_5_4_3]
899 d(:,:,4) =
900 [           L_1_1_3 - 6*L_1_1_4 + L_1_1_5 + L_1_2_4 + L_2_1_4,           L_1_1_4 + L_1_2_3 - 6*L_1_2_4 + L_1_2_5 + L_1_3_4 + L_2_2_4,           L_1_2_4 + L_1_3_3 - 6*L_1_3_4 + L_1_3_5 + L_1_4_4 + L_2_3_4,           L_1_3_4 + L_1_4_3 - 6*L_1_4_4 + L_1_4_5 + L_1_5_4 + L_2_4_4]
901 [ L_1_1_4 + L_2_1_3 - 6*L_2_1_4 + L_2_1_5 + L_2_2_4 + L_3_1_4, L_1_2_4 + L_2_1_4 + L_2_2_3 - 6*L_2_2_4 + L_2_2_5 + L_2_3_4 + L_3_2_4, L_1_3_4 + L_2_2_4 + L_2_3_3 - 6*L_2_3_4 + L_2_3_5 + L_2_4_4 + L_3_3_4, L_1_4_4 + L_2_3_4 + L_2_4_3 - 6*L_2_4_4 + L_2_4_5 + L_2_5_4 + L_3_4_4]
902 [ L_2_1_4 + L_3_1_3 - 6*L_3_1_4 + L_3_1_5 + L_3_2_4 + L_4_1_4, L_2_2_4 + L_3_1_4 + L_3_2_3 - 6*L_3_2_4 + L_3_2_5 + L_3_3_4 + L_4_2_4, L_2_3_4 + L_3_2_4 + L_3_3_3 - 6*L_3_3_4 + L_3_3_5 + L_3_4_4 + L_4_3_4, L_2_4_4 + L_3_3_4 + L_3_4_3 - 6*L_3_4_4 + L_3_4_5 + L_3_5_4 + L_4_4_4]
903 [ L_3_1_4 + L_4_1_3 - 6*L_4_1_4 + L_4_1_5 + L_4_2_4 + L_5_1_4, L_3_2_4 + L_4_1_4 + L_4_2_3 - 6*L_4_2_4 + L_4_2_5 + L_4_3_4 + L_5_2_4, L_3_3_4 + L_4_2_4 + L_4_3_3 - 6*L_4_3_4 + L_4_3_5 + L_4_4_4 + L_5_3_4, L_3_4_4 + L_4_3_4 + L_4_4_3 - 6*L_4_4_4 + L_4_4_5 + L_4_5_4 + L_5_4_4]
904 d(3,3,3)
905 ans =
906 L_2_3_3 + L_3_2_3 + L_3_3_2 - 6*L_3_3_3 + L_3_3_4 + L_3_4_3 + L_4_3_3
907 d(3,3,1)
908 ans =
909 L_2_3_1 + L_3_2_1 - 6*L_3_3_1 + L_3_3_2 + L_3_4_1 + L_4_3_1
910 u=grad3z(L,[1,1,1])
911 u =
912   1×3 <a href="matlab:helpPopup cell" style="font-weight:bold">cell</a> array
913     {5×4×4 sym}    {4×5×4 sym}    {4×4×5 sym}
914 u{3}(:,:,1)
915 ans =
916 [ L_1_1_1, L_1_2_1, L_1_3_1, L_1_4_1]
917 [ L_2_1_1, L_2_2_1, L_2_3_1, L_2_4_1]
918 [ L_3_1_1, L_3_2_1, L_3_3_1, L_3_4_1]
919 [ L_4_1_1, L_4_2_1, L_4_3_1, L_4_4_1]
920 u=grad3z(L,[1,1,1],0)
921 u =
922   1×3 <a href="matlab:helpPopup cell" style="font-weight:bold">cell</a> array
923     {5×4×4 sym}    {4×5×4 sym}    {4×4×5 sym}
925 X =
926   1×3 <a href="matlab:helpPopup cell" style="font-weight:bold">cell</a> array
927     {6×6×6 double}    {6×6×6 double}    {6×6×6 double}
929 L(:,:,1) =
930 [ L_1_1_1, L_1_2_1, L_1_3_1, L_1_4_1, L_1_5_1]
931 [ L_2_1_1, L_2_2_1, L_2_3_1, L_2_4_1, L_2_5_1]
932 [ L_3_1_1, L_3_2_1, L_3_3_1, L_3_4_1, L_3_5_1]
933 [ L_4_1_1, L_4_2_1, L_4_3_1, L_4_4_1, L_4_5_1]
934 [ L_5_1_1, L_5_2_1, L_5_3_1, L_5_4_1, L_5_5_1]
935 L(:,:,2) =
936 [ L_1_1_2, L_1_2_2, L_1_3_2, L_1_4_2, L_1_5_2]
937 [ L_2_1_2, L_2_2_2, L_2_3_2, L_2_4_2, L_2_5_2]
938 [ L_3_1_2, L_3_2_2, L_3_3_2, L_3_4_2, L_3_5_2]
939 [ L_4_1_2, L_4_2_2, L_4_3_2, L_4_4_2, L_4_5_2]
940 [ L_5_1_2, L_5_2_2, L_5_3_2, L_5_4_2, L_5_5_2]
941 L(:,:,3) =
942 [ L_1_1_3, L_1_2_3, L_1_3_3, L_1_4_3, L_1_5_3]
943 [ L_2_1_3, L_2_2_3, L_2_3_3, L_2_4_3, L_2_5_3]
944 [ L_3_1_3, L_3_2_3, L_3_3_3, L_3_4_3, L_3_5_3]
945 [ L_4_1_3, L_4_2_3, L_4_3_3, L_4_4_3, L_4_5_3]
946 [ L_5_1_3, L_5_2_3, L_5_3_3, L_5_4_3, L_5_5_3]
947 L(:,:,4) =
948 [ L_1_1_4, L_1_2_4, L_1_3_4, L_1_4_4, L_1_5_4]
949 [ L_2_1_4, L_2_2_4, L_2_3_4, L_2_4_4, L_2_5_4]
950 [ L_3_1_4, L_3_2_4, L_3_3_4, L_3_4_4, L_3_5_4]
951 [ L_4_1_4, L_4_2_4, L_4_3_4, L_4_4_4, L_4_5_4]
952 [ L_5_1_4, L_5_2_4, L_5_3_4, L_5_4_4, L_5_5_4]
953 L(:,:,5) =
954 [ L_1_1_5, L_1_2_5, L_1_3_5, L_1_4_5, L_1_5_5]
955 [ L_2_1_5, L_2_2_5, L_2_3_5, L_2_4_5, L_2_5_5]
956 [ L_3_1_5, L_3_2_5, L_3_3_5, L_3_4_5, L_3_5_5]
957 [ L_4_1_5, L_4_2_5, L_4_3_5, L_4_4_5, L_4_5_5]
958 [ L_5_1_5, L_5_2_5, L_5_3_5, L_5_4_5, L_5_5_5]
959 size(L)
960 ans =
961      5     5     5
962 diary off