Mathematics Stack Exchange is a question and answer site for people studying math at any level and professionals in related fields. Join them; it only takes a minute:

Sign up
Here's how it works:
  1. Anybody can ask a question
  2. Anybody can answer
  3. The best answers are voted up and rise to the top

I was looking speedy and simple algorithm to solve second order differential equation in numarical way. I thought an idea to solve it via geometrical equalities.

See my approach to the problem $$A''(x)=f(x) A(x)$$ I will try to solve $A'(1)$ via numerical method. Initial condition $A'(0)=y(0)=k$ and $A(0)=0$

$$A(x)=\int _0^x y(x) dx$$

$$A''(x)=y'(x) $$

$$y'(x)=f(x) A(x)$$

$A(x)$ is total area of $y(x)$ from $0$ to $x$

$y'(x)$ is $tan(\alpha)$ on $x$

If we divide n point between $0$ and $1$

for x=1/n



if $y'(x)=f(x)A(x)$ then $\frac{y(\frac{1}{n})-y(0))}{(1/n)}=y(0)f(\frac{1}{n})\frac{1}{n}$


for x=2/n



if $y'(x)=f(x)A(x)$ then $\frac{y(\frac{2}{n})-y(\frac{1}{n}))}{(1/n)}=f(\frac{2}{n})y(0)\frac{1}{n}+f(\frac{2}{n})y(\frac{1}{n})\frac{1}{n}$


for x=3/n

$y(3/n)=y(\frac{2}{n})[\frac{1}{n^2}f(\frac{3}{n})+1]+f(\frac{3}{n})y(\frac{1}{n})\frac{1}{n^2}+f(\frac{3}{n})y(0)\frac{1}{n^2}$ .





for x=n/n=1


if we increase n, I think we will be near to value.

Is my approach correct or not? I don't have matlab or other tool to simulate how successfull the idea to solve the equation. Could you please also share your experiences and theories about numerical methods to solve such second order differantial equations?

Thanks in advice and times for me and your shares.

share|cite|improve this question
Please exercise restraint while using $$...$$. (I have now edited it!) – user21436 Mar 3 '12 at 22:10

Your Answer


By posting your answer, you agree to the privacy policy and terms of service.

Browse other questions tagged or ask your own question.