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

Complex eigenvalues, It is the last case that we need to take a look at. Th...

It is the last case that we need to take a look at. Throughout this section we will look at solutions to the system, x?' = A x? Here the eigenvalues of the matrix A are compl

What does required to earn on his further science test in 93, Justin earned...

Justin earned scores of 85, 92, and 95 on his science tests. What does he required to earn on his further science test to have an average (arithmetic mean) of 93%? To earn an a

Green function, greens function for x''''=0, x(1)=0, x''(0)+x''(1)=0 is G(t...

greens function for x''''=0, x(1)=0, x''(0)+x''(1)=0 is G(t,s)= {1-s for t or equal to s

Find integer if consecutive even integers is the number 126, The sum of two...

The sum of two consecutive even integers is the number 126. What are the integers? Two consecutive even integers are numbers in sequence, such as 4 and 6 or -30 and -32, that a

Measures of central tendency-graphical method , Illustration In a soci...

Illustration In a social survey whether the main reason was to establish the intelligence quotient or IQ of resident in a provided area, the given results were acquired as tab

Phase plane, Before proceeding along with in fact solving systems of differ...

Before proceeding along with in fact solving systems of differential equations there's one topic which we require to take a look at. It is a topic that's not at all times taught in

Comparison-types of word problems related to subtraction, Comparison - the...

Comparison - the difference between two groups or numbers, namely, how much one is greater than the other, how much more is in one group than in the other. (e.g., if Munna has

How to solve systems of equations, How to solve Systems of Equations ? ...

How to solve Systems of Equations ? There's a simple method that you can use to solve most of the systems of equations you'll encounter in Calculus. It's called the "substitut

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