Please fix my code
3 views (last 30 days)
Show older comments
to solve y'=t(y.^2-4)
I try like this.
but something was wrong... Amend pls......
clear all
[x,y]=meshgrid(linspace(0,3,100));
u=ones(100);
v=(y.^2-4)*(1/2).*(log(abs(2-y/2+y)).^(1/2));
quiver(x,y,u,v,0.8)
hold on
t=0, w=1;
dt=0.01;
N=1/dt;
for n=1:N
t=t+dt;
w=w+dt*(y.^2-4)*(1/2).*(log(abs(2-y/2+y)).^(1/2));
st=[st t
end
plot(st,sw,'r')
what is wrong?
4 Comments
DGM
on 17 May 2021
What's the error message say?
Error using horzcat
Dimensions of arrays being concatenated are not consistent.
sw=[sw w];
so look at the size of sw and w
>> size(sw)
ans =
1 1
>> size(w)
ans =
100 100
why is sw the size it is? It's initialized as 1. You could initialize it as an empty vector [], but pay attention to the size of w.
w=1;
% ...
sw=w;
why is w the size it is?
[x,y]=meshgrid(linspace(0,3,100));
% ...
w=w+dt*(y.^2-4)*(1/2).*(log(abs(2-y/2+y)).^(1/2));
w is 2D. You're going to figure out how you want to plot these. Once you figure that out, you'll have to figure out how you need to store them. If you need to store all instances of them, you can concatenate them on dim3.
Answers (0)
See Also
Categories
Find more on Axis Labels in Help Center and File Exchange
Community Treasure Hunt
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!