Docsity
Docsity
Prepare for your exams
Prepare for your exams

Study with the several resources on Docsity

Earn points to download
Earn points to download

Earn points by helping other students or get them with a premium plan

Guidelines and tips
Guidelines and tips
Sell on Docsity
Sell on Docsity
Docsity AI
Docsity AI
Log inSign up

Prepare for your exams

Study with the several resources on Docsity

Find documents

Prepare for your exams with the study notes shared by other students like you on Docsity

Search for your university

Find the specific documents for your university's exams

Docsity AINEW

Summarize your documents, ask them questions, convert them into quizzes and concept maps

Explore questions

Clear up your doubts by reading the answers to questions asked by your fellow students

Earn points to download

Earn points by helping other students or get them with a premium plan

Share documents
download point icon20 Points

For each uploaded document

Answer questions
download point icon5 Points

For each given answer (max 1 per day)

All the ways to get free points
Get points immediately

Choose a premium plan with all the points you need

Study Opportunities

Choose your next study program

Get in touch with the best universities in the world. Search through thousands of universities and official partners

Community

Ask the community

Ask the community for help and clear up your study doubts

Free resources

Our save-the-student-ebooks!

Download our free guides on studying techniques, anxiety management strategies, and thesis advice from Docsity tutors

Fortran Program: Calculate K11-K42 and W1, W2 with A, B, N, H, TK, ALPHA, X0, Y, U1, U2, Exercises of Computer Fundamentals

Birla Institute of Technology and ScienceComputer Fundamentals

A fortran program that calculates the values of k11 to k42 and w1, w2 using given values of a, b, n, h, tk, alpha, x0, y, and u1, u2. The program uses a do loop to iterate through each value of x and calculates the corresponding k and w values using a series of calculations involving h, w2, x, and the previously calculated k values.

Typology: Exercises

2011/2012

Uploaded on 07/17/2012

pafmavasa
pafmavasa 🇮🇳

4.7

(24)

76 documents

1 / 7

Toggle sidebar

This page cannot be seen from the preview

Don't miss anything!

bg1
SOURCE CODE:
REAL A, B, N, H, TK, ALPHA, X0, Y
REAL K11, K12, K21, K22, K31, K32, K41, K42
REAL dK11, dK12, dK21, dK22, dK31, dK32, dK41, dK42
REAL W1, W2, U1, U2, dW1, dW2
A = 1.0
B = 3.0
N = 10
H = (B - A) / N
WRITE(*,*)' '
WRITE(*,*)' '
WRITE(*,*)'H =', H
WRITE(*,*)' '
WRITE(*,*)' '
TK = -14.000203
ALPHA = 17
W1 = ALPHA
W2 = TK
U1 = 0
U2 = 1
DO I=1, N
WRITE(*,*)'I =', I
WRITE(*,*)' '
WRITE(*,*)' '
dW1 = W1
dW2 = W2
X = A + (I - 1) * H
WRITE(*,*)'X =', X
WRITE(*,*)' '
WRITE(*,*)' '
K11 = H * W2
WRITE(*,*)'K11 =', K11
WRITE(*,*)' '
WRITE(*,*)' '
K12 = H * ((1.0/8.0) * (32.0 + (2.0*X**3) - (W1*W2)))
WRITE(*,*)'K12 =', K12
WRITE(*,*)' '
WRITE(*,*)' '
K21 = H * (W2 + (K12/2.0))
WRITE(*,*)'K21 =', K21
WRITE(*,*)' '
WRITE(*,*)' '
K22 = H * ((1.0/8.0) * (32.0 + (2.0*(X+H/2.0)**3)
$ - ((W1+K11/2) * (W2+K12/2)) ))
WRITE(*,*)'K22=', K22
WRITE(*,*)' '
docsity.com
pf3
pf4
pf5

Partial preview of the text

Download Fortran Program: Calculate K11-K42 and W1, W2 with A, B, N, H, TK, ALPHA, X0, Y, U1, U2 and more Exercises Computer Fundamentals in PDF only on Docsity!

SOURCE CODE:

REAL A, B, N, H, TK, ALPHA, X0, Y

REAL K11, K12, K21, K22, K31, K32, K41, K

