Initial value problems, math, Mathematics

Assignment Help:
Write a Matlab function MyIVP that solves an initial-value problem (IVP) for a system of ordinary
differential equations (ODEs) of the form
x ?(t) = f (t, x(t)),
where f : R × Rn ? Rn is an arbitrary function with one one-dimensional input (for time t) and one n-dimensional input x, and n-dimensional output. The function should implement a Runge-Kutta formula (for example, the rk34 formula or the Dormand & Prince formula).
The initial value x0 is provided by the user of MyIVP. The first line of MyIVP (saved in a file MyIVP.m) should look like this
function [xend,t,xt]=MyIVP(f,x0,tspan,N) Inputs
• f: function defining the right-hand side of the ODE. f should accept two arguments: t (a number) and x (an n-dimensional vector). The function f should return an n-dimensional vector y (the time derivative). Typical calling sequence: y=f(t,x), returning the value of f at time t in position x.
• x0: initial value where integration starts from (n-dimensional vector).
• tspan: Starting time and end time for integration. Integration has to run from time t =tspan(1)
to time t =tspan(2).
• N: number of steps for integration. The integration stepsize h=(tspan(2)-tspan(1))/N should
be small.
Outputs
• xend: result of integration at t =tspan(2).
• t: vector of times at which intermediate values have been computed (this should have N + 1
entries).
• xt: intermediate values (n × (N + 1)-array). xt(:,k) should be the solution at t(k).
You can check the built-in variable nargout inside your function to see if the user wants to get three outputs or only the end value xend. If nargout==1 you don’t need to store the intermediate values.

Related Discussions:- Initial value problems, math

Math project , Topic 1: Statistical Studies Find two different news storie...

Topic 1: Statistical Studies Find two different news stories in a mainstream media source (CNN, FoxNews, Newsweek, etc.), that cite data from a recognized poling agency. Locate th

The mean value theorem, The Mean Value Theorem : In this section we will ...

The Mean Value Theorem : In this section we will discuss the Mean Value Theorem.  Before we going through the Mean Value Theorem we have to cover the following theorem. Ro

Find the value of ((a+b)/(a-b)) , If arg (a/b) = pi/2, then find the value ...

If arg (a/b) = pi/2, then find the value of ((a+b)/(a-b)) where a,b are complex numbers. Ans) Arg (a/b) =Pi/2 Tan-1   (a/b)=   Pi/2 A/B = tanP/2 ,therefore a/b=infinity.

What is stem-and-leaf plots, Q. What is Stem-and-Leaf Plots? Ans. ...

Q. What is Stem-and-Leaf Plots? Ans. A stem-and-leaf plot is a table that provides a quick way to arrange a set of data and view its shape, or distribution. Each data val

Rounding, what is the result if 816.537 is rounded to the nearest tenth

what is the result if 816.537 is rounded to the nearest tenth

Write algorithm for the multiplication of a 3-digit number, E1) Why do we s...

E1) Why do we shift the place by one, of the result in the second row of the calculation, when we multiply, say, 35 by 237 E2) Write down the algorithm for the multiplication of

Shares and divend, a company of 10000 shares of rs 100 each declares a annu...

a company of 10000 shares of rs 100 each declares a annual dividend of 5 %.what is the total amount dividend paid by the company

Prove that seca+tana=2x, If secA= x+1/4x, prove that secA+tanA=2x or  1/2x....

If secA= x+1/4x, prove that secA+tanA=2x or  1/2x. Ans:    Sec? = x +  1/4x ⇒ Sec 2 ? =( x + 1/4x) 2                             (Sec 2 ?= 1 + Tan 2 ?) Tan 2 ? = ( x +

Calculate the area and circumference of a circle, Calculate the area and ci...

Calculate the area and circumference of a circle: Calculate the area and circumference of a circle with a 3" radius.  Solution: A =      πr2

Determine the area of the matting, A circular print is being matted in a sq...

A circular print is being matted in a square frame. If the frame is 18 in by 18 in, and the radius of the print is 7 in, what is the area of the matting? (π = 3.14) a. 477.86 in

Write Your Message!

Captcha
Free Assignment Quote

Assured A++ Grade

Get guaranteed satisfaction & time on delivery in every assignment order you paid with us! We ensure premium quality solution document along with free turntin report!

All rights reserved! Copyrights ©2019-2020 ExpertsMind IT Educational Pvt Ltd