Not enough input arguments
Afficher commentaires plus anciens
Hi there, I was running the spiral.m script in my computer and getting the following error.
Not enough input arguments.
"Error in spiral (line 19)
h = 2*L/n;"
My code:
function [T,W] = Spiral(n,L,Tend,nplots,D,e)
% function [T,Y] = Spiral(n,L,Tend,nplots,D,e)
% nominal usage: [T,Y] = Spiral(128,8,50,20,0.1,0.1)
% method of lines finite difference approximation for
% Fitzhugh-Nagumo equations, periodic domain
% Input: n is number of subintervals across domain
% L is half-width of domain, -L<x,y<L
% Tend is the end time, 0<t<Tend
% nplots is the number during the time interval (nplots+1 if t=0
% included)
% D is diffusivity (try 1,0.5,0.25,0.1)
% Output: T is time levels for answer (nplots+1)
% Y is the answer at T, needs to be reshaped to n+1 by n+1 grid
b = 0.67;
%e = 1;
%D = 0.1;
h = 2*L/n;
x = -L + h*(0:n)'; % includes boundary
y = x;
[X,Y] = ndgrid(x,y);
Xi = X(1:n,1:n); Yi = Y(1:n,1:n); % interior points + one boundary
% set up ICS
u = 5*ones(n,1)*([1:n]-0.4*n)/n;
v = 5*(([1:n]-0.4*n)')*ones(1,n)/n;
y0(1:n^2,1) = reshape(u,n^2,1);
y0(n^2+1:2*n^2,1) = reshape(v,n^2,1);
w=y0;
% Set up the one-dimensional 2nd-derivative matrices.
Dxx = toeplitz( [-2 1 zeros(1,n-2)]/h^2 );
Dxx(1,n) = 1; Dxx(n,1) = 1; % periodic box case
Dxx = sparse(Dxx);
Dyy = Dxx;
% Form discrete Laplacian.
I = speye(n);
A = sparse(kron(Dxx,I) + kron(I,Dyy));
MyTols = odeset('RelTol',1e-4,'AbsTol',1e-5);
[T,W] = ode45(@(t,u) rhs(t,u,A,D,b,e,n),[0:Tend/20:Tend],y0,MyTols);
for ii = 1:length(T)
u = reshape(W(ii,1:n^2),n,n);
max(max(u))
v = reshape(W(ii,n^2+1:2*n^2),n,n);
figure
% colormap(hot); axis image off;
%surface(Xi,Yi,u)
colormap(hsv), pcolor(u)
pause(1)
end
function dydt = rhs(t,z,A,D,b,e,n)
% right-hand side.
u1 = z(1:n^2,1);
u2 = z(n^2+1:2*n^2,1);
dydt = z;
dydt(1:n^2,1) = D*(A*u1)+(u1-u1.^3/3-u2)/e;
dydt(n^2+1:2*n^2,1) = e*(u1-u2/2+b);
end
end
What should I do now?
Best Regards
Zillur
1 commentaire
Tamir Suliman
le 4 Déc 2016
please upload the code or fix the formatting I dont know which is commented out and which is not
Réponse acceptée
Plus de réponses (0)
Catégories
En savoir plus sur Error Functions dans Centre d'aide et File Exchange
Community Treasure Hunt
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!