23/07/31 14:37:34.37 xgtz5I6M0.net
E=c(24+16/60+59/3600,153+59/60+12/3600) ; E
W=c(24+27/60+5/3600,122+55/60+57/3600) ; W
S=c(20+25/60+31/3600,136+4/11) ; S
N=c(45+33/60+26/3600,148+45/60+8/3600) ; N
# uE+(1-u)W == vS+(1-v)N
f = \(u,v) ((u*E+(1-u)*W) - (v*S+(1-v)*N))^2 |> sum()
opt=optim(runif(2),\(x) f(x[1],x[2]),method='L') ; opt
u0=opt$par[1]
v0=opt$par[2]
u0*E+(1-u0)*W
v0*S+(1-v0)*N