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

Find out the volume of the solid- method of rings, Find out the volume of t...

Find out the volume of the solid obtained by rotating the region bounded by y = x 2 - 2x and  y = x about the line y = 4 . Solution: Firstly let's get the bounding region & t

Determine the angle, In parallelogram ABCD, m∠A = 3x + 10 and m∠D = 2x + 30...

In parallelogram ABCD, m∠A = 3x + 10 and m∠D = 2x + 30, Determine the m∠A. a. 70° b. 40° c. 86° d. 94° d. Adjacent angles in a parallelogram are supplementary. ∠A a

Calculus, I need help with my calculus

I need help with my calculus

Trignometry: sin-3x, sin(2x+x)=sin2x.cosx+cos2x.sinx              =2sinxco...

sin(2x+x)=sin2x.cosx+cos2x.sinx              =2sinxcosx.cosx+(-2sin^2x)sinx              =2sinxcos^2+sinx-2sin^3x             =sinx(2cos^2x+1)-2sin^3x =sinx(2-2sin^2x+1)-2sin^3

Drawn to a circle with center o, From a point P, two tangents PA are drawn ...

From a point P, two tangents PA are drawn to a circle with center O.If OP=diameter of the circle show that triangle APB is equilateral. Ans:    PA=PB (length of tangents

Constructing tables versus rote learning maths, CONSTRUCTING TABLES VERSUS ...

CONSTRUCTING TABLES VERSUS ROTE LEARNING :  Ask any adult how she would help a child to acquire simple multiplication facts. There is a very strong possibility that she would say,

Topological spease, let X be a nonempty set. let x belong to X. show that t...

let X be a nonempty set. let x belong to X. show that the collection l={ union subset of X : union = empty or belong U

Maths for fun-mathematics- in our lives, Maths For Fun :  Often, when I ha...

Maths For Fun :  Often, when I have time on my hands, I try to solve interesting mathematical questions of the following kind. Sometimes my friends and I create the problems, and

Calculate the net amount and distance, 1. A train on the Bay Area Rapid Tra...

1. A train on the Bay Area Rapid Transit system has the ability to accelerate to 80 miles/hour in half a minute. A.   Express the acceleration in miles per hour per minute. B

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