SAGE er et gratis alternativ til dyre og kommersielle MATLAB. Syntaksen er lik den i Python, og funksjonaliteten er (for våre formal) omtrent lik. SAGE kan dessuten kjøres online på http://www.sagenb.org. Alt du trenger er en Google-konto for å logge inn.
Her skal jeg vise hvordan man kan plotte vektorfelt. SAGE har egne funksjoner for å løse diff.-likninger, men disse har jeg ikke satt meg så mye inn i.
Under er kode for å plotte vektorfeltet i Oblig 2.
x,y = var('x,y')
t = var('t')
p =plot_vector_field((x-5*y, 2*x+3*y), (x,-9,5), (y,-10,10))
p.show()
Resultatet ser slik ut:
Pilene er noe tykkere enn de man automatisk får i MATLAB. For å plotte en integralkurve bruker man denne koden:
y11 = (-cos(3*t)-3*sin(3*t))*e^(2*t) y12 = (2*cos(3*t))*e^(2*t) q = parametric_plot((y11,y12), (t,0,0.9)) p+q
Resultatet ser slik ut:
...