Solution Codes of Ordinary Diferenatial Equation
조회 수: 1 (최근 30일)
이전 댓글 표시
I wrote MATLAB codes for
equation with parameter n=450, a= -1.99 and x_0=3/4 as follows.
clc; clear; clear all;
a=-1.99; x0=3/4;
sumxx=[];
k=450;
for p=1:k
T=linspace(p-1,p,101);
sumx=1+a.*T;
for n=0:p-2
sumx=sumx+(a.^2).*(((1+a).^n).*(T-(n+1)));
end
sumx=sumx.*x0;
sumxx=[sumxx,sumx];
clear sumx
end
tt=linspace(0,k,k*101);
plot(tt,sumxx)
grid on
hold on
But with the same parameters and in the equation above at two places division by 4, i couldn't arrange MATLAB codes. Please, help me.
댓글 수: 2
답변 (1개)
Jan
2022년 5월 21일
편집: Jan
2022년 5월 21일
What is the purpose of the outer loop "for p=1:k"? Why does the inner loop "for n=0:p-2" run from 0 to p-2?
It is confusing that you've swapped n and k.
The purpose of your code is not clear to me. I have no idea, why inserting /4 twice is a problem and why you assume, that the graph is "not correct". I cannot guess what you expect instead.
n = 450;
a = -1.99;
x0 = 0.75;
T = linspace(0, 1000, 1000);
X = 1 + a * T;
for k = 0:n
X = X + a^2 * (1 + a / 4)^k * (T - (k + 1) / 4);
end
X = X * x0;
figure;
plot(T, X)
댓글 수: 0
참고 항목
Community Treasure Hunt
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!