Merge branch 'code' into main, just to not loose track of it.

Code/2Firm_vfi.jl

@@ -56,3 +56,8 @@ With methods which show own and other sizes as vectors?
 function loss(constellation_sizes::Vector{<:Int},debris::Float64,owner::Int)
     1-exp(-α1*sum(constellation_sizes) + (α1-α2)*constellation_sizes[owner] - α3*debris)
 end
+
+#=
+I think that I might need to consider the launch quantities  of the other operators
+as a state variable and use a nash equilibrium style idea of equilibrium.
+=#