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

Definition of differential equation, The first definition which we must cov...

The first definition which we must cover is that of differential equation. A differential equation is any equation that comprises derivatives, either partial derivatives or ordinar

Probability, The probability that a leap year will have 53 sunday is ? and ...

The probability that a leap year will have 53 sunday is ? and how please explain it ? (a)1/7    (b) 2/7    (c) 5/7    (d)6/7 Sol) A leap year has 366 days, therefore 52 weeks i.e

Show that the height of the aero plane, From  an  aero  plane  vertically  ...

From  an  aero  plane  vertically  above  a  straight  horizontal  road,  the  angles  of depression of two consecutive milestones on opposite sides of the aero plane are observed

Polynomial time algorithm - first order query, For queries Q 1 and Q 2 , w...

For queries Q 1 and Q 2 , we say Q 1 is contained in Q 2 , denoted Q 1 ⊆ Q 2 , iff Q 1 (D) ⊆ Q 2 (D) for every database D. The container problem for a fixed Query Q 0 i

Conclude the values of the six trigonometric functions, Conclude the values...

Conclude the values of the six trigonometric functions: Conclude the values of the six trigonometric functions of an angle formed through the x-axis and a line connecting the

Find the average, The center of a national park is located at (0,0). A spec...

The center of a national park is located at (0,0). A special nature preserve is bounded by by straight lines connecting the points A at (3,2), B at (5,1), C at (8,4) and D at (6,5)

Determine how many valid fortran identifiers, A valid identifier in the pro...

A valid identifier in the programming language FORTAN contains a string of one to six alphanumeric characters (the 36 characters A, B,...., Z, 0, 1,...9) starting with a letter. De

Wants to Join as expert, Hi.. This is dinesh kumar I just joined experminds...

Hi.. This is dinesh kumar I just joined experminds.com , i wamt to receive assignment in maths and want to complete students assignment within time. Please help me how i can become

Find out primes of each denominator, Find out primes of each denominator: ...

Find out primes of each denominator: Add 1/15 and 7/10 Solution: Step 1:             Find out primes of each denominator. 15 = 5 x 3 10 = 5 x 2 Step 2:

Cross product - vector, Cross Product In this last section we will loo...

Cross Product In this last section we will look at the cross product of two vectors.  We must note that the cross product needs both of the vectors to be three dimensional (3D

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