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

Algebra 1, Im having trouble with this word problem: The three Math Idol j...

Im having trouble with this word problem: The three Math Idol judges have been eliminating contestants all day! The number of one-step equations and two-step equations who have be

Find the constant height at which the jet is flying, The angle of ...

The angle of elevation of a jet fighter from a point A on the ground is 600. After a flight of 15 seconds, the angle of elevation changes to 300. If the jet is flying at a speed  o

Derivatives, application of derivatives in engg.

application of derivatives in engg.

Project, report on shares and dovidend using newspaer

report on shares and dovidend using newspaer

Circles, Two tangents TP and TQ are drawn to a circle with center O from an...

Two tangents TP and TQ are drawn to a circle with center O from an external point T.prove that angle PTQ=angle 2 OPQ

Rules of logarithms, Rule 1 The logarithm of 1 to any base is 0. Pro...

Rule 1 The logarithm of 1 to any base is 0. Proof We know that any number raised to zero equals 1. That is, a 0 = 1, where "a" takes any value. Therefore, the loga

Velocity of derivation, Velocity : Recall that it can be thought of as sp...

Velocity : Recall that it can be thought of as special case of the rate of change interpretation. If the situation of an object is specified by f(t ) after t units of time the vel

The quantity x + 6 is divided by negative four find number, Negative four i...

Negative four is multiplied through the quantity x + 8. If 6x is then added to this, the output is 2x + 32. What is the value of x? twice the quantity x + 6 is divided by negative

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