blob: 1820e90761243d1a8329ab0a1ab165bb19780893 (
plain)
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
|
// Example 9-7
// Finding e to the power At using laplace transforms
clear; clc;
xdel(winsid()); //close all windows
// please edit the path
// cd "<path for the dependencies>";
// exec("pf_residu.sci");
// exec("ilaplace.sci");
s = %s;
A = [0 1; 0 -2];
L = inv(s*eye(2,2) - A);
disp(L, 'inv(sI - A) =');
// Find the Inverse Laplace transform
for i = 1:2
for j = 1:2
phi(i,j) = ilaplace(L(i,j));
end;
end;
disp(phi,'e^At =');
|
