inclination.ks

//Adds a node to change inclination at equatorial ascending node.
declare parameter desired_i.
clearscreen.
function time_to_node {

	local w is ship:orbit:period/360.//calculates ship's speed in seconds per degree.
	local shiptolan is 360 - (orbit:argumentofperiapsis + orbit:trueanomaly).//returns number of degrees to lan.
	if shiptolan < 0 {
		set shiptolan to shiptolan + 360.
	}

	return shiptolan * w.
}
function dV {

	local delta_i is desired_i - orbit:inclination.
	local v is velocityat(ship,time:seconds + time_to_node()):orbit:mag.

	return 2 * v * sin(delta_i/2).
}
set mnv to node(time:seconds + time_to_node(),0,dV(),0).
add mnv.