source("http://macosa.dima.unige.it/R/curve.txt") source("http://macosa.dima.unige.it/R/curve2.txt") # Se hai caricato curve.txt puoi tracciare curve in forma parametrica # battendo param(X,Y, a,b, col) se hai già definito 2 funz. X e Y di un # param. che varia tra a e b, o in forma polare con polar(R, a,b, col) # se hai def. R funz. dell'angolo tra a e b. Es. param(f,g,0,2*pi,'red') # con f <- function(t) 2*cos(t), g <- function(t) sin(t) è un'ellisse; # polar(r,0,20,'blue') con r <- function(a) sqrt(a) è una spirale. # Puoi anche rototraslare una curva K(x,y)=0 di un angolo di A gradi e di # un vettore (h,k) e dilatarla orizzontalmente e verticalmente di due # fattori mx e my battendo rotramol(K, A, h,k, mx,my, col) # Usando rotra(K, A, h,k, col) puoi solo rototraslarla. boxm(-60,60, -60,60) ro1 <- function(a) a^1.5 polar(ro1,0,5*pi, "red") ro2 <- function(a) 30/(2+sin(a)-cos(a)) polar(ro2,0,5*pi, "green") asc1 <- function(t) 2*t^2-3*t ord1 <- function(t) t^3-t+1 param(asc1, ord1, -10,10, "blue") asc2 <- function(t) -30+20*cos(t) ord2 <- function(t) 10+30*sin(t) param(asc2, ord2, 0,2*pi, "magenta") # # H <- function(x,y) x*y-1 boxm(-7,7, -7,7); curva(H, "blue") rotra(H, -45, -2,-4, "red") rotra(H, atan(1/2)*180/pi, 4,2, "green4") # # boxm(-7,7, -7,7); curva(H, "blue") rotramol(H, 45, 0,0, 1.5,1, "red") #