Info
此问题已关闭。 请重新打开它进行编辑或回答。
can anyone help me ?
1 次查看(过去 30 天)
显示 更早的评论
here is a code matlab of one degree differential equation, i want to plot the error method for different values of N, can anyone help me
2 个评论
Geoff Hayes
2019-6-2
diadalina - please clarify what you mean by i want to plot the error method for different values of N. Does your code do some or all of this already? If not, where should this happen in your attached code?
回答(1 个)
Stephan
2019-6-2
编辑:Stephan
2019-6-2
clear all;
clc;
close all;
N=[100,200,300,400];
hold on
for k=1:numel(N)
%%les données mathématiques
t = 0:1/(N(k)-1):1;
%les conditions initiales
c0=0;
%%résolution de notre système
[~,Y]=ode45(@systest,t,c0 );
sol_exacte=t.*exp(t);
subplot(2,1,1)
plot(t,Y(:,1),t,sol_exacte,'+')
legend('sol-app','sol-exacte')
err(k)=norm(Y(:,1)'-sol_exacte)
e=err(k);
end
hold off
subplot(2,1,2)
semilogy(N,err,'ob--')
legend('error')
function F=systest(t,Y)
F(1,1)=Y(1)+exp(t);
end
0 个评论
此问题已关闭。
另请参阅
Community Treasure Hunt
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!