Dismiss Notice
Join Physics Forums Today!
The friendliest, high quality science and math community on the planet! Everyone who loves science is here!

Euler's Method in MAtlab

  1. Apr 8, 2012 #1
    Hello
    I have a program for Eulers method >>

    % Euler's Method for dy/dt = cost
    k = 1;
    y0 = 0;


    npoints = 500;
    dt = 0.01;

    y = zeros(npoints,1); % this initializes the vector y to being all zeros
    t = zeros(npoints,1);

    y(1) = y0; % the initial condition
    t(1) = 0.0;

    for step=1:npoints-1 % loop over the timesteps
    y(step+1) = y(step) + dt*k*(cos(y(step)));
    t(step+1) = t(step) + dt;
    end

    plot(t,y,'r'); %plots the numerical solution in red
    hold on; %keep the previously plotted lines
    plot(t,yexact2(t)); %plots the exact solution (default plot is in blue, solid line)

    The graph asymptotically approaches pi/2, and I cant think of a way to deal with this.
     
  2. jcsd
  3. Apr 8, 2012 #2

    Office_Shredder

    User Avatar
    Staff Emeritus
    Science Advisor
    Gold Member

    It sounds like you're getting the correct graph
     
Know someone interested in this topic? Share this thread via Reddit, Google+, Twitter, or Facebook




Similar Discussions: Euler's Method in MAtlab
Loading...