Ringleb_geom_2D.m File Reference

Variables
syms gama q	k
	gama
	x

Variable Documentation

gama

Initial value:

= 1.4

a(gama, q, k) = sqrt( 1 - (gama - 1) * q * q / 2  )

r(gama, q, k) = a^(2/(gama-1))

p(gama, q, k) = 1/gama * a^(2*gama/(gama-1))

J(gama, q, k) = 1/a + 1/(3*a^3) + 1/(5*a^5) - 1/2*log((1+a)/(1-a))

% --------------------------------------- %

x (gama, q, k) = 1/(2*r) * (2/k^2 - 1/q^2) - J/2

y (gama, q, k) = 1/(k*r*q) * sqrt( 1 - (q/k)^2 )

% -------- gama = 1.4 ------------------- %

x_14 = subs(x, gama, 1.4)
y_14 = subs(y, gama, 1.4)

% -------- inflow --------%
%  q = q_0 = 0.5
%  k_min <= k <= k_max

figure
hold on

syms gama q k

x