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:
...