# Proiezione di z sulla retta z0+tw. Proi = function (z,z0,w) {t=Scalare(z-z0,w)/Scalare(w,w) z0+t*w} # Riflessione di z nel punto m. Riflpunto = function (z,m) {2*m-z} # Riflessione di z nella retta z0+tw. Riflretta = function (z,z0,w) {Riflpunto(z,Proi(z,z0,w))} # Rotazione di z per alfa (in gradi) # attorno al centro. Rot = function (z,alfa,centro=0i) {centro+exp(1i*alfa*CMPid180)*(z-centro)} # Prodotto scalare di punti del piano. Scalare = function (z,w) {Re(z*Conj(w))}