It looks like you're new here. If you want to get involved, click one of these buttons!

- All Categories 139.8K
- Programming Languages 104K
- Assembler Developer 6.2K
- Basic 1.8K
- C and C++ 39.7K
- C# 4.2K
- Delphi and Kylix 7.9K
- Haskell 3
- Java 9.5K
- Pascal 4.1K
- Perl 1.3K
- PHP 1.9K
- Python 500
- Ruby 47
- VB.NET 4.3K
- VBA 1.6K
- Visual Basic 20.8K
- Game programming 2.6K
- Console programming 308
- DirectX Game dev 87
- Minecraft 1
- Newbie Game Programmers 107
- Oculus Rift 2
- Applications 8.8K
- Computer Graphics 1.8K
- Computer Hardware 721
- Database & SQL 3.4K
- Electronics development 517
- Matlab 1.6K
- Sound & Music 625
- XML Development 253
- Classifieds 3.2K
- Co-operative Projects 184
- For sale 161
- FreeLance Software City 188
- Jobs Available 1.9K
- Jobs Wanted 597
- Wanted 192
- Microsoft .NET 2.8K
- ASP.NET 1.7K
- .NET General 1.1K
- Miscellaneous 2.8K
- Join the Team 2
- User Profiles 1
- Comments on this site 351
- Computer Emulators 54
- General programming 1.7K
- New programming languages 114
- Off topic board 591
- Mobile & Wireless 146
- Android 20
- Palm Pilot 124
- Multimedia 334
- Demo programming 151
- MP3 programming 183
- Bash scripts 0
- Cloud Computing 10
- FreeBSD 52
- LINUX programming 1.7K
- MS-DOS 361
- Shell scripting 0
- Windows CE & Pocket PC 317
- Windows programming 4.1K
- Software Development 870
- Algorithms 400
- Object Orientation 67
- Project Management 80
- Quality & Testing 87
- Security 230
- WEB-Development 7.3K
- Active Server Pages 1.8K
- AJAX 61
- Bootstrap Themes 1
- CGI Development 55
- ColdFusion 19
- Flash development 221
- HTML & WEB-Design 1.4K
- Internet Development 1.3K
- JavaScript 2.2K
- JQuery 33
- WEB Servers 275
- WEB-Services / SOAP 100

We have migrated to a new platform! Please note that you will need to reset your password to log in (your credentials are still in-tact though). Please contact lee@programmersheaven.com if you have questions.

Tweets by @pheaven
Welcome to the new platform of Programmer's Heaven! We apologize for the inconvenience caused, if you visited us from a broken link of the previous version. The main reason to move to a new platform is to provide more effective and collaborative experience to you all. Please feel free to experience the new platform and use its exciting features. Contact us for any issue that you need to get clarified. We are more than happy to help you.

lapachacha
Posts: **3**Member

hello people...i have beeen given this projectile problem to be solved numerically in FORTRAN..i have coded it using runge0kutta 4th order ......now the thing is that i have not been given the end points of time,,;ie. the range.. and i am asked to

a)the maximum height attained by the projectile

what is to b added to the existing code to achive this?

b)the time required to reach the maximum height

what is to b added to the existing code to achive this?

c)the time required to return to the original elevation..

what is to b added to the existing code to achive this?

i did fortran 4 years ago...havent beeen in practice..this was the lone effort i cud recall...please help...

i have attached the question and the code..

heres the code

CODE

!************************************************* *********************

! *

! projectile question *

! *

!************************************************* **********************

DATA M,G,C/10.0,9.80665,0.1/

!FIRST i HAVE REDUCED THE GIVEN EQUATION IN TWO LINEAR 1ST-ORDER EQUATIONS

!FOR THE SOLUTION TO PROCEED

F1(T,U1,Y2) = U1

F2(T,U1,Y2) = (-M*G-C*U1**2)/M

!ASSIGNMENT OF VALUES TO CONSTANTS

WRITE(*,*) 'Input left and right endpoints separated by'

WRITE(*,*) 'blank'

WRITE(*,*) ' '

READ(*,*) A, B

WRITE(*,*) 'Input the two initial conditions.'

WRITE(*,*) ' '

READ(*,*) ALPHA1, ALPHA2

WRITE(6,*) 'Input a positive integer for the number'

WRITE(6,*) 'of subintervals '

WRITE(6,*) ' '

