Read about fplot/ ezplot.
clc
clear all
close all
syms t
A = [-1, 1, 0;
1,-2, 1;
0, 1, -1];
[P,D] = eig(A)
Phi = P*expm(D*t)*inv(P)
thetanot=[0.1;0.3;0.2]
t_t=Phi*thetanot
dt = 1/1000;
t = 0:dt:1000*dt;
S1=sum(t_t(1,1)) ;
fplot(S1)
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!