function [t,Y]=lindys(A,tvali,y0)
% Ratkaistaan lineaarinen DYS matriisiexp-funktiolla
% Syötteet:
%    matriisi A (n x n)
%    aikaväli tvali [a,b]
%    alkuarvovektori  y0 (hetkellä t=0)
% Tulos: T:aikapisteet 100 kpl.
%        Y (100,n) : ratkaisusarakkeet
% Esim: A=[1 0 0;1 3 0;0 1 1]; y0=[1;2;3];
%       [T,Y]=lindys(A,[-2 5],y0);
%       plot(T,Y)  % aikaparvi
%        plot(Y(:,1),Y(:,2)  % Faasitaso (1 2)
%        plot3(Y(:,1),Y(:,2),Y(:,3)); grid on % Faasiavaruus

t=linspace(tvali(1),tvali(2));
[n,n]=size(A);
Y=zeros(n,100);
for j=1:100
    Y(:,j)=expm(A*t(j))*y0;
end;
Y=Y';

Input argument "tvali" is undefined.

Error in ==> lindys at 15
t=linspace(tvali(1),tvali(2));

Published with MATLAB® 7.6