Register to reply

Need in solving a programm in FOTRAN

by Sampath Emani
Tags: fotran, programm, solving
Share this thread:
Sampath Emani
#1
Feb11-12, 01:04 PM
P: 2
Friends,
i have written a program in FOTRAN for solving a equation [f(x0,y0)= 2*(y**2)+5]
I'm in the learning stage of FOTRAN. So please don't mind if my program has got any silly mistakes.
Please help in executing the program.

program rungekutta


print *,' DIFFERENTIAL EQUATION WITH 1 VARIABLE OF ORDER 1'

print*,"enter the value of x0"
read*,x

print*,"enter the value of y0"
read*,y

print*,"enter the value of step size"
read*,h

print*,"solution for the equation is",rk(x0,y0,h)


call runge(m,f1,n,h,x1)

end

function rk(y) result(r)

!Example: y'=(2*y^2+5)

real y,r

r = (2*(y**2)+5)
end function rk(y)

Subroutine runge(m,f1,n,h,x1)

integer f1,i,j,m,n,h

real x,w1,w2,w3,w4,x0,y;

do i=1, m
ni = (i - 1) * f1 - 1
do j=1, f1
x = x1 + h * (n + j)
w1 = h * f(x0,y0)
w2 = h * f(x0+h/2, y0+w1/2)
w3 = h * f(x0+h/2, y0+w2/2)
w4 = h * f(x0+h, y0+w3)
x = x0 + h

x = y + (w1 + w2 + w2 + w3 + w3 + w4) / 6

end do

end do

end
Phys.Org News Partner Science news on Phys.org
Wildfires and other burns play bigger role in climate change, professor finds
SR Labs research to expose BadUSB next week in Vegas
New study advances 'DNA revolution,' tells butterflies' evolutionary history
gsal
#2
Feb11-12, 04:17 PM
P: 873
what does the subroutine runge do for you? it is not being used anywhere.

why if function rk is declared with 1 argument, you pass 3 to it?
Sampath Emani
#3
Feb11-12, 10:55 PM
P: 2
yes sir,
i do agree. but I'm not able not correct that .
Can u help me in writing the correct step.
I took an equation (2*(y**2)+5)
i want to solve this using Runge kutta.
Thanks in advance

Louisa
#4
Mar15-12, 06:12 AM
P: 24
Need in solving a programm in FOTRAN

Quote Quote by Sampath Emani View Post
Friends,
i have written a program in FOTRAN for solving a equation [f(x0,y0)= 2*(y**2)+5]
I'm in the learning stage of FOTRAN. So please don't mind if my program has got any silly mistakes.
Please help in executing the program.

program rungekutta


print *,' DIFFERENTIAL EQUATION WITH 1 VARIABLE OF ORDER 1'

print*,"enter the value of x0"
read*,x

print*,"enter the value of y0"
read*,y

print*,"enter the value of step size"
read*,h

print*,"solution for the equation is",rk(x0,y0,h)


call runge(m,f1,n,h,x1)

end

function rk(y) result(r)

!Example: y'=(2*y^2+5)

real y,r

r = (2*(y**2)+5)
end function rk(y)

Subroutine runge(m,f1,n,h,x1)

integer f1,i,j,m,n,h

real x,w1,w2,w3,w4,x0,y;

do i=1, m
ni = (i - 1) * f1 - 1
do j=1, f1
x = x1 + h * (n + j)
w1 = h * f(x0,y0)
w2 = h * f(x0+h/2, y0+w1/2)
w3 = h * f(x0+h/2, y0+w2/2)
w4 = h * f(x0+h, y0+w3)
x = x0 + h

x = y + (w1 + w2 + w2 + w3 + w3 + w4) / 6

end do

end do

end
There are at least several problems that you need to rectify.
(In future, include all the error messages with your post, along with the
source.)

1. You need "IMPLICIT NONE" in the program and in each subroutine
and function.

2. What is "f"?

3. when you invoke rk, it is with 3 arguments.
Yet the function is defined with only one dummy argument.

4. You requested the input of "x0", but you read the value into "x", etc.


Register to reply

Related Discussions
Anyone heard of Autonomous University of Madrid theoretical physics programm msc? Academic Guidance 0
Help on conjugate gradient programm in matlab Math & Science Software 0
Legendre poly in Fotran 90 Programming & Computer Science 2
Opening an fotran .exe file from Visual Basic Computing & Technology 3
Digital controller: how to programm it? Electrical Engineering 3