Hvordan plotte i Sage

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:

...

Publisert 2. nov. 2013 18:36 - Sist endret 23. feb. 2023 14:36