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 the return period in years of the critical event, 1. A drainage system...

1. A drainage system overflows into a river causing pollution when the rainfall falling in an hour exceeds 50 mm. During any given year, the maximum hourly rainfall follows an expo

Find out a vector that is orthogonal to the plane, A plane is illustrated b...

A plane is illustrated by any three points that are in the plane.  If a plane consists of the points P = (1, 0,0) , Q = (1,1,1) and R = (2, -1, 3) find out a vector that is orthogo

SHARES AND DIVIDEND, i am a student of class 10 and need help for making my...

i am a student of class 10 and need help for making my project on shares and dividend

Algebra, prove That J[i] is an euclidean ring

prove That J[i] is an euclidean ring

Find out indegree, Question: Consider a digraph D on 5 nodes, named x0...

Question: Consider a digraph D on 5 nodes, named x0, x1,.., x4, such that its adjacency matrix contains 1's in all the elements above the diagonal A[0,0], A[1,1], A[2,2],.., e

Tower of hanoi, how to create an activity of tower of hanoi

how to create an activity of tower of hanoi

Matrix, cramers rule introduction

cramers rule introduction

Invertible Matrices - Logical equivalence, 10 statements must be shown to b...

10 statements must be shown to be logically equivalent to the Statement the nxn matrix is invertible.

Probability, I have a question that hurts my head to work out. It is really...

I have a question that hurts my head to work out. It is really confusing for me. It sais " By the start of the 21st century, only 1 in 6 babies in America was born with blue eyes.

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