Fix issue moiseevigor#23 nonformant arguments in ellipj

I made no attempt to figure out the origin of this bug and so there may be a "cleaner" fix. However, this fix does the job.
cffk · Oct 18, 2023 · bd8a3b4 · bd8a3b4
1 parent 2b1bb98
commit bd8a3b4
Showing 1 changed file with 1 addition and 1 deletion.
diff --git a/ellipj.m b/ellipj.m
@@ -92,7 +92,7 @@
 
  mmax = length(I);
  phin = zeros(1,mmax);
- phin(:) = (2 .^ double(n(K))).*a(i,K).*u(I);
+ phin(:) = (2 .^ double(n(K)')).*a(i,K).*u(I);
  while i > 1
  i = i - 1;
  in = uint32( find(n(K) >= i) );