# Voglio risolvere (x-2*y)^2+3*y = 4 rispetto ad x (fissato y=k) o ad y (fissato x=k) # Caricando questo file posso, in genere, farlo semplicemente: source("http://macosa.dima.unige.it/r.R") # # Prima cerco una approssimazione "grafica" delle soluzioni: f <- function(x,y) (x-2*y)^2+3*y - 4 Piano(-3,4,-4,3) curva(f,"brown") # Poi procedo risolvendo numericamente le equazioni, stimando prima dal grafico # intervalli per x, o per y, in cui f(x,k), o f(k,y), cambia segno: # (ricordo che soluz(F,7, 1,3) risolve F(x)=7 se tra 1 e 3 F scavalca il valore 7) # Per quale x f(x,y) passa per y=0 a sinsitra dell'asse y? h <- function(x) f(x,0); soluz(h,0, -3,-1) # -2 # Per quale x passa per y=1, vicino all'asse y? h <- function(x) f(x,1); ; soluz(h,0, 0,2) # 1 # Per quale y f(x,y) passa per x=0 sotto all'asse x? h <- function(y) f(0,y); soluz(h,0, -2,-1) # -1.443 # # Infine controllo le soluzioni graficamente: PUNTO(0,-1.443,"blue"); PUNTO(1,1,"red"); PUNTO(-2,0,"black")