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

Angles, samuel left mauritius at 22:30 on saturday and travelled to london ...

samuel left mauritius at 22:30 on saturday and travelled to london (GMT) for 14h30min he had a stopover for 4 h in london and he continued to travel to toronto for another 6h20min

Matrix, find the matrix of the linear transformations T:R2->R2 defined by T...

find the matrix of the linear transformations T:R2->R2 defined by T(x,y,z)=(x+2y,x-3z).

Proof of various derivative facts formulas properties, PROOF OF VARIOUS DER...

PROOF OF VARIOUS DERIVATIVE FACTS/FORMULAS/PROPERTIES Under this section we are going to prove several of the different derivative facts, formulas or/and properties which we en

Applied mathematics, I have a journal article in applied mathematics and wa...

I have a journal article in applied mathematics and want to analyze the solutions step by step. Is there anyone specialize in this file?

Geometry, triangular with base AB = 48cm and height CH=16cm is inscribed a ...

triangular with base AB = 48cm and height CH=16cm is inscribed a rectangle MNPQ in which MN: MQ = 9:5 Find MN and MQ

Potency of a drug , An experiment designed to test the potency of a drug on...

An experiment designed to test the potency of a drug on 20 rats. Last animal studies have shown that a 10 mg dose of the drug is lethal 5% of the time within the first 4 hours; of

Direction cosines - vector, Direction Cosines This application of the ...

Direction Cosines This application of the dot product needs that we be in three dimensional (3D) space not like all the other applications we have looked at to this point.

Compute standard and variance deviation, A firm is manufacturing 45,000 uni...

A firm is manufacturing 45,000 units of nuts. The probability of having a defective nut is 0.15 Compute the given i. The expected no. of defective nuts ii. The standard an

Example of partial fraction decomposition, Example of Partial Fraction Deco...

Example of Partial Fraction Decomposition Evaluate the following integral. ∫ (3x+11 / x 2 -x-6) (dx) Solution: The 1 st step is to factor the denominator so far as

Calculate the score of contestant on a tv game show, Danny is a contestant ...

Danny is a contestant on a TV game show. If he gets a question right, the points for that question are added to his score. If he gets a question wrong, the points for that question

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