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 final position of point by rotation -translation matric, Question: ...

Question: A point in 3D is first rotated anticlockwise by 45 degrees about x axis,then translated along y axis by 2 units.Find the final position of the point if its initial po

Number line, I need to graph rational numbers on the number line Point A-.6...

I need to graph rational numbers on the number line Point A-.60, point B-1/4, point C-.4,point D-7/8

Technique of teching, What is a review technique? What are its advantages a...

What is a review technique? What are its advantages and disadvantages?

Mathematics is all around us-mathematics- in our lives, Mathematics Is All ...

Mathematics Is All Around Us :  What is the first thing you do when you get up? Make yourself a nice cup of tea or coffee? If so, then you're using mathematics! Do you agree? Cons

Example of word problems involving money, Example of Word Problems Involvin...

Example of Word Problems Involving Money: A collection of coins consists of nickels, dimes & quarters. The number of quarters is double the number of nickels, and the number o

Calculate the fourier cosine series, The Fourier series expansion for the p...

The Fourier series expansion for the periodic function, f ( t ) = |sin  t | is defined in its fundamental interval. Taking π = 3.142, calculate the Fourier cosine series app

Addition involving negative numbers, Q. Addition Involving Negative Numbers...

Q. Addition Involving Negative Numbers? Ans. When you add together positive and negative numbers, there are essentially three possibilities that you can encounter. Let's e

ALGEBRA, FIND PRODUCT (-41)*(102)

FIND PRODUCT (-41)*(102)

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