| draw n-th eigenfunction (n=4 ... the lowest positive eigenvalue) |
>> N=160
>> mu=0.3
>> a=plate_f1(N,mu);
>> [v,d]=eigs(a,200,0);
>> format long g
>> sort(diag(d))
ans =
2.62483490092012e-07
3.23730085560577e-07
3.97631977037723e-07
181.367464056553
383.879231655339
588.857298220385
1210.6936945447
1210.69369457591
3729.71656195089
3729.716561981
(後略)
>> n=4
>> plot_n(v(:,201-n),N,N)
|