Flat truss to solve in pyhton

$E=2{10}^{8}k\frac{N}{m^2}$

I want you to add the following explicitly step by step, indicating the equation that was used to solve in comments. It is to solve the exercise using the stiffness method of a truss in Python and the instructions to solve are a$)$ print on the screen the system of equations in local coordinates of each element using Python.

b$)$ print on the screen the system of equations of each element in global coordinates in python.

c$)$ Print on the screen in python the Calculation of local displacements indicating node $1$ local displacement in x is:

node $1$ local displacement in y is:

d$)$ prints on the screen the displacement in global coordinates of each bar. indicating the $1$ x global node is:

node $1$ and global is:

e$)$ Print the reactions on the screen indicating the restricted degree of freedom in python.

f$)$ Print on the screen the calculation of axial forces in each element per bar in local coordinates indicating the axial forces of bar number $1$ is: $.$

g$)$ Print on the screen the calculation of axial forces in each element per bar in global coordinates indicating the axial forces of bar number $1$ in global x is:

axial force of bar $1$ in y global is:

h$)$ Plot the axial force of the complete structure

import numpy as np

import matplotlib as mpl

import matplotlib.pyplot as plt

'Code to solve $2$D trusses'

'Units are consistent and depend on the user'

#$\%\%$ Geometry of the structure

#Input nodes and coordinates#

#Initializing node names matrix

ID $=[1,2,3,4,5]$

#Enter the coordinates

CoordinatesX $=[0,5,15,5,15]$

CoordinatesY $=[0,8\mathrm{.}5,8\mathrm{.}5,0,0]$

#$\%\%$ Input supports

U$\_$X $=[1,0,0,0,0]$ #####X Restriction

U$\_$Y $=[1,0,0,0,1]$ #####Y Restriction

#$\%\%$ Input connectivity

Node$\_$i $=[1,2,2,2,1,4,3]$

Node$\_$j $=[2,3,5,4,4,5,5]$

#$\%\%$ Input properties

A $=[0\mathrm{.}03,0\mathrm{.}03,0\mathrm{.}03,0\mathrm{.}03,0\mathrm{.}03,0\mathrm{.}03,0\mathrm{.}03]$

E $=[2*10**11,2*10**11,2*10**11,2*10**11,2*10**11,2*10**11,2*10**11]$

#$\%\%$ Input nodal loads

Fx $=[4250,4250,10000,0,0]$

Fy $=[-2500,-3333\mathrm{.}333,-1666\mathrm{.}67,-20000,0]$

#$\%\%$ Scale factor for deformations

esc $=10000$

#$\%\%$

############################ Function definitions ###################

#$\%\%$ Plot function

def Plot$($Node$\_$i$,$ Node$\_$j$,$ ID$,$ CoordinatesX, CoordinatesY$)$:

fig $=$ plt$.$figure$()$

for i in range$($len$($Node$\_$i$))$:

ni $=$ Node$\_$i$[$i$]$

nf $=$ Node$\_$j$[$i$]$

XX $=[]$

YY $=[]$

for j in range$($len$($ID$))$:

if ni $==$ ID$[$j$]$:

XX$.$append$($CoordinatesX$[$j$])$

YY$.$append$($CoordinatesY$[$j$])$

if nf $==$ ID$[$j$]$:

XX$.$append$($CoordinatesX$[$j$])$

YY$.$append$($CoordinatesY$[$j$])$

plt$.$plot$($XX$,$ YY$,$ marker$=\text{'}$o$\text{'},$ color$=$'black'$)$

plt$.$show$()$

#$\%\%$

################################################################

############## Calculations begin ############################

################################################################

#$\%\%$ Plotting the structure

Plot$($Node$\_$i$,$ Node$\_$j$,$ ID$,$ CoordinatesX, CoordinatesY$)$

#$\%\%$ GENERATING THE GLOBAL STIFFNESS MATRIX

# Generate the global DOF $($Degrees of Freedom$)$ vector

n$\_$nodes $=$ len$($CoordinatesX$)$

Indices $=[]$ # Vector of indices ordered by DOF and ID in alphanumeric, used for assembling

for i in range$($n$\_$nodes$)$:

Indices.append$(\text{'}$x$\text{'}+$ str$($ID$[$i$]))$

Indices.append$(\text{'}$y$\text{'}+$ str$($ID$[$i$]))$

print$(\text{'}$The DOFs of the structure are', Indices$)$

Indices$\_$num $=[]$ # Numeric indices vector for each DOF, used to identify rows and columns to remove

for i in range$($n$\_$nodes $*2)$:

Indices$\_$num.append$($i$)$

# Identify the restricted DOFs

Restriction $=[]$

for i in range$($n$\_$nodes$)$:

if U$\_$X$[$i$]==1$:

Restriction.append$(\text{'}$x$\text{'}+$ str$($ID$[$i$]))$

if U$\_$Y$[$i$]==1$:

Restriction.append$(\text{'}$y$\text{'}+$ str$($ID$[$i$]))$

print$(\text{'}$The restricted DOFs are', Restriction$)$

# Create the same restriction vector but in numeric form

Restriction$\_$num $=[]$

for i in range$($len$($Restriction$))$:

for j in range$($n$\_$nodes $*2)$:

if Restriction$[$i$]==$ Indices$[$j$]$:

Restriction$\_$num.append$($Indices$\_$num$[$j$])$

# Calculating the local and assembled matrices

# A loop is created to iterate through the number of elements,

# since local matrices are element$-$specific

L $=[]$ # Initialize the length vector

num$\_$el $=$ len$($Node$\_$i$)$ # The length of connectivity equals the number of elements

# Initialize the global stiffness matrix

K $=$ np$.$asmatrix$($np$.$zeros$([$n$\_$nodes$*2,$ n$\_$nodes$*2]))$

for i in range$($num$\_$el$)$:

# Extracting coordinates

xi $=$ CoordinatesX$[$Node$\_$i$[$i$]-1]$

xj $=$ CoordinatesX$[$Node$\_$j$[$i$]-1]$

yi $=$ CoordinatesY$[$Node$\_$i$[$i$]-1]$

yj $=$ CoordinatesY$[$Node$\_$j$[$i$]-1]$

# Calculating length

L$.$append$((($xi $-$ xj$)**2+($yi $-$ yj$)**2)**0\mathrm{.}5)$

# Generating local DOF alphanumeric indices

IL $=[]$

IL$.$append$(\text{'}$x$\text{'}+$ str$($Node$\_$i$[$i$]))$

IL$.$append$(\text{'}$y$\text{'}+$ str$($Node$\_$i$[$i$]))$

IL$.$append$(\text{'}$x$\text{'}+$ str$($Node$\_$j$[$i$]))$

IL$.$append$(\text{'}$y