var('t') p0=1 phi = pi/6 r= 5 p(t) = p0*exp(sin(phi) + sin(r*t-phi)) C = plot(p, (t,0,6*pi), color = "blue") show(C)
var('L,R') f(L,R) = L^(0.02)*R^(0.08)/(exp(0.00002*L + 0.001*R)) f(1000, 400).n()
Courbe = implicit_plot(f(L,R) == f(1000, 40), (R,0,160),(L,0,4000), color="blue") show(Courbe, aspect_ratio=1/50)