# qd - %For loop to recursively compute the solution for i =...

This preview shows page 1. Sign up to view the full content.

%This function plots the solution to the linearly dispersive wave equation. First, save the file in a suitable folder and change the matlab environment to that particular folder. One can call the function by giving a suitable value for the time and position, note here that the position is an array of numbers. function [result] = qd(t,x) %Name of the function. dmn = size(x);%Size of the position array. x_d = dmn(2); result = zeros(size(x));%Providing the size of ‘result’. sum_t = zeros(size(x));%Providing the size of ‘sum_t’
This is the end of the preview. Sign up to access the rest of the document.

Unformatted text preview: %For loop to recursively compute the solution. for i = 1:x_d for j = 0:1000 sum_t(i) = sum_t(i) + (sin(((2*j+1)*x(i))-((2*j +1)^3)*t)/(2*j+1)); end end %For loop to compute the final solution. for k = 1:x_d result(k) = .5 - (2/pi)*sum_t(k); end plot(x,result); % Plot the solution xlabel('x ','FontSize',14); % Labeling the ‘x’ axis ylabel('u(t,x)','FontSize',14); % Labeling the ‘y’ axis title('Solution to the linearly dispersive wave equation','FontSize',14);% Providing the title to the graph....
View Full Document

## This note was uploaded on 11/19/2011 for the course MATH 101 taught by Professor Wormer during the Spring '08 term at UCSD.

Ask a homework question - tutors are online