1. Limited time only! Sign up for a free 30min personal tutor trial with Chegg Tutors
    Dismiss Notice
Dismiss Notice
Join Physics Forums Today!
The friendliest, high quality science and math community on the planet! Everyone who loves science is here!

Homework Help: MatLab Approximations

  1. Oct 10, 2012 #1
    1. The problem statement, all variables and given/known data
    I need to answer the attached question, but we haven't done anything similar in class and the book isn't proving to be helpful for me. I'm not sure if I have to use matlab, we were given this code but I can't even understand what it is calculating.

    % f(x), the function to integrate
    % f= @(x) x^4-2*x ;
    % f= @(x) exp(x);
    f=@(x) sin(x);
    % a, the lower limit of integration
    a=0 ;
    % b, the upper limit of integration
    b=pi ;
    % b=1.0;
    % n, the number of segments. Note that this number must be even.
    % n=20 ;
    %**********************************************************************
    format long g
    h=(b-a)/n ;
    % Sum the odd index function values, and multiply by 4
    sumOdd=0 ;
    for i=1:2:n-1
    sumOdd=sumOdd+f(a+i*h) ;
    end
    % Sum the even index function values, and multiply by 2
    sumEven=0 ;
    for i=2:2:n-2
    sumEven=sumEven+f(a+i*h) ;
    end
    sum=4*sumOdd+2*sumEven+f(a)+f(b) ;
    % Then multiply by h/3
    approx=h/3*sum ;
    %exact = quad(f,a,b) ;
    %exact=exp(b)-exp(a);
    exact=-cos(b)-(-cos(a));
    error=abs(approx-exact);
    disp(approx);
    disp(exact);
    disp(error);

    Any help would be greatly appreciated!
     

    Attached Files:

  2. jcsd
  3. Oct 10, 2012 #2

    Ray Vickson

    User Avatar
    Science Advisor
    Homework Helper

    It looks like a Simpson's rule approximation to the integral ##\int_0^{\pi} \sin(x) \, dx, ## using 20 equal intervals of length π/20. You are supposed also to compare the approximate and exact answers.

    Note: I am not a Matlab user and so am unfamiliar with all the syntax and commands, etc., but it is more-or-less obvious what this code is attempting.

    RGV
     
  4. Oct 11, 2012 #3
    Thank you for the help, but you don't need to make me feel bad for not understanding something you see is obvious. I do appreciate you pointing it out to me, I was having hard time seeing if it was the composite version.
     
Share this great discussion with others via Reddit, Google+, Twitter, or Facebook




Loading...
Similar Threads for MatLab Approximations Date
Rewriting ODE's into lower orders Nov 15, 2016
Matrix riccati differential equation using matlab Aug 5, 2016
Convolution in MATLAB Sep 28, 2015
Matlab help Jul 12, 2015
MATLAB Taylor Approximation Feb 21, 2011