Matlab code problem with differential equations

[pugtm]

Homework Statement

For a following differential equation
d^2y/dx^2-4y=(e^x)/x
 Find the solution using numerical methods

Homework Equations

d^2y/dx^2-4y=(e^x)/x

The Attempt at a Solution

%num
dx=0.01;
x=1:dx:3;
l=zeros(1,length(x));
m=zeros(1,length(x));
l(1)=1;
m(1)=0.25;
for ii=2:length(x)
    l(ii) = l(ii-1)+m(ii-1)*dx;
    m(ii) = m(ii-1) -(1/4) * (l(ii-1)) * dx +(exp(2*x*(ii-1)/(x(ii-1))))*dx;<----- this line keeps throwing an error
end
plot(x,l);

any assistance would be greatly appreciated

 

[jedishrfu]

What error are you getting?

If its during execution then perhaps a print statement inside your for loop printing out the results of L(ii) and M(ii) and some of your expressions will give you some insight into what going on.

You should be able to use the computer to help you find the error in your code by tracing its operations.

 

[Delta2]

The way the code is it adds at every step $e^{\frac{2x}{x}}dx$ while I believe you want to add $\frac{e^{2x}}{x}dx$. Check carefully the parentheses in the (exp(2*x)...) term at that line.

Also what it seems to be a typo it says exp(2*x*(ii-1)... while I believe it should be exp(2*x(ii-1))...

Also not sure since I am rusty on numerical methods and MATLAB code but this looks like the Euler Method for solving ODEs correct?

 

[jedishrfu]

Good catch that x*(ii-1) is clearly wrong as the program needs x(ii-1)