REAL dK11, dK12, dK21, dK22, dK31, dK32, dK41, dK REAL W1, W2, U1, U2, dW1, dW A = 1. B = 3. N = 10 H = (B - A) / N WRITE(,)' ' WRITE(,)' ' WRITE(,)'H =', H WRITE(,)' ' WRITE(,)' ' TK = -14. ALPHA = 17 W1 = ALPHA W2 = TK U1 = 0 U2 = 1 DO I=1, N WRITE(,)'I =', I WRITE(,)' ' WRITE(,)' ' dW1 = W dW2 = W X = A + (I - 1) * H WRITE(,)'X =', X WRITE(,)' ' WRITE(,)' ' K11 = H * W WRITE(,)'K11 =', K WRITE(,)' ' WRITE(,)' ' K12 = H * ((1.0/8.0) * (32.0 + (2.0X3) - (W1W2))) WRITE(,)'K12 =', K WRITE(,)' ' WRITE(,)' ' K21 = H * (W2 + (K12/2.0)) WRITE(,)'K21 =', K WRITE(,)' ' WRITE(,)' ' K22 = H * ((1.0/8.0) * (32.0 + (2.0(X+H/2.0)3) $ - ((W1+K11/2) * (W2+K12/2)) )) WRITE(,)'K22=', K WRITE(,*)' '

WRITE(,)' '

K31 = H * (W2 + K22/2.0)

WRITE(,)'K31 =', K

WRITE(,)' '

WRITE(,)' '

K32 = H * ((1.0/8.0) * (32.0 + (2.0(X+H/2.0)*3)

$ - ((W1+K21/2) * (W2+K22/2)) ))

WRITE(,)'K32 =', K

WRITE(,)' '

WRITE(,)' '

K41 = H * (W2 + K32)

WRITE(,)'K41 =', K

WRITE(,)' '

WRITE(,)' '

K42 = H * ((1.0/8.0) * (32.0 + (2.0(X+H)*3)

$ - ((W1+K31) * (W2+K32)) ))

WRITE(,)'K42 =', K

WRITE(,)' '

WRITE(,)' '

W1 = W1 + (K11 + 2K21 + 2K31 + K41)/

WRITE(,)'W1 =', W

WRITE(,)' '

WRITE(,)' '

W2 = W2 + (K12 + 2K22 + 2K32 + K42)/

WRITE(,)'W2 =', W

WRITE(,)' '

WRITE(,)' '

dK11 = H * U WRITE(,)'dK11 =', dK WRITE(,)' ' WRITE(,)' ' dK12 = H * ((U1 * ((-1.0/8.0) * (dW2))) $ + (U2 * ((-1.0/8.0) * (dW1))) ) WRITE(,)'dK12 =', dK WRITE(,)' ' WRITE(,)' ' dK21 = H * (U2 + dK12/2) WRITE(,)'dK21 =', dK WRITE(,)' ' WRITE(,)' ' dK22 = H * (((U1+dK11/2) * ((-1.0/8.0) * (dW2))) $ + ((U2+dK12/2) * ((-1.0/8.0) * (dW1))) ) WRITE(,)'dK22 =', dK WRITE(,)' ' WRITE(,)' ' dK31 = H * (U2 + dK22/2)

