octave:2> A=rand(4) A = 0.772248 0.031450 0.680464 0.875687 0.012036 0.280289 0.579630 0.255253 0.480991 0.105398 0.636785 0.536776 0.209713 0.180108 0.587152 0.576423 octave:3> phi1m(A) ans = 1.689263 0.093810 0.761984 0.867341 0.108362 1.185367 0.472916 0.282392 0.458312 0.112443 1.628996 0.574905 0.267512 0.148123 0.555622 1.532848 octave:4> det(A) ans = 0.0025380 octave:5> A(expm(A)-eye(4)) error: subscript indices must be either positive integers less than 2^31 or logicals octave:5> A\(expm(A)-eye(4)) ans = 1.689263 0.093810 0.761984 0.867341 0.108362 1.185367 0.472916 0.282392 0.458312 0.112443 1.628996 0.574905 0.267512 0.148123 0.555622 1.532848 octave:6> phi1m(0) ans = 1 octave:7> z=0.1;(exp(z)-1)/z ans = 1.0517 octave:8> format long e octave:9> z=0.1;(exp(z)-1)/z ans = 1.05170918075648e+00 octave:10> z=1e-8;(exp(z)-1)/z ans = 9.99999993922529e-01 octave:11> z=1e-12;(exp(z)-1)/z ans = 1.00008890058234e+00 octave:12> z=1e-14;(exp(z)-1)/z ans = 9.99200722162641e-01 octave:13> z=1e-16;(exp(z)-1)/z ans = 0.00000000000000e+00 octave:14> z=1e-16;(exp(z)-1)/log(exp(z)) warning: division by zero ans = NaN octave:15> z=1e-14;(exp(z)-1)/log(exp(z)) ans = 1.00000000000001e+00 octave:16> z=1e-14;(exp(z)-1)/z ans = 9.99200722162641e-01 octave:17> ese11ODE octave:18> ese11ODE octave:19> y y = 6.67694573647680e-05 9.44263526544512e-05 6.67694573647680e-05 octave:20> ese11ODE error: 'yExp' undefined near line 12 column 12 error: called from: error: /home/accounts/personale/clrmrc90/aa1617/equazioni_differenziali/ese11ODE.m at line 12, column 10 octave:20> ese11ODE octave:21> yEexp yEexp = 6.67694573647680e-05 9.44263526544512e-05 6.67694573647680e-05 octave:22> yERexp yERexp = 6.61285911419947e-05 9.35200215987930e-05 6.61285911419947e-05 octave:23> ese11ODE octave:24> ese1PDE octave:25> A A = Compressed Column Sparse (rows = 10, cols = 10, nnz = 28 [28%]) (1, 1) -> -6.56561270002349e+01 (2, 1) -> 3.28280635001174e+01 (1, 2) -> 3.28280635001174e+01 (2, 2) -> -6.56561270002349e+01 (3, 2) -> 3.28280635001174e+01 (2, 3) -> 3.28280635001174e+01 (3, 3) -> -6.56561270002349e+01 (4, 3) -> 3.28280635001174e+01 (3, 4) -> 3.28280635001174e+01 (4, 4) -> -6.56561270002349e+01 (5, 4) -> 3.28280635001174e+01 (4, 5) -> 3.28280635001174e+01 (5, 5) -> -6.56561270002349e+01 (6, 5) -> 3.28280635001174e+01 (5, 6) -> 3.28280635001174e+01 (6, 6) -> -6.56561270002349e+01 (7, 6) -> 3.28280635001174e+01 (6, 7) -> 3.28280635001174e+01 (7, 7) -> -6.56561270002349e+01 (8, 7) -> 3.28280635001174e+01 (7, 8) -> 3.28280635001174e+01 (8, 8) -> -6.56561270002349e+01 (9, 8) -> 3.28280635001174e+01 (8, 9) -> 3.28280635001174e+01 (9, 9) -> -6.56561270002349e+01 (10, 9) -> 3.28280635001174e+01 (9, 10) -> 3.28280635001174e+01 (10, 10) -> -6.56561270002349e+01 octave:26> full(A) ans = Columns 1 through 3: -6.56561270002349e+01 3.28280635001174e+01 0.00000000000000e+00 3.28280635001174e+01 -6.56561270002349e+01 3.28280635001174e+01 0.00000000000000e+00 3.28280635001174e+01 -6.56561270002349e+01 0.00000000000000e+00 0.00000000000000e+00 3.28280635001174e+01 0.00000000000000e+00 0.00000000000000e+00 0.00000000000000e+00 0.00000000000000e+00 0.00000000000000e+00 0.00000000000000e+00 0.00000000000000e+00 0.00000000000000e+00 0.00000000000000e+00 0.00000000000000e+00 0.00000000000000e+00 0.00000000000000e+00 0.00000000000000e+00 0.00000000000000e+00 0.00000000000000e+00 0.00000000000000e+00 0.00000000000000e+00 0.00000000000000e+00 Columns 4 through 6: 0.00000000000000e+00 0.00000000000000e+00 0.00000000000000e+00 0.00000000000000e+00 0.00000000000000e+00 0.00000000000000e+00 3.28280635001174e+01 0.00000000000000e+00 0.00000000000000e+00 -6.56561270002349e+01 3.28280635001174e+01 0.00000000000000e+00 3.28280635001174e+01 -6.56561270002349e+01 3.28280635001174e+01 0.00000000000000e+00 3.28280635001174e+01 -6.56561270002349e+01 0.00000000000000e+00 0.00000000000000e+00 3.28280635001174e+01 0.00000000000000e+00 0.00000000000000e+00 0.00000000000000e+00 0.00000000000000e+00 0.00000000000000e+00 0.00000000000000e+00 0.00000000000000e+00 0.00000000000000e+00 0.00000000000000e+00 Columns 7 through 9: 0.00000000000000e+00 0.00000000000000e+00 0.00000000000000e+00 0.00000000000000e+00 0.00000000000000e+00 0.00000000000000e+00 0.00000000000000e+00 0.00000000000000e+00 0.00000000000000e+00 0.00000000000000e+00 0.00000000000000e+00 0.00000000000000e+00 0.00000000000000e+00 0.00000000000000e+00 0.00000000000000e+00 3.28280635001174e+01 0.00000000000000e+00 0.00000000000000e+00 -6.56561270002349e+01 3.28280635001174e+01 0.00000000000000e+00 3.28280635001174e+01 -6.56561270002349e+01 3.28280635001174e+01 0.00000000000000e+00 3.28280635001174e+01 -6.56561270002349e+01 0.00000000000000e+00 0.00000000000000e+00 3.28280635001174e+01 Column 10: 0.00000000000000e+00 0.00000000000000e+00 0.00000000000000e+00 0.00000000000000e+00 0.00000000000000e+00 0.00000000000000e+00 0.00000000000000e+00 0.00000000000000e+00 3.28280635001174e+01 -6.56561270002349e+01 octave:27> format octave:28> full(A) ans = Columns 1 through 7: -65.65613 32.82806 0.00000 0.00000 0.00000 0.00000 0.00000 32.82806 -65.65613 32.82806 0.00000 0.00000 0.00000 0.00000 0.00000 32.82806 -65.65613 32.82806 0.00000 0.00000 0.00000 0.00000 0.00000 32.82806 -65.65613 32.82806 0.00000 0.00000 0.00000 0.00000 0.00000 32.82806 -65.65613 32.82806 0.00000 0.00000 0.00000 0.00000 0.00000 32.82806 -65.65613 32.82806 0.00000 0.00000 0.00000 0.00000 0.00000 32.82806 -65.65613 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 32.82806 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 Columns 8 through 10: 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 32.82806 0.00000 0.00000 -65.65613 32.82806 0.00000 32.82806 -65.65613 32.82806 0.00000 32.82806 -65.65613 octave:29> ese1PDE error: ese1PDE: operator +: nonconformant arguments (op1 is 100x100, op2 is 100x1) error: called from: error: /home/accounts/personale/clrmrc90/aa1617/equazioni_differenziali/ese1PDE.m at line 18, column 5 octave:29> ese1PDE octave:30> ese1PDE octave:31> ese1PDE octave:32> ese1PDE octave:33> ese1PDE octave:34> ese1PDE octave:35> y(1) ans = 5.0064e-14 octave:36> ese1PDE octave:37> ese1PDE octave:38> ese1PDE octave:39> ese1PDE octave:40> ese1PDE octave:41> ese1PDE octave:42> ese1PDE octave:43> ese1PDE octave:44> ese1PDE octave:45> ese1PDE octave:46> ese1PDE octave:47> quit