Computational Physics Programming using fortran g95

[jhosamelly]

Homework Statement

Create a program that computes for the value of sine function. Then compute for its
integral from 0 to pi with N intervals, where N=4,8,16,256 and 1024 and compare the
result for the trapezoid and simpson method.

Homework Equations

Trapezoid rule of Integration is (b-a)((f(a)+(f(b))/2)

The Attempt at a Solution

here is the code I made

but I'm getting "Unclassifiable statement" at the i1 equation.
Is there something wrong with my code? We are using fortran 90

I'm trying the trapezoid rule 1st.

 

[SteamKing]

Drop the '*' in your print statements for the unit specifier
The format part of the print statement should be enclosed in single quotes ['] instead of double quotes ["].

Thus, the first output should be:
print 'SIN(0) = ', sina

 

[SteamKing]

OOPS: forgot format descriptor. First output should be:

print '("SIN(0) = ", F6.2)', sina

 

[jhosamelly]

its still not working.. and this is how our prof taught us. we use
print *, "ddd" , ddddd

something wrong with my equation I think.

 

[Mark44]

Your i1 "equation" is not an equation - it's an assignment statement. An equation expresses the equality of two expressions. That's not what = is used for in programming. Instead, = is used to assign the value of the expression on the right side to the variable on the left side.

Your compile error "unclassifiable statement" comes from not having a multiplication operator. This assignment statement should be

i1 = (b - a) [B]*[/B][/color]((sina + sinb)/2)

You could also write this as

i1 = (b - a) * (sin(angrada) + sin(angradb))/2

In mathematics you can write a product of two numbers by putting them next to each other, with $2b$ meaning 2 times b.

There's no such shortcut in programming, though. If you need to perform a multiplication, you MUST use *.

 

[jhosamelly]

Thank you so much ))

 

[jhosamelly]

Hmmmm.. I can't seem to finish this program I'm doing. :( What should I do with N?

 

[Mark44]

You are a very long way from being done with this program, I'm afraid.

Create a program that computes for the value of sine function. Then compute for its integral from 0 to pi with N intervals, where N=4,8,16,256 and 1024 and compare the result for the trapezoid and simpson method.

The program will need to compute the value of sine(x) at a number of points in the interval [0, $\pi$]. One part will need to split this interval into four subintervals, and calculate the integral using the trapezoid method and Simpon's method. The next part will need to do the same thing using eight subintervals.

A third part will need to do the same thing, using 16 subintervals. A fourth part will need to do the same thing using 256 subintervals. Finally, the last part will need to do the same thing using 1024 subintervals.