• API
• FAQ
• Tools
• Archive
SHARE
TWEET # Untitled a guest Dec 7th, 2019 137 Never
Not a member of Pastebin yet? Sign Up, it unlocks many cool features!
1. i = n+1
2. h = i
3. k = 0
4. z = 2
5.
6. for j in range(1,n):
7.     while i <= h + n - k-1:
8.         if i == h:
9.             if i == node_number-1:
10.                     y[i] = 0
11.                     theta[i] = 0
12.                     M[i] = Me
13.                     nu[i] = PM(M[i], gamma)
14.                     mu[i] = m.asin(1 / M[i]) * RTOD
15.                     mI= np.tan(m.radians(mu[i - n + k] - theta[i - n + k] + mu[i]) / 2)
16.             else:
17.                     y[i] = 0
18.                     theta[i] = 0
19.                     nu[i] = theta[i - n + k] + nu[i - n + k] + (1 / (np.sqrt((M[i - n + k]**2) - 1) - mp.cot(m.radians(theta[i - n + k]))))
20.                     M[i] = InversePrandtlMeyer(nu[i])
21.                     mu[i] = m.asin(1 / M[i]) *RTOD
22.                     mI= np.tan(m.radians(mu[i - n + k] - theta[i - n + k] + mu[i]) / 2)
23.                     print(mI)
24.         i = i + 1
25.     k = k + 1
26.     h = i
27.
28. 2.184370031112393
29. 1.7380673494238976
30. 1.5419781887088249
31. 1.4039585635698797
32. 1.2989038769190393
33. 1.214861123121338
34. 1.1452669347438043
35. 1.0861568995207378
36. 1.0349662574998555
37. 0.98994725077706
38. 0.9498593079129821
39. 0.9137924304933154
40. 0.8810607737523871
41. 0.8511355525702466
42. 0.8236011144240196
43. 0.7981252463827491
44. 0.7744385459289631
45. 0.752319743774369
46. 0.7315850413556333
47. 0.7120802207977384
48. 0.6936747098613351
49. 0.6762570513386782
50. 0.6597313984083366
51. 0.6440147708625135
52. 0.6290348834011446
53. 0.6147284094475984
54. 0.6010395803489821
55. 0.5879190455792305
56. 0.5753229380400084
RAW Paste Data
We use cookies for various purposes including analytics. By continuing to use Pastebin, you agree to our use of cookies as described in the Cookies Policy.
Not a member of Pastebin yet?