MathGroup Archive 2007

[Date Index] [Thread Index] [Author Index]

Search the Archive

Write a computer program to using Simpson's rule (Have to revise)

  • To: mathgroup at smc.vnet.net
  • Subject: [mg74121] Write a computer program to using Simpson's rule (Have to revise)
  • From: "Evanescence" <origine26 at yahoo.com.tw>
  • Date: Mon, 12 Mar 2007 04:32:29 -0500 (EST)

Hello Dear all:
My questions are as follows:
First I definite as follows functions.
B[c_, m_] function of c, m
A[c_, m_] function of c, m
F[c_] function of c
G[c_] function of c
K[n_, c_] function of n, c
Z[c_] function of c
U[c_] function of c

Then I want to write a computer program about Simpson's 1/3 rule that
can express below function
L[m_, n_] :=(1/3)*((16-0)/32)*(0
+4(((1/2)*K[n, 1/2] )*(B[1/2, m]*F[1/2]*(Z[1/2])+2*( A[1/2, m]-B[1/2,
m]*G[1/2])*U[1/2]))
+2 (((1)*K[n, 1] )*(B[1, m]*F[1]*(Z[1])+2*( A[1, m]-
B[1,m]*G[1])*U[1]))
+4(((3/2)*K[n, 3/2] )*(B[3/2, m]*F[3/2]*(Z[3/2])+2*( A[3/2, m]-B[3/2,
m]*G[3/2])*U[3/2]))
+2(((2)*K[n, 2] )*(B[2, m]*F[2]*(Z[2])+2*( A[2, m]-B[2,
m]*G[2])*U[2]))
+4.......................................................................
+(((16)*K[n, 16] )*(B[16, m]*F[16]*(Z[16])+2*( A[16, m]-B[16,
m]*G[16])*U[16])))


the 16 is upper limit of integration
and the 0 is lower limit of integration
subdivision is 1/2


I want to write a computer program using the Simpson's rule that can
choose the subdivision by myself(maybe can choose 0.001 or 0.002
etc...)and upper and lower limit can choose by myself too,final can
express the above function.


Thank you

Evanescence 2007,3,12



  • Prev by Date: Re: Factorizing...
  • Next by Date: Re: accessing a kernel on a network
  • Previous by thread: Re: ComplexityFunction
  • Next by thread: PN junction Simulation with Mathematica