function[] = plot_format()
    //Get the handle of current axes
    g = gca()
    //Give labels and set label properties
    g.labels_font_color=5
    g.font_size=3
    g.grid=[1,1]
    g.box="off"
endfunction

//Obtain path of solution file
path = get_absolute_file_path('solution7_18.sce')
//Obtain path of data file
datapath = path + filesep() + 'data7_18.sci'
//Obtain path of function file
funcpath = path + filesep() + 'functions7_18.sci'
//Clear all
clc
//Execute the data file
exec(datapath)
exec(funcpath,[-1])
//Calculate the endurance limit stress for bolt Sdash (N/mm2)
Sdash = (50/100)*Sut
//Calculate Ka, Kb and Kc
[Ka, Kb, Kc] = fluctuate(op, d, reliability)
//Calculate Kd 
Kd = 1/Kf
//Calculate the corrected endurance limit stress Se (N/mm2)
Se = Kb * Kc * Kd * Sdash
//Calculate the stiffness of bolt kb (N/mm)
kb = ((%pi/4)*(d^2))*(E1/l)
//Calculate the area of the two plates Ac (mm2)
Ac = (%pi/4)*(((2*d)^2) - (d^2))
//Calculate the stiffness of the plates kc (N/mm)
kc = (Ac * E2)/l
//Calculate the maximum force in the bolt PMAX (N)
PMAX = (Pi * 1000) + ((kb/(kb + kc))*(Pmax * 1000))
//Calculate the minimum force in the bolt PMIN (N)
PMIN = (Pi * 1000) + ((kb/(kb + kc))*(Pmin * 1000))
//Calculate the mean force and force amplitude
Pm = (PMAX + PMIN)/2
Pa = (PMAX - PMIN)/2
//Plot modified Goodman diagram
//The common quadrilateral in the plot is the area of concern
y1 = {Se 0}
x1 = {0 Sut}
y2 = {Syt 0}
x2 = {0 Syt}
plot(x1,y1,'--*')
plot(x2,y2,'-*')
plot_format()
title('Modified Goodman diagram (Example 7.18)')
xlabel('sigmaM (N/mm2)')
ylabel('sigmaA (N/mm2)')
//Calculate the actual core cross-section area of the bolt A (mm2)
A = (Pa + ((Pi * 1000)/((1 + (Sut/Se)) * fs)))/(Sut/((1 + (Sut/Se))*fs))
//Choose proper diameter from Table 7.1
//Print results
printf('\nCore cross-section area of the bolt(A) = %f mm2\n',A)