# Usando source("http://macosa.dima.unige.it/r.R")
# primo modo:
source("http://macosa.dima.unige.it/r.R")
PLANE(-11,11, -11,11)
circle(0,0, 10, "blue")
# secondo modo:
PLANE(-11,11, -11,11)
f1 <- function(x) sqrt(100-x^2); f2 <- function(x) -f1(x)
graph(f1,-10,10, "blue"); graph(f2,-10,10, "blue")
# terzo modo:
PLANE(-11,11, -11,11)
f <- function(x,y) x^2+y^2-100
CURVE(f, "blue")
#
# Ovvero, senza source(...):
# primo modo:
plot(c(-11,11),c(-11,11),type="n",xlab="",ylab="",asp=1)
axis(1,pos=0,label=FALSE,col="blue"); axis(2,pos=0,label=FALSE,col="blue")
symbols(0, 0, circles=10, inches=FALSE, add=TRUE)
# secondo modo:
f1 <- function(x) sqrt(100-x^2)
plot(c(-11,11),c(-11,11),type="n",xlab="",ylab="",asp=1)
axis(1,pos=0,label=FALSE,col="blue"); axis(2,pos=0,label=FALSE,col="blue")
curve(f1(x),-10,10,add=TRUE); curve(-f1(x),-10,10,add=TRUE)
# terzo modo:
a <- -11; b <- 11; h <- b-a; e <- 1/2
f <- function(x,y) x^2+y^2-100
plot(c(a,b),c(a,b),type="n",xlab="",ylab="",asp=1)
axis(1,pos=0,label=FALSE,col="blue"); axis(2,pos=0,label=FALSE,col="blue")
n <- 500; for(i in 0:n) for(j in 0:n) {
x <- a+h/n*i; y <- a+h/n*j; if(abs(f(x,y)) < e) points(x,y,pch=".")}