READ(5,*) N

WRITE(*,6)

6 FORMAT(12X,'t(i)',11X,'w1(i)',11X,'w2(i)')

H=(B-A)/N

T=A

! the initiaal conditions

W1=ALPHA1

W2=ALPHA2

WRITE(*,1) T,W1,W2

!RK PARAMETER EVALUATIONS HERE

DO 110 I=1,N

X11=H*F1(T,W1,W2)

X12=H*F2(T,W1,W2)

X21=H*F1(T+H/2,W1+X11/2,W2+X12/2)

X22=H*F2(T+H/2,W1+X11/2,W2+X12/2)

X31=H*F1(T+H/2,W1+X21/2,W2+X22/2)

X32=H*F2(T+H/2,W1+X21/2,W2+X22/2)

X41=H*F1(T+H,W1+X31,W2+X32)

X42=H*F2(T+H,W1+X31,W2+X32)

W1=W1+(X11+2*X21+2*X31+X41)/6

W2=W2+(X12+2*X22+2*X32+X42)/6

T=A+I*H

WRITE(*,1) T,W1,W2

110 CONTINUE

STOP

1 FORMAT(3(1X,E15.8))

END

a)the maximum height attained by the projectile

what is to b added to the existing code to achive this?

b)the time required to reach the maximum height

what is to b added to the existing code to achive this?

c)the time required to return to the original elevation..

what is to b added to the existing code to achive this?

i did fortran 4 years ago...havent beeen in practice..this was the lone effort i cud recall...please help...

i have attached the question and the code..

heres the code

CODE

!************************************************* *********************

! *

! projectile question *

! *

!************************************************* **********************

DATA M,G,C/10.0,9.80665,0.1/

!FIRST i HAVE REDUCED THE GIVEN EQUATION IN TWO LINEAR 1ST-ORDER EQUATIONS

!FOR THE SOLUTION TO PROCEED

F1(T,U1,Y2) = U1

F2(T,U1,Y2) = (-M*G-C*U1**2)/M

!ASSIGNMENT OF VALUES TO CONSTANTS

WRITE(*,*) 'Input left and right endpoints separated by'

WRITE(*,*) 'blank'

WRITE(*,*) ' '

READ(*,*) A, B

WRITE(*,*) 'Input the two initial conditions.'

WRITE(*,*) ' '

READ(*,*) ALPHA1, ALPHA2

WRITE(6,*) 'Input a positive integer for the number'

WRITE(6,*) 'of subintervals '

WRITE(6,*) ' '

READ(5,*) N

WRITE(*,6)

6 FORMAT(12X,'t(i)',11X,'w1(i)',11X,'w2(i)')

H=(B-A)/N

T=A

! the initiaal conditions

W1=ALPHA1

W2=ALPHA2

WRITE(*,1) T,W1,W2

!RK PARAMETER EVALUATIONS HERE

DO 110 I=1,N

X11=H*F1(T,W1,W2)

X12=H*F2(T,W1,W2)

X21=H*F1(T+H/2,W1+X11/2,W2+X12/2)

X22=H*F2(T+H/2,W1+X11/2,W2+X12/2)

X31=H*F1(T+H/2,W1+X21/2,W2+X22/2)

X32=H*F2(T+H/2,W1+X21/2,W2+X22/2)

X41=H*F1(T+H,W1+X31,W2+X32)

X42=H*F2(T+H,W1+X31,W2+X32)

W1=W1+(X11+2*X21+2*X31+X41)/6

W2=W2+(X12+2*X22+2*X32+X42)/6

T=A+I*H

WRITE(*,1) T,W1,W2

110 CONTINUE

STOP

1 FORMAT(3(1X,E15.8))

END

About & Contact / Terms of use / Privacy statement / Publisher: Lars Hagelin

Programmers Heaven articles / Programmers Heaven files / Programmers Heaven uploaded content / Programmers Heaven C Sharp ebook / Operated by CommunityHeaven LLC

© 1997-2013 Programmersheaven.com - All rights reserved.

## Comments

59Member[code] F1(T,U1,Y2) = U1

F2(T,U1,Y2) = (-M*G-C*U1**2)/M

[/code]

2nd, all of the other FORTRAN looks correct, but I suggest another computer language to solve the problem, if you know any, like, C, MatLab, or Pascal.

- Spam

0 · Vote Down Vote Up ·