Press any key to continue

  • W2 = -8.
  • dK11 = 2.000000E-
  • dK12 = -4.250000E-
  • dK21 = 1.575000E-
  • dK22 = -2.996870E-
  • dK31 = 1.700313E-
  • dK32 = -3.337536E-
  • dK41 = 1.332493E-
  • dK42 = -2.236429E-
  • U1 = 1.647187E-
  • U2 = 6.807460E-
  • Yi = 13.
  • I =
  • X = 1.
  • K11 = -1.
  • K12 = 4.
  • K21 = -1.
  • K22= 3.
  • K31 = -1.
  • K32 = 3.
  • K41 = -1.
  • K42 = 2.
  • W1 = 13.
  • W2 = -5.
  • dK11 = 1.361492E-
  • dK12 = -2.160768E-
  • dK21 = 1.145415E-
  • dK22 = -1.614673E-
  • dK31 = 1.200025E-
  • dK32 = -1.738886E-
  • dK41 = 1.013715E-
  • dK42 = -1.259320E-
  • U1 = 2.824868E-
  • U2 = 5.119592E-
  • Yi = 12.
  • I =
  • X = 1.
  • K11 = -1.
  • K12 = 2.
  • K21 = -7.537050E-
  • K22= 2.
  • K31 = -8.009763E-
  • K32 = 2.
  • K41 = -5.469890E-
  • K42 = 2.
  • W1 = 12.
  • W2 = -2.
  • dK11 = 1.023918E-
  • dK12 = -1.352225E-
  • dK21 = 8.886959E-
  • dK22 = -1.058938E-
  • dK31 = 9.180246E-
  • dK32 = -1.116921E-
  • dK41 = 8.005342E-
  • dK42 = -8.583542E-
  • U1 = 3.731183E-
  • U2 = 4.025876E-
  • Yi = 12.
  • I =
  • X = 2.
  • K11 = -5.531506E-
  • K12 = 2.
  • K21 = -3.457285E-
  • K22= 1.
  • K31 = -3.734006E-
  • K32 = 1.
  • K41 = -1.841270E-
  • K42 = 1.
  • W1 = 12.
  • W2 = -9.367088E-
  • dK11 = 8.051752E-
  • dK12 = -1.014542E-
  • dK21 = 7.037210E-
  • dK22 = -8.263632E-
  • dK31 = 7.225388E-
  • dK32 = -8.596111E-
  • dK41 = 6.332529E-
  • dK42 = -6.928701E-
  • U1 = 4.446341E-
  • U2 = 3.179316E-
  • Yi = 12.
  • H = 2.000000E-
  • I =
  • X = 1.
  • K11 = -2.
  • K12 = 6.
  • K21 = -2.
  • K22= 5.
  • K31 = -2.
  • K32 = 5.
  • K41 = -1.
  • K42 = 4.
  • W1 = 14.
  • W2 = -8.
  • dK11 = 2.000000E-
  • dK12 = -4.250000E-
  • dK21 = 1.575000E-
  • dK22 = -2.996870E-
  • dK31 = 1.700313E-
  • dK32 = -3.337536E-
  • dK41 = 1.332493E-
  • dK42 = -2.236429E-
  • U1 = 1.647187E-
  • U2 = 6.807460E-
  • Yi = 13.
  • I =
  • X = 1.
  • K11 = -1.
  • K12 = 4.
  • K21 = -1.
  • K22= 3.
  • K31 = -1.
  • K32 = 3.
  • K41 = -1.
  • K42 = 2.
  • W1 = 13.
  • W2 = -5.
  • dK11 = 1.361492E-
  • dK12 = -2.160768E-
  • dK21 = 1.145415E-
  • dK22 = -1.614673E-
  • dK31 = 1.200025E-
  • dK32 = -1.738886E-
  • dK41 = 1.013715E-
  • dK42 = -1.259320E-
  • U1 = 2.824868E-
  • U2 = 5.119592E-
  • Yi = 12.
  • I =
  • X = 1.
  • K11 = -1.
  • K12 = 2.
  • K21 = -7.537050E-
  • K22= 2.
  • K31 = -8.009763E-
  • K32 = 2.
  • K41 = -5.469890E-
  • K42 = 2.
  • W1 = 12.
  • W2 = -2.
  • dK11 = 1.023918E-
  • dK12 = -1.352225E-
  • dK21 = 8.886959E-
  • dK22 = -1.058938E-
  • dK31 = 9.180246E-
  • dK32 = -1.116921E-
  • dK41 = 8.005342E-
  • dK42 = -8.583542E-
  • U1 = 3.731183E-
  • U2 = 4.025876E-
  • Yi = 12.
  • I =
  • X = 2.
  • K11 = -5.531506E-
  • K12 = 2.
  • K21 = -3.457285E-
  • K22= 1.
  • K31 = -3.734006E-
  • K32 = 1.
  • K41 = -1.841270E-
  • K42 = 1.
  • W1 = 12.
  • W2 = -9.367088E-
  • dK11 = 8.051752E-
  • dK12 = -1.014542E-
  • dK21 = 7.037210E-
  • dK22 = -8.263632E-
  • dK31 = 7.225388E-
  • dK32 = -8.596111E-
  • dK41 = 6.332529E-
  • dK42 = -6.928701E-
  • U1 = 4.446341E-
  • U2 = 3.179316E-
  • Yi = 12.