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

The parallelogram, love is a parallelogram where prove that is a rectangle...

love is a parallelogram where prove that is a rectangle

Wavy curve method, In order to compute the inequalities of the form ...

In order to compute the inequalities of the form   where n 1 , n 2 , ....... , n k , m 1 , m 2 , ....... , m p are natural and real numbers and a 1 , a 2 , ... , a k ,

Triangle treat, what letters to fill in the boxes

what letters to fill in the boxes

Find the time it takes to slide down, The time in seconds that it takes for...

The time in seconds that it takes for a sled to slide down a hillside inclined at and angle θ is where d is the length of the slope in metres. Find the time it takes to sli

How to solve lim 1-cos(x)/1-cos(4x) as x tends to zero, Use L''hopital''s r...

Use L''hopital''s rule  since lim X-->0  1-cos(x)/1-cos(4x)  is in the indeterminate form 0/0 when we apply the limt so by l''hoptital''s rule differentiate the numerator and den

Queuing Theory, A telephone exchange has two long distance operators.The te...

A telephone exchange has two long distance operators.The telephone company find that during the peak load,long distance calls arrive in a poisson fashion at an average rate of 15 p

Numerical analysis and computer techniques, write a fortan programme to gen...

write a fortan programme to generate prime number between 1 to 100

Compare and contrast african immigrants, Compare and contrast African immig...

Compare and contrast African immigrants with our immigrant groups? How are they different? What are the implications of these differences for their adjustment to the larger society

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