# 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")